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This publication gives details of the 
computational subroutines available in the 
PL/I Library. These subroutines are used 
by the PL/I (F) compiler in the implementa- 
tion of PL/I built-in functions and of the 
operators used in the evaluation of PL/I 
expressions. Not all PL/I built-in func- 
tions and expression operators are support- 
ed by the PL/I Library; the compiler gener- 
ates in-line code for a small number of 
them. The details provided include timing 
figures, summaries of the mathematical 
methods used, and (where appropriate) fig- 
ures for range and accuracy. This informa- 
tion is intended to be of interest chiefly 
to those programmers concerned with the 
performance of computational subprograms. 



PREFACE 



This publication provides the PL/I pro- 
grammer with detailed information about the 
computational subroutines which are part of 
the OS/360 PL/I Library. 

The reader is assumed to be a programmer 
with a particular concern for performance 
information associated with individual 
modules. The numerical analyst is provided 
with a description of the algorithms, and a 
specification of accuracy and range, where 
these are considered to be significant. 

Useful background reading is provided in 
the following IBM publications: 



IBM System/ 360 Prin c iples of Operation , 
Form A2 2-6 821 



IBM System/360 Oper a ting System PL/I; 
Language Specificatio ns, Form C28-6571 



IBM System/3 60 Operat i ng System; PL/I (F ) 
Programmer's Guide , Form C28-659U 



IBM System/ 360 Oper a ting System; Assem- 
bler Language , Form C28-6514 



Copies of this and other IBM publications can be obtained through IBM 
Branch Offices. 
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It may be mailed directly to IBM. Address any additional comments 
concerning this publication to the IBM Corporation, Department D39,1271 
Avenue of the Americas, New York, N.Y. 10020. 
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COMPUTATIONAL SUBROUTINES 



INTRODUCTION 

The PL/I Library computational subrou- 
tines provide support for the operators and 
functions of the PL/I language in four 
major categories: 

1. Bit and Character Strings 

2. Arithmetic 

3. Mathematical 

4 . Arrays 

These subroutines have been designed to 
allow their use in a multi-tasking environ- 
ment. 

This publication gives detailed informa- 
tion in each of the four sections mentioned 



above with respect to performance, accura- 
cy, choice of algorithm, and range of 
values handled (where appropriate) . 



A number of exceptional conditions may 
arise in the execution of the library 
subroutines. Many of these are not direct- 
ly related to PL/I ON conditions. The 
method of treatment in these cases is to 
write a diagnostic message and raise the 
ERROR condition. This allows the user the 
opportunity of investigating the error by 
use of the ONCODE built-in function in his 
ON ERROR unit and of making a choice on the 
action which he wishes to take. Full 
details of the diagnostic messages printed 
at object time and of the ONCODE values 
associated with them are specified in IBM 
System/360 Operating System; PL/I (F) 
Programmer's Guide , Form C28-6594. 



CHAPTER 1: BIT AND CHARACTER STRING OPERATIONS AND FUNCTIONS 



The Library string package contains 
modules for handling bit and character 
strings. Generally, a string function or 
operator is supported by only one module, 
but in the interests of efficiency some of 
the bit string operators are provided with 
additional modules to deal with byte- 
aligned input data. 

Execution times are based on information 
in IBM System/360 Instruction Timing 



Information , Form A22-6825. They are 
intended to be simple enough to give a good 
general guide to performance while averag- 
ing out many logical variations, the effect 
of which is comparatively small. 



A complete list of the modules provided 
in the Library string package is given in 
Figure 1. 
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Figure 1. Bit and Character String Operations and Functions 



BIT STRING OPERATIONS 

The 'And' Operator (&) 

Module Name: IHEBSA 
Entry Point: IHEBSAO 

Function: 

To implement the 'and' operator between 
two byte-aligned bit strings, placing the 
result in a byte-aligned target field. 



Method: 

The current length of the target string 
is set equal either to the maximum of 
those of the operands,, or to the maximum 
length of the target field (when trunca- 
tion is necessary to avoid exceeding the 
length of this field) . The strings are 
' and'ed together for a length equal to 
the minimum of the lengths of the oper- 
ands, and the result is extended with 
zeros, if necessary, up to the current 
length calculated for the target field. 



Implementation : 

• Module size: 296 bytes 

• Execution time: 

Let L A = the length of the shorter 
string 

L 2 = the difference between the 
string lengths 

Mj_ = FLOORULi - D/2048) 

Ni = MODCLjl, 2048) - 1 

Pj_ = FLOOR (N-j/8) 

S = SIGN(L 2 ) 

i = 1,2 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*Mi + c*Pi + S*(d + e*M 2 + f*P 2 ) 



I I 



30 



40 



50 



"T 

| 65 



1 a 


2184 | 


1 b 


2569 | 


1 c 


9 1 


1 <3 


363 | 


1 e 


1171 | 


| f 


4 1 


1 a' 


1437 | 



693 


296 


86 


1543 


757 


243 


5.6 


2.8 





117 


48 


15 


696 


316 


107 


2.5 


1.1 





466 


200 


60 




Note: If the two strings are equal in 
length and end on a byte boundary, 
then a is replaced by a ' . 



Method : 

The current length of the target string 
is set equal to either the maximum of 
those of the operands or to the maximum 
length of the target field (when trunca- 
tion is necessary to avoid exceeding the 
length of this field) . The strings are 
'or'ed together for a length equal to the 
minimum of the lengths of the operands 
and the remainder of the longer string is 
moved into the target field up to the 
current length calculated for it; the 
remainder of the target field is left 
unchanged. 



Implementation : 

• Module size: 312 bytes 

• Execution time: 

Let Li. = the length of the shorter 
string 

L 2 = the difference between the 
string lengths 

M ± = FLOORULi- D/2048) 

N A = MOD (LjJ 2048) - 1 

Pi = FLOOR (Ni/8) 

S = SIGN(L 2 ) 

i = 1,2 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*M., + c*P, + S*(d + e+M-a + f*P 2 ) 



The 'Or* Operator (|) 

Module Name: IHEBSO 
Entry Point: IHEBSO0 
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Function: 

To implement the 'or* operation between 
two byte-aligned bit strings,, placing the 
result in a byte-aligned target field. 



Note: If the two strings are equal in 
length and end on a byte boundary, 
then a is replaced by a*. 
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The 'Not' Operator (■>) 



Bit String Concatenate/REPEAT 



Module Name: IHEBSN 
Entry Point: IHEBSNO 

Function: 

To implement the 'not' operator for a 
byte-aligned bit string, placing the 
result in a byte-aligned target field. 



Method: 

The current length of the target string 
is set equal to either the current length 
of the operand or to the maximum length 
of the target field (when truncation is 
necessary to avoid exceeding the length 
of this field) . The target field is set 
to a string of l's for a length equal to 
its calculated current* length and the 
result is obtained by an 'exclusive or' 
with the operand. The remainder of the 
target field beyond the calculated cur- 
rent length is left unchanged. 



Implementation : 

• Module size: 192 bytes 

• Execution time: 

Let Li = the string length 

Mi = FLOOR((L x - D/2048) 
N ± = MOD (1,1,2048) - 1 
Pi = FLOOR (Ni/8) 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*Mi + c*Pi 

or, if MOD(Li,8) = 0, from: 

a' + b*Mi + c*Pi 



Entry 
point 

IHEBSKK 
IHEBSKR 



Module Name: IHEBSK 
Entry Points: 

Operation 

Concatenate ( | | ) 
REPEAT (Bit string, n) 

Function: 



IHEBSKK: to concatenate two bit strings 
into a target field. 

IHEBSKR: to concatenate n + 1 instances 
of the single source string into a 
target field. If n £ 0, the result is 
the string itself. 



Method : 

The current length of the target field is 
made equal to the smaller of two values: 
the sum of the current lengths of the 
source strings, and the maximum length of 
the target field. Both entry points use 
a loop which obtains data from the source 
fields, aligns it correctly, and moves it 
to the target field. The remainder of 
the target field beyond the calculated 
current length is left unchanged. 



Implementation : 

• Module size: 328 bytes 

• Execution time: 

Let Li,L a = the lengths of the two 
strings 

Fi = FLOOR (Li/32) 

R =n+l, ifnSO 
=1, if n < 

i = 1,2 
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Then the approximate execution tiroes in 
microseconds for the System/360 models 
given below are obtained from the for- 
mula for the appropriate entry point: 

IHEBSKK: a + b* (F ± + F a > 

IHEBSKR: c + R*(d + b+F^) 

r t t ■ — t t -t n 

I I 30 | 40 | 50 | 65 | 75 | 
j. + x x + _ + ., 

| a | 3497 | 1311 | 445 | 126 | 74.7 | 
| b | 345 | 111 | 42 | 11.8 | 8.1 | 
| c | 594 | 187 | 73 | 16.4 | 11.0 | 
| d | 1505 | 479 | 196 | 56.3 | 33.2 | 
l x J. x x „x J 



Implementation : 

• Module size: 272 bytes 

• Execution time: 

Let L x = the number of bytes compared 
up to the first inequality 

L 2 = the number of bytes in the 
additional part of the longer 
string, compared to zeros if 
necessary 

M ± = FLOOR ((L x - D/256) 

N x = MOD (Ljl - 1,256) 

S = SIGN (L 2 ) 



Bit String Byte-aligned Comparison 

Module Name: IHEBSC 
Entry Point: IHEBSCO 

Function: 

To compare two byte-aligned bit strings 
and to return a condition code as bits 2 
and 3 of a full-word target field as 
follows: 



00 
01 

10 



if strings are equal 

if first string compares low at the 

first inequality 
if first string compares high at the 

first inequality 



The shorter string is treated as though 
extended with zeros to the length of the 
longer. 

The first byte of the target field is 
also used to preserve the program mask in 
the PSW for the calling routine. This 
byte contains: 



Bits 



Contents 



to 1 Instruction length code 01 

2 to 3 Condition code as above 

4 to 7 Program mask (calling routine) 



Method: 

The two strings are compared up to the 
current length of the shorter string. 
The remainder of the longer string is 
compared with zeros. 



Then the appropriate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*Mi + c+Ni + S*(d + e*L 2 ) 
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Bit String General Comparison 

Module Name: IHEBSD 
Entry Point: IHEBSD0 

Function: 

To compare two bit strings and return a 
condition code as bits 2 and 3 of a 
full-word target field as follows: 

00 if strings are equal 

01 if first string compares low at the 

first inequality 
10 if first string compares high at the 
first inequality 

The shorter string is treated as though 
extended with zeros to the length of the 
longer. 
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The first byte of the target field is 
also used to preserve the program mask in 
the PSW for the calling routine. This 
byte contains: 



Bits 

to 1 
2 to 3 
4 to 7 

Method : 



Contents 

Instruction length code 01 
Condition code as above 
Program mask (calling routine) 



The two strings are compared up to the 
current length of the shorter string. 
The remainder of the longer string is 
compared with zeros. 

Implementation : 

• Module size: 192 bytes 

• Execution time: 

Let F ± = the number of 32-bit words 
compared up to the first in- 
equality 

F 2 = the number of 32-bit words 
compared to zero if necessary 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b+Fi + c*F 2 

r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
j. + 1 + x x .| 

| a | 1113 | 344 | 120 | 31.7 | 18.2 | 
| b | 624 | 194 | 76 j 21.3 j 14.9 | 
| c | 437 | 139 | 55 | 15.3 | 11.4 | 

L J. X X X X J 



Bit String Assign/Fill 

Module Name: IHEBSM 
Entry Points: 

Operation 

Fixed-length assign 
Variable-length assign 
Zero fill only 



Entry 
point 

IHEBSMF 
IHEBSMV 
IHEBSMZ 



Function: 

IHEBSMF: to assign a byte-aligned string 
to a byte-aligned fixed-length target, 
filling out with zero bits if neces- 
sary. 



IHEBSMV: to assign a byte-aligned string 
to a byte-aligned variable- length tar- 
get. 



IHEBSMZ: to fill out the target area from 
its current length to its maximum 
length with zero bits. 

Method : 

IHEBSMF: the minimum of the source cur- 
rent length and the target maximum 
length is calculated and the source 
string is moved to the target for a 
length equal to this length. Zero 
filling of the target is performed if 
necessary. The current length of the 
target is set equal to the maximum 
length. 

IHEBSMV: the source string is moved to 
the target field as above, but without 
zero filling. The current length of 
the target is set appropriately. 

IHEBSMZ: zeros are propagated in the 
target from the current length to the 
maximum length. The current length of 
the target is set equal to the maximum 
length. 

Implementation : 

• Module size: 384 bytes 

• Execution time: 

Let Bj. = the bit length of the field to 
be assigned 

B 2 = the bit length for zero fill- 
ing 

Li = FLOOR (Ej/8) 

Mi = FLOOR (Li/2 56) 

Ni = MOD (Li, 256) 

S = SIGN (L 2 ) 

i = 1,2 



12 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the for- 
mula for the appropriate entry point: 
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f + b+Mx + c*N ± + S*(g + e*M 2 
+ C*N 2 ) 

a + b*MjL + c*N ± 

d + e*M 2 + c*N 2 
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Method: 

Arithmetic is performed according to the 
function definition, using the current 
length of the argument string. The 
result describes a fixed-length string. 



Implementation : 

• Module size: 192 bytes 

• Execution times: 

Approximate execution time in micro- 
seconds for the System/ 360 models given 
below is shown for each entry point: 

T T T T T T T 

|Entry | 30 | 40 | 50 | 65 | 75 | 

I Point | I I I I I 
j. x x x x x -I 

|IHEBSS2| 1184 | 378 | 152 | 42. 6| 28.8 | 
|. 1 x x x x -I 

JIHEBSS3J 1315 | 415 | 168 | 46. 7 j 31.7 j 
L X X X X X J 



Note: a' and f replace a and f respective- 
ly if MOD (B llf 8) = 



Other Information: 

This routine supplies assignment of byte- 
aligned bit strings of both fixed and 
variable lengths. Non-aligned strings 
may be assigned by using the REPEAT entry 
IHEESKR with n equal to 0. Any filling 
required for fixed length strings can 
then be obtained using the IHEBSMZ entry 
described above. 



BIT STRING FUNCTIONS 



Bit Str ing SUBSTR 



Module Name: 



Entry Points: 



IHEBSS 



Operation 

SUBSTR (Bit-string, i) 
SUBSTR (Bit-string,, i , j ) 

Function: 



Entry 
point 

IHEBSS2 
IHEBSS3 



To produce a string dope vector 
describing the SUBSTR pseudo-variable and 
function of a bit-string. 



Bit String INDEX 

Module Name: IHEBSI 
Entry Point: IHEBSI 

Function: 

To compare two bit strings to see if the 
second is identical to a substring of the 
first, and, if it is, to produce a binary 
integer (the index) which indicates the 
first bit position in the first string at 
which such a substring begins. If no 
such index is found, or if either string 
is null, the function value is zero. 

Method: 

The index is found by shifting and com- 
paring portions of the two strings in 
registers. 
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Implementation : 

• Module size: 296 bytes 

• Execution time: 

Let F^ = the length (in words) actually 
processed at the ith compari- 
son (the length up to the 
first inequality or the length 
of the second string if no 
inequality is found) 

B = (length of first string) 
(length of second string) 

I = the resulting index, except 
when this is with B > 0, in 
which case I = B + 1 

J = MOD(I,32) 

i = l f n 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 



a + b*I + c*J + 



I 
d*E 
i=l 



(F< - 1) 



r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
J. + x + x + ^ 

| a | 2026 j 646 | 246 | 63.6 | 37.4 | 
| b | 422 | 124 | 48 | 13.9 j 9.2 | 
| c | 104 | 41 | 18 | 4.3 | 3.1 | 
| d | 543 | 180 | 73 j 18.9 | 13.2 | 
L J. — , x x . x x j 



Method: 

The current length of the target string 
is set equal to either the maximum of the 
current lengths of the source strings or 
to the maximum length of the target field 
(when truncation is necessary to avoid 
exceeding the length of this field) . The 
necessary operation is performed on the 
strings and the result stored in the 
target field. If one string is shorter 
than the other, it is regarded as being 
extended on the right with zeros up to 
the length of the longer. The field 
between the calculated current length and 
the maximum length of the target is left 
unchanged. 

Implementation : 

• Module size: 480 bytes 

• Execution time: 

Let B ± = the bit length of the shorter 
string 

B 2 = the difference in bit lengths 
of the strings 

F L = CEIL(B.j/32) 

i = 1,2 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*Fa. + c*F a 



Bit String BOOL (Boolean Function) 



Module Name: 
Entry Point: 
Function: 



IHEBSF 



IHEBSF0 



To take two source strings and perform 
one of the sixteen possible logical oper- 
ations between corresponding bits. The 
particular operation performed is defined 
by inserting the bit pattern - T\ x n a n 3 n H '- 
yielded by the third argument into the 
table below: 

r t t t t 1 

| First field | | | 1 | 1 | 
|. + + + x ^ 

| Second field | | 1 | | 1 | 
j. + x x x 1 

| Target field | nj. | n a | n 3 | n^\ 

t x x x x — _j 



r t t t t t 1 

| | 30 | 40 | 50 | 65 | 75 | 

L X X X X X J 

| a | 2275 | 767 | 291 | 75.0 | 46.5 | 
| b | 728 | 230 | 95 | 30.0 | 19.4 j 
| c | 501 | 162 | 68 | 21.8 | 15.2 | 

L X X X X X J 



CHARACTER STRING OPERATIONS 



Character String Concatenate/REPEAT 



Module Name: 


IHECSK 




Entry Points: 






Operation 




Entry 
point 


Concatenate 
REPEAT 


<ll> 


IHECSKK 
IHECSKR 



(Character string, n) 
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Functions 



Character String Compare 



IHECSKK: to concatenate two character 
strings into a target field. 



IHECSKR: to concatenate n + 1 instances 

of the single source string into a 

target field. If n < 0, the result is 
the string itself. 



Method: 

The current length of the target field is 
made equal to the smaller of two values: 
the sum of the current lengths of the 
source fields, and the maximum length of 
the target field. The source strings are 
then moved to the target. Characters 
beyond the range of the target current 
length remain unaltered. 



Implementation : 

• Module size: 208 bytes 

• Execution time: 

Let Li. f L 2 = the lengths of the source 
strings 

M A = FLOOR((Li~ D/256) 

N^ = MODCLi - 1,256) 

R = n + 1, if n £ 
=1, if n < 

i = 1,2 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

IHECSKK: a + b*(M ± + M a > + c*(N ± + N a ) 

IHECSKR: R*(b*Mj. + c+Nj. + d) + e 



Module Name: 



IHECSC 



1 




"T- 
1 


30 


■T" 
1 


40 




"T- 

1 


50 




"T" 
1 


65 


-T- 
1 


75 


"1 

1 


h 


a 

b 


■+■ 


1407 
1196 


■+- 


491 
706 




-+- 


209 
319 




-+- 


63.7 
108 


-+- 


46.4 
95.8 


H 




c 




4 




2. 


5 




1. 


1 




0.4 




0*3 






d 




703 




236 






101 






31.7 




22.0 






e 




108 




54. 


6 




26. 







5.4 




6.2 




L. 




.i. 




.1. 






..1. 






.X. 




.1.. 




.j 



Entry Point: IHECSC0 

Function: 

To compare two character strings and to 
return a condition code as bits 2 and 3 
of a full-word target field as follows: 

00 if strings are equal 

01 if first string compares low at the 

first inequality 
10 if the first string compares high at 
the first inequality 

The shorter string is treated as though 
extended with blanks to the length of the 
longer one. 

The first byte of the target field is 
also used to preserve the program mask in 
the PSW for the calling routine. This 
byte contains: 



Bits 



Contents 



to 1 Instruction length code 01 

2 to 3 Condition code as above 

4 to 7 Program mask (calling routine) 



Method; 

The two strings are compared in storage. 
If the strings are of different lengths 
and are identical up to the length of the 
shorter, the remainder of the longer is 
compared with blanks. 



Implementation : 

• Module size: 200 bytes 

• Execution time: 

Let Li = the length of the strings com- 
pared up to the first in- 
equality (proceeding left to 
right) 

L a = the length of the additional 
part of the longer string com- 
pared with blanks if necessary 

Mi = FLOOR ((Li- D/256) 

Ni = MOD(Li - 1,256) 

S = SIGN(L 2 ) 

i = 1,2 
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Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b+Mi. + c*Ni + S*(d + e*M 2 + c*N 2 ) 



I I 



30 



40 



50 



65 



I 75 | 



a 




849 




284 






116 






36 


8 


| 27. 


1 




b 




1469 




790 






289 






114 




|110 






c 




5 




2. 


8 




1. 










4 


1 o. 


4 




d 




620 




210 






88 






27 


3 


| 21. 


7 




e 


X. 


1474 


.x. 


794 




_X. 


290 




.X. 


114 




|110 

.X 




..l 



IHECSMV: moves the string as above, but 
without blank filling. The current 
length of the target is set appropri- 
ately. 



IHECSMB: propagates blanks and sets the 
current length of the target equal to 
its maximum length. 



IHECSMH, IHECSML: uses part of the blank 
fill routine to propagate the highest 
or lowest character in the collating 
sequence up to the current length of 
the target. 



Character String Assign/Fill/HIGH/LOW 

Module Name: IHECSM 
Entry Points: 





Entry 


Operation 


point 


Fixed-length assign 


IHECSMF 


Variable-length assign 


IHECSMV 


Blank fill only 


IHECSMB 


HIGH 


IHECSMH 


LOW 


IHECSML 



Function: 

IHECSMF: to assign a character string to 
a fixed-length target, filling out with 
blanks if necessary. 

IHECSMV: to assign a character string to 
a variable-length target. 

IHECSMB: to fill out the target field 
from its current length to its maximum 
length with blanks. 

IHECSMH: to fill a target field with the 
highest character in the collating 
sequence, up to its current length. 

IHECSML: to fill the target field with 
the lowest character in the collating 
sequence, up to its current length. 

Method: 

IHECSMF: The minimum of the source cur- 
rent length and the target maximum 
length is calculated and the source 
string is moved to the target for a 
length equal to this length. Filling 
of the target with blanks up to the 
target maximum length is performed if 
necessary. The current length of the 
target is set equal to its maximum 
length. 



Implementation : 

• Module size: 280 bytes 

• Execution time: 

Let L ± = either the specified length 
(for IHECSMH/L/B) or the 
length for blank filling 
(IHECSMF) 

L 2 = the length of the shorter of 
the source and target fields 
( IHECSMF/V) 

M-l = FLOOR((Li - 2)/256) 

Ni = MOD(Li - 2, 256) (where L ± > 2 
for both Mi and N A ) 

M 2 = FLOOR((L 2 - D/256) 

N 2 = MOD(L 2 - 1, 256) 

S = SIGN(Li) 
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Then the approximate execution times in 
microseconds for the System /360 models 
given below are obtained from the fol- 
lowing formulas: 

IHECSMF: h + g*M 2 + c*N 2 + S*(i + b*M x 
+ c+Nj.) 

IHECSMV: f + g*M 2 + c*N 2 

IHECSMB: e + b+M^. + c+N^. 

IHECSMH: a + b*M ± + c+N^ 

IHECSML: d + b*M x + c*N ± 



30 



40 



50 



65 



75 



1 a 


830 | 


1 fc 


1171 | 


1 c 


4 | 


1 d 


799 | 


1 e 


892 | 


| f 


790 | 


1 g 


1196 | 


1 h 


1298 | 


| i 


400 | 



282 


120 




37.0 


696 


316 




107 


2.5 


1. 


1 


0.4 


268 


116 




35.5 


300 


130 




40.5 


267 


114 




34.6 


706 


319 




108 


446 


190 




57.4 


129 


57. 


6 


18.5 



| 26. 


8 1 


| 95. 


2 1 


1 o. 


3 1 


| 25. 


6 1 


| 28. 


o 1 


| 24. 


9 | 


| 95. 


8 1 


| 40. 


6 1 


| 13. 


7 | 



Implementation : 

• Module size: 176 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/ 360 models given 
below are as shown for each entry 
point: 

r T T T T T 1 

|Entry | 30 | 40 | 50 | 65 | 75 | 

I Point | I I I I I 
j. + 1 + + + -i 

JIHECSS2J 887 j 310 j 127 j 36.9 j 26.2 j 
J. + + 1 + + .| 

|IHECSS3| 1018 | 347 | 143 | 41. 0| 29.1 | 

L ± X X X J. J 



Character String INDEX 

Module Name: IHECSI 
Entry Point: IHECSI0 



CHARACTER STRING FUNCTIONS 

Character String , SUBSTR 

Module Name: IHECSS 
Entry Points: 



Operation 

SUBSTR (Character-string,, i) 
SUBSTR (Character-string , i„ j ) 



Entry 
point 

IHECSS2 
IHECSS3 



Function: 

To produce a string dope vector describ- 
ing the SUBSTR pseudo-variable and func- 
tion of a character string. 



Function: 

To compare two chara 
if the second is ident 
of the first, and, if 
binary integer (the 
ates the first charact 
first string at whi 
begins. If no such in 
either string is null, 
is zero. 



Method: 



cter strings to see 
ical to a substring 
it is, to produce a 
index) which indic- 
er position in the 
ch such a substring 
dex is found, or if 
the function value 



The point required is located by compar- 
ing in storage the second string with a 
corresponding number of characters in the 
first string. 



Method: 

Arithmetic is performed according to the 
function definition, using the current 
length of the argument string. The 
result describes a fixed-length string. 
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Implementation : 

• Module size: 168 bytes 

• Execution time: 

Let L i = the length processed at the 
ith comparison (the length up 
to the first inequality, or 
the length of the second 
string if no inequality is 
found) 

M i = FLOOR (Li/256) 

Nj_ = MOD(Li, 256) 

B = (length of first string) 
(length of second string) 

I = the resulting index, except 
when this is with B £ 0, in 
which case I = B + 1 

i = l,#n 



Then the approximate execution times in 
microseconds for the System/360 models 
given below is obtained from the fol- 
lowing formula: 



a + b*I + 2 (c*Mi + d*Ni) 
i=l 



r T T T T T 1 

| | 30 | 40 | 50 | 65 j 75 | 

| a | 954 | 314 | 127 | 38.2 | 24.5 | 
i b | 160 | 59.9 | 28.8 | 10.1 | 9.3 | 
j c | 1427 | 776 | 280 | 110 |107 | 
| d | 5 | 2.8 | 1.0 | 0.4 | 0.4 | 
i J. J. x i i J 
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CHAPTER 2: ARITHMETIC OPERATIONS AND FUNCTIONS 



Litrary arithmetic modules support all 
those arithmetic generic functions and 
operators for which the compilers neither 
produce in-line code nor (as for the func- 
tions FIXED, FLOAT, BINARY and DECIMAL) use 
parts of the conversion package. I — 

Statistics for accuracy of floating- 
point modules are given where considered 
meaningful and helpful; an explanation of 
their use is given in the chapter on 
mathematical routines. Precise results are 
obtained from all fixed-point modules 
except complex division and complex ABS, 
where small truncation errors inevitably 
occur, and the ADD function (fixed 
decimal), in which the effect of truncation 
errors depends on the relative values of 
the scale factors of the arguments. — 

Any restrictions on the admissibility of 
arguments are noted under the headings 
•Range' and 'Error and Exceptional 
Conditions'. H — 

Range; This states any ranges of arguments 
which a module assumes to have been exclud- 
ed prior to its being called. 

Error and Exce p tional Conditions : These Speed 
cover conditions which may result from the 
use of a routine; they are listed in four 
categories: 



P — Programmed conditions in the module 
concerned. Programmed tests are 
made where this is not too costly 
and,, if an invalid argument is 
found, a branch is taken to the 
entry point IHEERRC of the execution 



error package (EXEP) . This results 
in the printing of an appropriate 
message and in the ERROR condition 
being raised. 

Interrupt conditions in the module 
concerned. For those routines where 
SIZE and FIXEDOVERFLOW are detected 
by programmed tests or where hard- 
ware interruptions may occur, the 
OVERFLOW, UNDERFLOW, FIXEDOVERFLOW, 
SIZE and ZERODIVIDE conditions pass 
to the ON handler (IHEERR) and are 
treated in the normal way. The 
machine is assumed to be enabled for 
all interruptions except signifi- 
cance, which is masked off. 

Programmed conditions in modules 
called by the module concerned. 
These occur when invalid arguments 
are detected in the module called. 

As I, but the interrupt conditions 
occur in the modules called by the 
module concerned. 



The average execution times given are 
based on the IBM System/ 360 Instruction 
Timing Information , Form A22-6825. These 
times include the times taken by the 
modules called. 

A summary of the Library arithmetic 
modules is given in Figures 2 and 3. 
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ARITHMETIC OPERATIONS 



Operation | Binary 

| fixed 

_ _ j. 


1 

1 
i 


Decimal | 
fixed | 


Short 
float 


1 

1 
i 


Long | 
float j 


Real Operations j 


T - 

Integer exponentiation: x**n | IHEXIB 
General exponentiation: x**y | 
Shift-and-assign, Shift-and-load j 

_x 


i 

1 

1 

1 
i 


IHEXID | 

1 
IHEAPD j 
. x_ 


IHEXIS 
IHEXXS 


i 

1 

1 

1 
l 


IHEXIL | 
IHEXXL j 


Complex Operations 


"T" 


T- 









Multiplication/division: Zi*Z2,Zi/z 2 | IHEMZU | IHEMZV | 

Multiplication: z ± *z 2 | - j j 

Division: z ± /z 2 j j - j 

Integer exponentiation: z**n | IHEXIU | IHEXIV j 

General exponentiation: z ± **z 2 j j j 

JL X X- 



IHEMZW | IHEMZZ | 

IHEDZW | IHEDZZ | 

IHEXIW j IHEXIZ j 

IHEXXW j IHEXXZ | 



Figure 2. Arithmetic Operations 



ARITHMETIC FUNCTIONS 



Function 


1 
1 

1 


T T 

Binary | Decimal | Short 
fixed | fixed j float 


1 

1 
i 


Long j 
float j 


Real Arguments j 


MAX, MIN 
ADD 


1 

1 
1 

1 


T T 

IHEMXB j IHEMXD j IHEMXS 

| IHEADD | 
X X 


i 

1 

1 
i 


IHEMXL | 






Complex Arguments 







ADD | - | IHEADV | - | - | 
MULTIPLY j IHEMPU j IHEMPV j j - | - | 
DIVIDE | IHEDVU j IHEDVV j - | - | 
ABS | IHEABU j IHEABV | IHEABW j IHEABZ j 
X X X X J 



Figure 3. Arithmetic Functions 



REAL OPERATIONS 



Method ; 



Positive Integer Exponentiation (fixed 
binary) 



Module Name: IHEXIB 



The result is set initially to the value 
of the argument. The final result is 
then obtained by repeated squaring of 
this value or squaring and multiplying by 
the argument. 



Entry Point: IHEXIBO 



Range: 



Function: 

To calculate x**n, where n is a positive 
integer. 



< n < 2**31 

The precision rules of PL/I impose a 
further restriction in that if x has a 
precision (p,q), this module will be 
called only if n*(p + 1) - 1 < 31. This 
implies that n < 32/ (p + 1) < 16 for all 
such cases. 
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Implementation : 

• Module size: 88 bytes 

• Execution time: 

Let M = number of significant bits in 
the exponent 

N = number of 1 bits in the expo- 
nent 

Then the approximate execution times in 
microseconds for the System /360 models 
given below are obtained from the fol- 
lowing formula: 

-a + b*M + c*N 



Implementation : 

• Module size: 136 bytes 

• Execution time: 

Let M = number of significant bits in 
the exponent 

N = number of 1 bits in the expo- 
nent 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

-a + b*M + c*N 



r T T — T T -T T 

| | 30 | 40 | 50 | 65 | 75 | 

|. + + — + x„— — .X 4 

| a | 238 | 9 | 5.5 | -1.46| -0.6 | 
| b | 708 | 188 | 63.8 | 15.6 | 9.8 j 
| c | 335 | 94 | 33.0 | 6.1 | 3.9 | 
i x x -_ J. x_ .x J 



r T T T T T T 

| | 30 | 40 | 50 | 65 | 75 | 
|. + x x X + .] 

| a j 580 | 48 | -16.3 | 37.1 | 44.5 j 
| b | 1113 | 279 | 91 | 45.5 | 40.3 | 
| c | 752 | 192 | 65 | 38.8 | 37.2 | 

L X X X X 1 X J 



Positive Integer Exponentiation (fixed 
decimal) 



Module Name: 



IHEXID 



Entry Point: IHEXID0 

Function: 

To calculate x**n,, where n is a positive 
integer. 

Method: 

The result is set initially to the value 
of the argument. The final result is 
then obtained by repeated squaring of 
this value or squaring and multiplying by 
the argument. 

Range: 

The precision rules of PL/I impose the 
restriction that if x has a precision 
<P»q)» this module will be called only if 
n*(p + 1) - 1 < 15. This implies that 
n < 16/ (p + 1) < 8 for all such cases 
and, in fact, this module will operate 
only for the range < n < 8. 



Integer Exponentiation (floating-point) 



Module Names and Entry Points : 



Argument 

Short float 
Long float 

Function: 



Module 
name 

IHEXIS 
IHEXIL 



Entry 
point 

IHEXIS0 
IHEXIL0 



To calculate x**n, where n is an integer 
between -2**31 and 2**31 - 1 inclusive. 

Method: 

If the exponent is zero and the argument 
non-zero, the result 1 is returned 
immediately. Otherwise the result is set 
initially to the value of the argument 
and the exponent is made positive. The 
argument is raised to this positive power 
by repeated squaring of the contents of 
the result field or squaring and multi- 
plying by the argument. Then, if the 
exponent was negative, the reciprocal of 
the result is taken, otherwise it is left 
unchanged. 
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Accuracy: 

The values given here are for the rela- 
tive error divided by the exponent for 
exponents between 2 and 1023; the argu- 
ments are uniformly distributed over the 
full range for each exponent for which 
neither OVERFLOW nor UNDERFLOW occurs. 
There are 2** (10 - k) arguments for each 
exponent in the range 2**k £ exponent £ 
2**(k + 1) - 1, where k has integral 
values from 1 to 9 inclusive. 

IHEXIS 

r t 1 

| R.M.S. relative | Maximum relative | 
J error/exponent | error/exponent j 
J *io**6 j *10**6 j 

j. + .| 

| 0.00871 | 0.692 | 

L X J 

IHEXIL 

r t 1 

| R.M.S. relative | Maximum relative | 
| error/exponent j error/exponent j 
| *10**15 | *10**15 | 

j! + 1 

| 0.0995 | 1.73 | 

L X J 



Error and Exceptional Conditions : 

P : x = with n < 

I : OVERFLOW, UNDERFLOW 

Since x**(-m), where m is a positive 
integer, is evaluated as l/(x**m), 
the OVERFLOW condition may occur when 
m is large, and the UNDERFLOW condi- 
tion when x is very small. 

Implementation : 

• Module size: IHEXIS 152 bytes 

IHEXIL 152 bytes 

• Execution time: 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

a + b*M + c*N for positive exponents 

a'+ b*M + c*N for negative exponents 

IHEXIS 

r T T * T T T 1 

| I 30 I 40 | 50 | 65 | 75 | 
h + x x + x -I 

| a | -104 | 23 | 29 | 10.7 | 7.6 | 
| b | 701 | 176 | 56 | 14.3 | 8.7 | 
| c | 342 | 90 | 26 | 5.7 | 3.2 | 
| a« | 552 | 171 | 57 | 18.6 | 12.7 | 

L X X X X X J 

IHEXIL 
T T T T T T 1 

| | 30 | 40 I 50 | 65 | 75 | 
j. x x x + x -i 

| a |-1535 |-322 | 0.7 | 4.1 | 3.7 | 
| b | 1441 j 355 | 73 | 17.5 | 10.7 | 
| c | 1082 | 269 | 42 | 8.9 | 5.2 | 
| a' | 1055 j 180 | 78 | 20.0 | 11.9 | 
L X X X X X J 



Other Information: 

IHEXIS: For large exponents, tor example, 
those greater than 1023, it is general- 
ly faster and more accurate to use the 
module IHEXXS rather than IHEXIS, pass- 
ing the exponent as a floating-point 
argument. However,, it should be noted 
that IHEXXS will not accept a negative 
first argument, and thus it is neces- 
sary to pass the absolute value of this 
argument, and also, in cases where the 
exponent is odd, to test the sign of 
the argument in order to be able to 
attach the correct sign to the numer- 
ical result returned. 



Let M = number of significant bits in 
the exponent 

N = number of 1 bits in the expo- 
nent 



General Floating-Point Exponentiation 



Module Names and Entry Points: 



Arqument 


Module 
name 


Entry 
point 


Short float 
Long float 


IHEXXS 
IHEXXL 


IHEXXS0 
IHEXXL0 


Function: 







To calculate x**y, where x and y 
floating-point numbers. 



are 
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Method : 



Method: 



When x = 0„ the result x**y = is given 
if y > 0.i and an error message if y :S 0. 
When x * and y = 0, the result x**y = 1 
is given. Otherwise x**y is computed as 
EXP(y*LOG(x)), r using the appropriate 
mathematical function routines. 



Error and Exceptional Conditions: 

P : x = with y £ 

: a. x < with y * : error caused in 
LOG routine 

b. y*LOG(x) > 174.673: error caused 
in EXP routine 



Implementation : 

• Module size: IHEXXS: 144 bytes 

IHEXXL: 152 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

r t t t t t— ■ n 

| Module j 30 J 40 j 50 j 65 j 75 j 

I Name | | | III 
(. + __x + x x .) 

| IHEXXS | 9809 | 2861 j 902 | 236 | 143 j 
| IHEXXL |30444 | 7453 | 1579 | 358 ] 203 | 

L X J. X .X X J 



The argument scale factor is subtracted 
from the target scale factor. The argu- 
ment is converted to precision 31 in a 
field with a shift equal to the magnitude 
of the difference between the scale fac- 
tors; the shift is to the left if the 
difference is positive and to the right 
if negative. 

If entry point IHEAPDB is used, the field 
is moved unchanged to the target. If 
entry point IHEAPDA is used, the result 
is checked for FIXEDOVERFLOW and then 
assigned to the target with the specified 
precision. The assignment may cause the 
SIZE condition to be raised. 



Error and Exceptional Conditions: 
I : FIXEDOVERFLOW or SIZE 

Implementation : 

• Module size: 360 bytes 

• Execution times: 

Let S = (Target scale factor) 
(argument scale factor) 

fa. = if S £ 31 
= 1 if S < 31 

fa = if ABS(S) > p 
= 1 if ABS(S) < p 

(where p is the precision of 
the argument) 



Shift-and-assiqn, Shift-and-load (fixed 
decimal) 



Module Name: IHEAPD 
Entry Points: 

Operation 

Shift and assign 
Shift and load 



Entry 
point 

IHEAPDA 
IHEAPDB 



Function: 

IHEAPDA: To convert a real fixed decimal 
number with precision (pi,qi) to preci- 
sion (Pa/qa), where pj. S 31 and 
Pa * 15. 

IHEAPDB: To convert a real fixed decimal 
number with precision (pi, q A ) to preci- 
sion (31, q 2 ), where p*. < 31. 
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Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 



IHEAPDA : 

(i) S > 

(ii) S = 

(iii) S < 
IHEAPDB : 

(i) S > 

(ii) S = 

(iii) S < 



t« + t a *f ± 

t 7 

i a + t *f a 

ti + t 2 *fi 

t 3 

t* + t B + f a 



r 


" T - 

1 


30 


-T- 


40 


"T" 


50 




-T - 

1 


65 


■T" 


1 

75 | 


1 ti 


1 


1314 


"+- 
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■+- 


241 




-+- 
1 


68.4 


-+- 


53.5 | 


1 t 3 


-+- 

1 
1 


451 


I 
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-+- 
1 


64. 


3 


-+- 

1 

1 


17.1 


•+- 

1 


15.1 | 


1 t 3 


1 

1 
— X- 


1117 


1 
| 


446 


1 
1 
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i 

1 
| 
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1 
1 


47.2 | 
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t 
1 

1 


1220 
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386 


1 
-+- 
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1 
1 

1 


49.5 


1 

■+- 


36.9 | 
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502 


221 


94. 


2 


37.3 


21.8 | 




I 
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| 
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1 




1 te 


1 
1 


1996 


1 


736 


1 


345 




1 
1 


95.6 


1 


70.9 | 


1 t 7 


1 


1799 


-+- 


670 


"+- 


305 




_X. 
1 


87.3 


-+- 


64.5 | 


1 te 


1 


1902 


-+- 


610 


-+- 


274 




1 


76.7 


-+- 


54.4 | 








.X. 
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COMPLEX OPERATIONS 

Multiplication/Division (fixed binary) 

Module Name: IHEMZU 
Entry Points: 



Mathematical 
Operation 


Entry 
point 


z ± *z a 

Zi/Z 2 


IHEMZUM 
IHEMZUD 


net ion: 





To calculate Zj.*z a or Zi/z aif where z x and 
z a are fixed-point binary complex num- 
bers. 



Method: 

Let z ± = a + bl and z a = c + dl. Then, 
for multiplication, an incorporated sub- 
routine is used to compute a*c - b*d and 
b*c + a*d; these are tested for FIXED- 
OVERFLOW and then stored as the real and 
imaginary parts of the result. 

For division, the subroutine is used to 
compute a*c + b*d and b*c - a*d. The 
expression c**2 + d**2 is computed and 
the real and imaginary parts of the 
result are then obtained by division. 

The subroutine computes the expressions 
u*x + v*y and v*x - u*y. 



Error and Exceptional Conditions: 

I : FIXEDOVERFLOW in either routine, 
ZERODIVIDE in the division routine. 



Implementation : 

• Module size: 240 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the following 
formulas and tables: 

r ■ T T T T T T 

I Entry | 30 | 40 | 50 | 65 | 75 | 
I Point | I I I I I 

j. + 1 x x x -J 

|IHEMZUM| 2421 | 689 | 256 | 56. 6| 37.1 | 

L X X X X X J 

IHEMZUD 

Let M = number of significant bits in 
Z a *CONJG(Z a ) 

N = FLOOR (M/4 - 8) 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b + c + N*d 

r T T T T T T 

| | 30 | 40 | 50 | 65 | 75 | 
^ x x x + x .| 

j a J 3021 J1340 j 420 j 94.1 j 64.2 j 

I b I 338 j 79 j 30 j 5*7 | 2.7 | 

| c | 78 | 24 | 18 | 6.6 | 2.1 | 

| d | 213 | 56 | 22 | 6.2 | 3.8 | 

L X X X X X J 

Note: b = if M < 31 
c = d = if M ^ 32 
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Multiplication/Division (fixed decimal) 



Module Name: 


IHEMZV 


Entry Points: 




Mathematical 
Operation 


Entry 
point 


Zi*Z 2 

z ± /z 2 


IHEMZVM 
IHEMZVD 


Function: 





To calculate z i *z 2 or z x /z 2 where z ± and 
z 2 are fixed-point decimal complex num- 
bers. 



Method: 

Let z x = a + bl and z a = c + dl. The 
products a*c f b*c, a*d and b*d are com- 
puted. Then the required result is 
obtained as follows: 

Multiplication : 

Real part a*c - b*d 
Imaginary part b*c + a*d 

Division: 

Real part (a*c + b*d)/(c*c + d*d) 
Imaginary part (b*c - a*d)/(c*c + d*d) 



Then approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

IHEMZVM: a + b+Lj. + c*L 2 + d*Li*L 2 



IHEMZVD: 



30 



e + f+Lj. + g*L 2 + h*Li*L 2 + 
j*L 2 **2 + T*(k + m*L 2 ) 



40 



50 



65 



I 75 | 



a 


2605 


956 


457 


128 


103 | 


b 


246 


105 


49 


17 


13 | 


c 


30 


33 


19 


5.3 


3.5 | 


d 


112 


15 





4.0 


4.7 | 


e 


14525 


4203 


1144 


626 


508 | 


f 


246 


105 


49 


17 


12 | 


q 


168 


103 


41 


16 


11 | 


h 


112 


15 





4.0 


4.7 | 


i 


56 


7.5 





2.0 


2.3 | 


k 


731 


251 


131 


41.2 


30.2 | 


m 


12 


16 


4.3 


0.8 


0.3 1 





Other Information: 

It should be noted from the timings for 
multiplication that where the operands 
differ in precision, it is faster to 
present the longer operand as the second 
argument rather than the first. 



Error and Exceptional Conditions: 

I : FIXEDOVERFLOW in either routine, 
ZERODIVIDE in the division routine. 



Multiplication (floating-point) 



Implementation : 

• Module size: 672 bytes 

• Execution time: 

Let (p,q)„ (r,s) = precisions of the 
operands 

L x = FLOOR (p/2) + 1 

L 2 = FLOOR (r/2) + 1 

T = c**2 + d**2 2: 10**(15 - 2*s) 

(i.e., T = 1 or depending on 
whether the relation is true 
or false) 



Module Names and Entry Points: 



Argument 



Module 
name 



Short float IHEMZW 
Long float IHEMZZ 



Function: 



Entry 
point 

IHEMZW0 
IHEMZZ0 



To compute z 1 *z 2 in floating-point, when 
Zi = a + bl and z 2 = c + dl. 



Method : 

The real and imaginary parts of the 
result are computed as a*c - b*d and 
b*c + a*d, respectively. 



Error and Exceptional Conditions: 
I : Exponent OVERFLOW and UNDERFLOW 
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Implementation : 

• Module size: IHEMZW 64 bytes 

IHEMZZ 64 bytes 

• Execution times: 

Approximate execution time in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

(Module | 30 | 40 | 50 | 65 | 75 | 
I Name I | | I I I 
j. x x x + x 1 

| IHEMZW | 1979 | 550 | 172 | 41.9 | 23.3 | 

|. + + + + + ^ 

| IHEMZZ | 5115 | 1307 | 251 j 62.3 | 31.3 | 

L X X X X X J 



Implementation : 

• Module size: IHEDZW 104 bytes 

IHEDZZ 104 bytes 

• Execution times : 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

r T T T T T 1 

j Module j 30 j 40 j 50 j 65 j 75 j 
I Name I I I I I | 
|. x x x — x x .] 

|IHEDZW | 3546| 875 | 221 | 60. 8| 35.7 | 
|. + x x x + .| 

|IHEDZZ | 11741| 2515 | 234 | 92. 5| 51.1 | 

L X X X X X J 



Division (floating-point) 



Module Names and Entry Points: 



Argument 



Module 
name 



Entry 
point 

IHEDZW0 
IHEDZZ0 



Short float IHEDZW 
Long float IHEDZZ 

Function: 



To compute Zi/z 2 in floating-point, when 
z ± = a + bl and z 2 = c + dl. 



Method: 

1. ABS(c) Z ABS(d) 

Compute g = d/c 

then REAL (z ± /z 2 ) = (a + b*q)/(c + d*q) 
IMAG (Z;,./z 2 ) = (b - a*q)/(c + d*q) 

2. ABS(c) < ABS(d) 

Ca + bl)/(c + dl) = (b - al)/(d - cl) , 
which reduces to the first case. 

The comparison between ABS(c) and ABS(d) 
is adequately performed in short preci- 
sion in both modules. 



Error and Exceptional Conditions : 

I : OVERFLOW, UNDERFLOW and ZERODIVIDE 



Positive Integer Exponentiation (fixed 
binary) 



Module Name: IHEXIU 

Entry Point: IHEXIU0 

Function: 

To" calculate z**n, where n is a positive 
integer less than 2**31. 

Method: 

The contents of the target field are set 
to the value of z. The final result is 
obtained by repeated squaring of the 
contents of the target field or squaring 
and multiplying by z. Multiplication is 
performed by the complex multiplication 
routine IHEMZU. 

Range: 

< n < 2**31. 

The precision rules of PL/I impose a 
further restriction in that if z has a 
precision (p, q) , this module may only be 
called if n*(p + 1) - 1 < 31. This 
implies that n < 32/ (p +1) < 16 for all 
such cases. 

Implementation : 

• Module size: 120 bytes 

• Execution time: 

Let M = number of significant bits in 
the exponent 

N = number of 1 bits in the expo- 
nent 
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Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

-a + b*M + c*N 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

-a + b*M + c*N 



r T T T" T T 1 

| | 30 | 40 | 50 | 65 | 75 | 

| a | 4169 | 1183 | 405 | 101 | 69.5 | 
| b | 2409 j 822 | 306 i 73.3 | 48.3 j 
| c | 2553 | 738 | 276 | 62.7 j 42.1 | 
L X X X X X J 



r T T T T T 1 

| | 30 | 40 I 50 | 65 | 75 | 
|~— +— + +— + + H 

j a j 9200 | 3100 j 1300 j 450 | 370 | 
| b | 5000 | 1700 | 730 | 250 | 200 | 
| c | 5500 | 1800 | 750 | 250 | 200 | 
L X X X X X J 



Positive Integer Exponentiation (fixed 
decimal) 



Module Name: IHEXIV 
Entry Point: IHEXIV0 

Function: 

To calculate z**n, where n is a positive 
integer less than 2**31. 

Method: 

The contents of the target field are set 
to the value of the argument. The final 
result is obtained by repeated squaring 
of the contents of the target field or 
squaring and multiplying by the argument. 
Multiplication is performed by the com- 
plex multiplication routine IHEMZV. 

Range: 

The precision rules of PL/I impose the 
restriction that if z has a precision 
(p#<3)# this module may only be called if 
n*(p + 1) - 1 < 15. This implies that 
n £ 16/ (p + 1) < 8 for all such cases 
and, in fact, this module will operate 
only for the range < n < 8. 

Implementation : 

• Module size: 192 bytes 

• Execution time: 

Let M = number of significant bits in 
the exponent 

N = number of 1 bits in the expo- 
nent 



Integer Exponentiation (floating-point) 



Module Names and Entry Points: 



Argument 

Short float 
Long float 

Function: 



Module 
name 

IHEXIW 
IHEXIZ 



Entry 
point 

IHEXIW0 
IHEXIZ0 



To calculate z**n, where n is an integer 
between -2**31 and 2**31 - 1 inclusive. 

Method: 

If the exponent is and the argument 
non-zero, the answer 1 is returned 
immediately. If the exponent is non- 
zero, the contents of the target field 
are set to the argument value. The 
exponent is made positive and the 
argument raised to this positive power by 
repeated squaring of the contents of the 
target field or squaring and multiplying 
by the argument. Multiplication is per- 
formed by a branch to the complex multi- 
plication routine. Then, if the exponent 
was negative, the reciprocal of the 
result is taken, otherwise it is left 
unchanged. 

Error and Exceptional Conditions: 

P : z = with n < 

I : OVERFLOW, UNDERFLOW 

Since x** (-m) , where m is a positive 
integer, is evaluated as l/(x**m), 
the OVERFLOW condition may occur when 
m is large and the UNDERFLOW condi- 
tion when x is very small. 

H : OVERFLOW or UNDERFLOW in complex 
multiplication routine (IHEMZW or 
IHEMZZ) 
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Implementation : 



Method : 



• Module size: IHEXIW 256 bytes 

IHEXIZ 256 bytes 

• Execution time: 

Let M = number of significant bits in 
the exponent 

N = number of 1 bits in the expo- 
nent 



When Zjl = , the result is returned if 
REAL(z 2 ) > and IMAG(z 2 ) = 0. 
Otherwise, z ± **z 2 is computed as 

EXP(z 2 *LOG(z 1 )) f 

with the proviso that if IMAG(zj.) = 
then LOGCABS (za.) > is calculated by a call 
to the real LOG routine, not to the 
complex LOG routine. 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

-a + b*M + c*N for positive exponents 

-a' + b*M + c*N for negative exponents 

IHEXIW 



a 
b 
c 

a' 

L X 



30 



40 



50 



65 



1446 
2125 
1782 

-142 



582 
565 
484 

-31 



128 
174 
147 

-12 



35.2 
45.8 
36.1 



I 75 | 
H 



17.6 
26.9 
21.1 



4.6 | -6.3 
x J 



Error and Exceptional Conditions: 

P : z x = with either REAL(z 2 ) < or 
IMAG(z 2 > * 

O : a. REAL(z 2 *LOG(Zi.) ) > 174.673: error 
caused in IHEEXS or IHEEXL 

b. IHEXXW: 

ABS(IMAG(z a *LOG(Zi))) > 2**18*pi: 
error caused in SIN routine 
(IHESNS) 

IHEXXZ: 

ABSdMAGCza+LOGCzj.))) > 2**50*pi: 

error caused in SIN routine 

(IHESNL) 



IHEXIZ 



30 



40 



50 



65 



75 | 



a | 8524 | 2042 




275 




60.0 | 30.8 | 


b | 5393 | 1397 




289 




59.9 | 34.9 | 


c | 4918 | 1245 




226 




50.2 | 29.1 | 


a'| 95 | -374 




-53 




-5.1 | -4.8 | 













General Floating-Point Exponentiation 



Module Names and Entry Points : 



Argument 



Module 
name 



Entry 
point 



Short float IHEXXW IHEXXW0 
Long float IHEXXZ IHEXXZ0 



Function: 



To calculate Zi.**z 2 , where z x and z 2 are 
complex numbers of the same precision. 
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Implementation : 



Method : 



• Module size: IHEXXW 280 bytes 

IHEXXZ 280 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for System/360 models given 
below are obtained from the table: 

a = IMAG(zj.) b = IMAG(z 2 ) 

C = REAL(Zi) d = REAL(z 2 ) 



If both arguments are non-zero, a call to 
the module IHEAPD is used to shift the 
one with the larger scale factor to give 
it the scale factor of the other, and 
convert it to precision 31. The argu- 
ments are added together, and IHEAPD is 
used to convert the sum to the specified 
precision and to assign it to the target 
field. 

If one of the arguments is zero, the 
other is treated as the sum above. 



r T T T "T T 

| | 30 | 40 | 50 J 65 | 75 

L X X X X X 

IHEXXW 

r T T T T 'T 

I a = j 20606 j 5834 j 1816 j 480 j 291 

|c > | | I I | 

j. x +_„ x + _ + 

|a - | 21750| 6171| 1929 | 509 | 311 
I d < j I | | j 

|. + x_„ x „x x 

| a * j 27448 j 8022 | 2414 | 687 | 417 

|b = | || | | 

j. x x x x x 

| a * j 28263| 8229 | 2576 j 711 j 424 

|b * I I I I I 

l x x x „x x 

IHEXXZ 

r ' T T T ~T ^T 

|a = | 62440| 15325J 3206 j 745| 426 

l c > I I I I I 

j. x x x x _ + 

|a = | 65208| 16062| 3366 | 781| 450 

|d <- | j | | | 

j. + x x -x x 

|a * | 90418| 21611| 4524 | 1056 | 604 
|b = I I I I I 

J. X X X _x _i 

I a * | 92809| 22197 1 4623 | 1077 1 616 
|b * | | | | | 

L X X X X X 



FUNCTIONS WITH REAL ARGUMENTS 



Error and Exceptional Conditions: 

H : FIXEDOVERFLOW or SIZE may occur in 
IHEAPD. 



Implementation : 

• Module size: 216 bytes 

• Execution time: 

Let S ± = time for IHEAPDA with argument 
equal to the sum of the two 
arguments to IHEADD, precision 
31 and scale factor equal to 
the minimum of the scale fac- 
tors of the two arguments 

S 2 = time for IHEAPDB with argument 
equal to the argument to 
IHEADD with the larger scale 
factor 

Then the approximate execution time in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

(i) Both arguments non-zero: 

ti + s ± + s 2 

(ii) At least one argument zero: 

t 2 + S ± 



ADD (Fixed decimal) 

Module Name: IHEADD 
Entry Point: IHEADD0 



r ' — T ■ T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
| t ± | 1895 | 648 | 320 | 85.0 | 65.6 | 
| t 2 | 1619 | 534 | 260 | 71.5 | 53.0 | 

L X X X X X J 



Function: 

ADD (x lf x 2f p,q) where x ± and x a are real 
fixed-point decimal numbers, and (p,q) is 
the required precision of the result. 
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MAX. MIN 



• Execution time: 



Module Names and Entry Points: 



Arqument 


PL/I 
function 


Module 
name 


Entry 
point 


Fixed binary 


MAX 

MIN 


IHEMXB 


IHEMXBX 
IHEMXBN 


Fixed decimal 


MAX 

MIN 


IHEMXD 


IHEMXDX 
IHEMXDN 


Short float 


MAX 

MIN 


IHEMXS 


IHEMXSX 
IHEMXSN 


Long float 


MAX 

MIN 


IHEMXL 


IHEMXLX 

IHEMXLN 



Function: 

To find the maximum or the minimum 
group of arithmetic values. 



of 



have the same base, 



All arguments must 
scale and precision. 

Method: 



IHEMXB, IHEMXS, IHEMXL: The value of the 
current maximum or minimum is set to 
the value of the first argument; it is 
then compared algebraically with the 
next argument and replaced by it if 
appropriate. The process is repeated 
until a test on the argument list 
indicates that all source items have 
been processed, when the current value 
is stored as the result. 

IHEMXD: The address of the current maxi- 
mum or minimum is set to the address of 
the first argument; this argument is 
then compared algebraically with the 
next argument, and the address of the 
latter replaces that of the former if 
appropriate. The process is repeated 
until a test on the argument list 
indicates that all source items have 
been processed,, when the result is 
moved into the target field. 



Implementation : 
• Module sizes: 



IHEMXB 96 bytes 
IHEMXD 120 bytes 
IHEMXS 96 bytes 
IHEMXL 96 bytes 



IHEMXB 

Let N = the number of source arguments 

Then the average execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 
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7 




51 
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186 | 


70 


7 




30 


.0 | 11 


.0 




8.06| 




c 




54 | 


21 


3 




8 


.0 | 2 


.5 




0.9 | 




a 


i | 


356 | 


123 


8 




47 


.3 | 10 


.7 




6.3 | 



























IHEMXD 

Let N = the number of source arguments 

<P#g> = the precision 

L = FLOOR (p/2) + 1 

Then the average execution times in 
microseconds for the System 360 models 
given below are obtained from the fol- 
lowing formulas: 

IHEMXDX 



a + (b + d*L)*N + c* £ <l/i> 

i=2 



IHEMXDN 



a* + (b + d*L)*N + c* £ (1/i) 

i=2 



30 



40 



50 



65 



■t 1 

I 75 | 



a 


670 


217 


84.1 


20.9 


12.3 


b 


232 


92*5 


46.4 


16.4 


13.5 


c 


44 


16.9 


6.5 


1.9 


1.5 


d 


5 


2.5 


2.1 


0.4 


0.4 


a' 


635 


205 


79.6 


19.6 


11.3 
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IHEMXS, IHEMXL 

Let N = the number of source arguments 

Then the average execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 



IHEMXSX,IHEMXLX 

a + b*N + c* 

IHEMXSN,IHEMXLN 

a' + b*N + c* 

IHEMXS 



N 

s 

i=2 



N 

E 

i=2 



(1/i) 



(1/i) 
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1 
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— X- 
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-+- 




-+- 





a 
b 
c 

a' 

L X 
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219 

55 



131.2 
73.3 
21.3 



332 | 118.6 
x X 



44.6 

29.1 

7.3 



11.5 

10.5 

2.5 



40.1 | 10.2 
x . 



7.3 
7.7 
1.9 

6.2 



IHEMXL 




40 



50 



65 



I 75 | 



133.2 


46.4 


11.7 


81.2 


31.4 


10.6 


26.3 


9.3 


2.7 


120.6 


41.9 


10.4 



1 7. 


3 1 


1 7. 


71 1 


1 1- 


9 1 


1 6. 


2 | 



FUNCTIONS WITH COMPLEX ARGUMENTS 



ADD (Fixed decimal) 



Module Name: 



IHEADV 



Entry Point: IHEADV0 

Function: 

ADD(z i# z a , f p,q) where z± and z a are com- 
plex fixed-point decimal numbers, and 
(p,,q) is the required precision of the 
result. 

Method: 

The real parts of each argument are added 
and the sum is assigned to the target 
field by using the real fixed decimal ADD 
module (IHEADD).. The imaginary parts are 
treated similarly. 



Error and Exceptional Conditions : 

H : FIXEDOVERFLOW or SIZE may occur in 
IHEAPD. 



Implementation : 

• Module size: 96 bytes 

• Execution time: 

Let T ± = execution time for IHEADD with 
the real parts of the argu- 
ments for IHEADV as arguments 

T a = execution time for IHEADD with 
the imaginary parts of the 
arguments for IHEADV as argu- 
ments 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

t + Ti + T 2 

r T T T T T 1 

J | 30 | 40 | 50 | 65 | 75 | 
\ — x x x x x H 

j t j 1094 J 396 | 156 | 42.6 j 26.1 | 

L X X X X X J 



MULTIPLY (fixed binary) 

Module Name: IHEMPU 

Entry Point: IHEMPU0 

Function: 

MULTIPLY (z lf z ar p f q) where z, and z a are 
complex fixed- point binary numbers, and 
(p#q> is the required precision of the 
result. 

Method: 

Let the arguments be z ± = a + bl and 
z a = c + dl. 

Then REAL(Zi*z a ) = a*c - b*d 
IMAG(z ± *z a ) = b*c + a*d 

The real and imaginary parts of the 
product are computed. These numbers are 
then shifted to give them the required 
scale factor (q). 

The results of the shifts are tested for 
FIXEDOVERFLOW and truncated by left 
shifts. 
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Error and Exceptional Conditions: 

I : FIXEDOVERFLOW 
Implementation : 

• Module size: 240 bytes 

• Execution time: 

Let q x = scale factor of the 
arguirent 



first 



q 2 = scale factor of the second 
arguirent 

Q = scale factor of the target 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the 
table: 

k = q* + q 2 
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Implementation : 

• Module size: 288 bytes 

• Execution time: 

Let (Piiqa.) (P2,q2> = the precisions of 
the arguments - 

Li = FLOOR ( Pi/2) 

L 2 = FLOOR (p 2 /2) 

T = Time to shift-and-assign the 
result using IHEAPDA, with 
argument of precision 
<Px + Pa + l»<3i + qa> and 
result (p,q) 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

2*T + a + b+Li + c*L 2 + d*Li*L 2 
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MULTIPLY (fixed decimal) 



Module Name: 



IHEMPV 



DIVIDE (fixed binary) 



Entry Point: IHEMPV0 

Function: 

MULTIPLY (Zi i# z 2 ,p f q) where z x and z 2 are 

complex fixed-point decimal numbers, and 

(p,q) is the required precision of the 
result. 

Method: 

Let z x - a + bl and z 2 = c + dl, then: 

REAL(Zi*z 2 ) = a*c - b*d. 
IMAG(Zi*z 2 ) = b*c + a*d. 

The real and imaginary parts are calcu- 
lated and then each is assigned to the 
target with precision (p,q) by separate 
calls to the entry point IHEAPDA of the 
decimal shift and assign module IHEAPD. 

Error and Exceptional Conditions : 

H : FIXEDOVERFLOW or SIZE in IHEAPD. 



Module Name: 
Entry Point: 
Function: 



IHEDVU 



IHEDVU0 



DIVIDE (z 1( z 2 f p»q) where z x and z 2 are 

complex fixed-point binary numbers, and 

(p,q) is the required precision of the 
result. 

Method: 

Let z A = a + bl, and z 2 = c + dl, then: 

REAL(Zi/z 2 ) = (a*c + b*d)/(c**2 + d**2) 
IMAG(Zi/z 2 ) = (b*c - a*d)/(c**2 + d**2) 

The expressions a*c + b*d, b*c - a*d, and 
c**2 + d**2 are computed with a precision 
of 63. The denominator, c**2 + d**2 is 
shifted to precision 31 by either a right 
or left shift. 



32 



Two calls are then made to an incorporat- 
ed subroutine which accepts a numerator 
and shifts it so that it has two insigni- 
ficant leading digits. It then divides 
ty c**2 + d**2 and shifts the quotient to 
the required scale factor (q) . 



Error and Exceptional Conditions: 
I : FIXEDOVERFLOW or ZERODIVIDE 

Implementation : 

• Module size: 408 bytes 

• Execution time: 

Let N A = number of significant bits in 
the expression a*c + b*d 

N 2 = number of significant bits in 
the expression b*c - a*d 

N 3 = number of significant bits in 
the expression c**2 + d**2 

F ± = FLOOR ((61 - N ± )/4) if N ± < 61 
= if N ± > 61 

F 2 = FLOOR ((61 - N 2 )/4) if N 2 ^ 61 
= if N 2 > 61 

F 3 = FLOOR ((N 3 - 32)/4) 

S x = (61 - Nj.) - F ± *4 if N ± < 61 
= if N ± > 61 

S 2 . = (61 - N 2 ) - F 2 *4 if N 2 < 61 
= if N 2 > 61 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

(i) n 3 > 32: ta. + t 2 *(F x + F 2 ) + 
t 3 *(Si + S 2 > + t tt *F 3 

(ii) n 3 = 32: t 5 + t 2 *(F ± + F 2 ) + 
t 3 *(Si + S 2 ) 

(iii) n 3 < 32: t e + t 2 *(F ± + F 2 ) + 
t 3 *(S ± + S 2 ) + t 7 *N 3 
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|t*| 213| 55.6 | 21.5 | 6.0 | 3.8 
|t 5 | 5714| 1780 | 582 | 138 \ 100 
|t 6 | 3712| 1288 | 473 | 115 | 86.5 



|t 7 | -111| -38.8 | -13.3 | -3.4 | -2.6 

L X J. X X X J 



DIVIDE (fixed decimal) 



Module Name: 



IHEDVV 



Entry Point: IHEDW0 

Function: 

DIVIDE (z ± ,z 2 ,p r q) where z ± and z 2 are 
complex fixed- point decimal numbers, and 
(p, f q) is the required precision of the 
result. 

Method: 

Let z ± = a + bl, and z 2 = c + dl, then 

REAL(Zj./z 2 ) = (a*c + b*d)/(c**2 + d**2) 
IMAG(Zi/z 2 ) = (b*c - a*d)/(c**2 + d**2) 

The expressions a*c + b*d, b*c - a*d, and 
c**2 + d**2 are computed. Leading zeros 
are removed from the denominator 
(c**2 + d**2) by truncation on the left 
and a left shift if necessary. If the 
denominator is still more than 15 digits 
long it is truncated on the right to 15 
digits. 

Two calls are then made to an incorporat- 
ed subroutine which accepts a numerator 
and shifts it to precision 31 with 2 
leading zeros by calling IHEAPD (via 
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entry point IHEAPDB) . It then divides by 
c **2 + d**2 and calls IHEAPD (via entry 
point IHEAPDA) to assign the quotient to 
the target field with the required preci- 
sion (p,q) . 

Error and Exceptional Conditions: 

I : ZERODIVIDE 

H : FIXEDOVERFLOW or SIZE in IHEAPD 
Implementation : 

• Module size: 576 bytes 

• Execution time: 

Let Pi,p 2 = the precisions of the argu- 
ments 
L x = FLOOR (pa/2) 

L a = FLOOR (p a / 2) 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formula: 

a + b*L x + c*L 2 + d*Li*L 2 + e*L 2 **2 
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Then ABS(z) is computed as 

X1+SQRTC1 + (Yl/Xl)**2), 

where the fixed binary calculation of 
SQRT(g) for 1 < g < 2 is included within 
the module. 

The first approximation to the square 
root is taken as 

g/(l+g) + (l+g)/4, 

with maximum relative error 1.8*2**-10. 
One Newton- Raphs on iteration gives maxi- 
mum relative error 1.6*2**-20, and suffi- 
ces if XI < 2**(15-q) where q is the 
scale factor of z. 

Otherwise a second iteration is used, 
with theoretical maximum relative error 
of 1.3*2**-40. 



Error and Exceptional Conditions : 

I : FIXEDOVERFLOW 
Implementation : 

• Module size: 184 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for System/360 models given 
below are obtained from the table: 

a = 2**(15-q) 

r T T T T T 1 

I XI I 30 I 40 I 50 I 65 I 75 I 
,. x x 1 x x .j 

| < a j 3809 j 1218 j 320 j 79. 4 j 52.3 j 
^ 1 1 x x x ., 

| > a | 4601 | 1473 | 372 | 93.1 | 59.8 | 
l x x x x x J 



AES (fixed binary) 



Module Name: IHEABU 

Entry Point: IHEABU0 

Function: 

To calculate ABS(z) = SQRT(x**2 + y**2) , 
where z = x + yl. 

Method: 

If x = y, result is x*SQRT(2). Other- 
wise, 

let XI = MAX(ABS(x),ABS(y)) 

Yl = MIN(ABS(X) ,ABS(y)) . 



ABS (fixed decimal) 

Module Name: IHEABV 

Entry Point: IHEABV0 

Function: 

To calculate ABS(z) = SQRT (x**2 + y**2) 
where z = x + yl. 

Method: 

x and y are converted to binary, with 
appropriate scaling if either exceeds 9 
significant decimal digits. 
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Let XI be the maximum,, and Yl the mini- 
mum, of the absolute values of the two 
binary numbers thus obtained. 

Then if XI = Yl = 0, result is 
returned. Otherwise, an approximation to 
ABS(z) is computed as 

X1*SQRT(1 + (Y1/X1)**2), 

where the fixed binary calculation of 
SQRT(g) for 1 < g < 2 is included within 
the module. 



The first approximation to 
root is taken in the form 



the square 



A + B*(l + g) - A/(l + g) 

with maximum relative error 2.17*10**-4, 
and one Newton-Raphson iteration then 
gives a value with maximum relative error 
2.35*10**-8. 

Multiplication by XI produces a value for 
ABS(z) which is rounded and converted to 
decimal, and this suffices if it has not 
more than 7 significant decimal digits. 
Otherwise,, this approximation is scaled 
if necessary and used in a final Newton- 
Raphson iteration for SQRT(x**2 + y**2) 
in decimal, with theoretical maximum 
relative error 2. 76*10**-16. 



Then the approximate execution times in 
microseconds for the System/360 models 
shown below are obtained from the fol- 
lowing formulas: 





L < 5 and 


D 2 


< 7 




a 






L < 5, 7 < 


Di 


< 10 










and 


D 2 


> 7 




b+f 


*L+g*L**2 




5 < L £ 8 and 


D 2 


< 7 




c 






5 < L < 8, 7 < 


D ± 


< 10 










and 


D a 


> 7 




d+f 


*L+g*L**2 




5 < L < 8 














and 10 < 


D± 


< 15 




e+f 


*L+g*L**2 




"T T 


"T~ 




I 




T 




| 30 | 40 
.4. ± 


1 
-X- 


50 


1 


65 


1 75 
X - 




t 1 


T 




1 




t 


a 


| 6220 |1971 


1 


656 




169 


|116 


b 


|13001 |3785 


J1101 




460 


1 352 


c 


| 6666 |2200 


1 


737 




190 


|132 


d 


|13447 |4014 


J1182 




481 


| 368 


e 


j 13918 J4194 


|1279 




509 


| 391 


f 


| 82 | 61.6 


1 


40.1 




7. 


9 | 5.2 


g 


| 56 | 7.5 


1 


0.0 




2. 


| 2.3 




-X x 


_x_ 








x_ 



Error and Exceptional Conditions: 
I : FIXEDOVERFLOW 

Implementation : 

• Module size: 544 bytes 

• Execution times: 

Let (p, q) = the precision of the argu- 
ment 

L = CEIL<(p+l)/2)„ i.e., the 
length in bytes of each of the 
real and imaginary parts of 
the argument 

Dj. = maximum number of significant 
digits in real and imaginary 
parts of the argument 

D a = number of significant digits 
in result 



ABS (floating-point) 



Module Names and Entry Points : 



Argument 



Module 
name 



Entry 
point 



Short float IHEABW IHEABW0 
Long float IHEABZ IHEABZ0 



Function: 



To calculate ABS(z) = SQRT (x**2 + y**2), 
where z = x + yl. 



Method: 

Let z=x+yl. If x = y = 0, answer is 
0. 

Otherwise let Zl = MAX (ABS(x) ,ABS(y) ) 
and Z2 = MIN (ABS(x) ,ABS(y) ) . 

Then the answer is computed as 

Z1*SQRT(1 + (Z2/Z1)**2). 
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Accuracy : 

IHEABW 

r t 11 

| Arguments | Relative Error | 
| | *io**6 | 
j. T + T .j 

| Range j Distribution! RMS j Maximum | 
j. + + + ^ 

| Full range | Exponential | 0.833 | 2.02 | 
| | radially, | | | 

| | uniform | | | 

| | round origin j j j 

l x x x J 

IHEABZ 

r t 1 

| Arguments | Relative Error | 
| | *10**15 j 

h T + T i 

| Range | Distribution! RMS | Maximum | 
,. x x x ^ 

| Full range | Exponential | 0.828 | 3.38 | 
| | radially, | | | 
| | uniform j j | 
| | round origin j j j 
l x x x J 



Error and Exceptional Conditions: 
I : OVERFLOW 



Implementation : 

• Module size: IHEABW 128 bytes 

IHEABZ 128 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

| Module | 30 | 40 | 50 | 65 | 75 | 
I Name I 1 I I I I 

j. X x 1 x x .| 

| IHEABW j 5595| 1493 j 447 j 129 j 79. 9 j 

| IHEABZ | 14318| 3191 | 695 | 174 | 104 | 

L X X X X X J 
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CHAPTER 3: MATHEMATICAL FUNCTIONS 



The Library supports all float arithmet- 
ic generic functions and has separate 
modules for short and long precision real 
arguments and also for short and long 
precision complex arguments where these are 
admissible. 



Speed 



The average execution times given are 
based on the IBM System/360 Instruction 
Timing Information , Form A22-6825. These 
times include times for the modules called. 



Since the calling sequence generated in 
compiled code is the same as that required 
for passing the same arguments to a PL/I 
procedure, it is permissible to pass the 
names of any of the float arithmetic gener- 
ic functions as arguments between proce- 
dures, according to the normal rules for 
entry names. 

Any restrictions on the admissibility of 
arguments are noted under the headings 
'Range' and 'Error and Exceptional Condi- 
tions. ' 

Range ; This states any ranges of arguments 
which a module assumes to have been exclud- 
ed prior to its being called. 

Error and Exceptional Conditions ; These 
cover conditions which may result from the 
use of a routine; they are listed in four 
categories: 

P — Programmed conditions in the module 
concerned. Programmed tests are 
made where this is not too costly 
and, if an invalid argument is 
found, a branch is taken to the 
entry point IHEERRC of the execution 
error package (EXEP) . This results 
in the printing of an appropriate 
message and in the ERROR condition 
being raised. 

I — Interrupt conditions in the module 
concerned. For those routines where 
SIZE and FIXEDOVERFLOW are detected 
by programmed tests or where hard- 
ware interruptions may occur, the 
OVERFLOW and UNDERFLOW conditions 
pass to the ON handler (IHEERR) and 
are treated in the normal way. The 
machine is assumed to be enabled for 
all interruptions except signifi- 
cance, which is masked off. 

O — Programmed conditions in modules 
called by the module concerned. 
These occur when invalid arguments 
are detected in the module called. 

H — As I, but the interrupt conditions 
occur in the modules called by the 
module concerned. 



Accuracy 



In order to appreciate properly the 
meaning of the statistics for accuracy 
given with each module, some consideration 
of the limits and implications of these 
statistics is required. Because the size 
of a machine word is limited, small errors 
may be generated by mathematical routines. 
In an elaborate computation, slight inaccu- 
racies can accumulate and become large 
errors. Thus, in interpreting final 
results, errors introduced during the var- 
ious intermediate stages must be taken into 
account. 

The accuracy of an answer produced by a 
routine is influenced by two factors: (1) 
the accuracy of the argument and (2) the 
performance of the routine. 

Most arguments contain errors. ,,„An error 
in a given argument may have accumulated 
over several steps prior to the use of the 
routine. Even data fresh from input con- 
version may contain slight errors. The 
effect of an argument error on the accuracy 
of an answer depends solely on the nature 
of the mathematical function involved and 
not on the particular coding by which that 
function is computed within a routine. In 
order to assist users in assessing the 
accumulation of errors,, a guide on the 
propagational effect of argument errors is 
provided for each function. Wherever pos- 
sible, this is expressed as a simple formu- 
la. 

The performance statistics supplied in 
this document are based upon the assumption 
that the arguments are perfect (i.e., with- 
out errors, and therefore having no argu- 
ment error propagation effect upon 
answers) . Thus the only errors in answers 
are those introduced by the routines them- 
selves. 

For each routine, accuracy figures are 
given for the valid argument range or for 
representative segments of this. In each 
case the particular statistics given are 
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those most meaningful to the function and 
range under consideration. 

For example, the maximum relative error 
and the root-mean-square of the relative 
error of a set of answers are generally 
useful and revealing statistics, but are 
useless for the range of a function where 
its value tecomes 0, since the slightest 
error of the argument value can cause an 
unbounded fluctuation in the relative mag- 
nitude of the answer. Such is the case 
with SIN(x) for values of x close to pi; in 
this range it is more appropriate to dis- 
cuss absolute errors. 



errors range from 2**-24 to 2**-20 
(5.96*10**-8 to 9.5*10**-7). Assuming 
exact operands , a product or quotient is 
correct to the 24th binary digit of the 
mantissa. Hence truncation errors contrib- 
uted by multiplication or division are no 
more than 2**- 20. The same is true for the 
sum of two operands of the same sign. 
Subtraction, on the other hand, is the 
commonest cause of loss of significant 
digits in any number system. For short- 
precision operations, therefore, a guard 
digit is provided which helps to reduce 
such loss. 



The results were derived from random 
distributions of 5000 arguments per seg- 
ment, generated to be either uniform or 
exponential, as appropriate. It must be 
emphasized that each value quoted for the 
maximum error refers to a particular test 
using the method described above, and 
should be treated only as a guide to the 
true maximum error. 

This explains, for example, why it is 
possible that the maximum error quoted for 
a segment may be greater than that found 
from a distribution of different arguments 
over a larger range which includes the 
former. 



Hexadecimal Truncation Errors 



While the use of hexadecimal numbers in 
System/360 has led to increased efficiency 
and flexibility, the effect of the variable 
number of significant digits carried by the 
floating-point registers must be noted in 
making allowance for truncation errors. In 
the production of the PL/I Library, special 
care was taken to minimize such errors, 
whenever this could be accomplished at 
minor cost. As a result, the relative 
errors produced by some of the Library 
routines may be considerably smaller than 
the relative error produced in some instan- 
ces by a single operation such as multi- 
plication. 

Representations of finite length entail 
truncation errors in any number system. 
With binary normalization, the effect of 
truncation is roughly uniform. With hexa- 
decimal normalization, however, the effect 
varies by a factor of 16 depending on the 
size of the mantissa; in a chain of compu- 
tations^, the worst error committed in the 
chain usually prevails at the end. 

In short-precision representation, a 
number has between 21 and 24 significant 
binary digits. Therefore, the truncation 



In long-precision representation, a 
number has between 53 and 56 significant 
binary digits. Therefore truncation errors 
range from 2**-56 to 2**-52 (1. 39*10**-17 
to 2. 22*10**-16) . Assuming exact operands, 
a quotient is correct to the 56th binary 
digit of the mantissa. Therefore, trunca- 
tion errors resulting from division are no 
more than 2**- 52. The accuracy of a pro- 
duct, on the other hand, depends on the 
necessity for post-normalization. If the 
mantissas of both operands are close to 1, 
the truncation error of a product is about 
2**-56. If the product of the mantissas is 
about 1/16, the truncation error is about 
2**-52. On the other hand, if the mantis- 
sas of both operands are close to 1/16, the 
intermediate product has 7 leading zeros, 
and post-normalization introduces 4 trail- 
ing zeros. In this case, the truncation 
error can be close to 2**-48 
(3.55*10**-15) . In particular, multi- 
plication by 1 in the long-precision form 
has the effect of erasing the last hexa- 
decimal digit of the multiplicand. 



Normal care in numerical analysis should 
be exercised for addition and subtraction. 
In particular, when two algorithms are 
theoretically equivalent, it usually pays 
to choose the one which avoids subtraction 
between operands of similar size. There is 
no guard digit for long- precision additions 
and subtractions. 



Hexadecimal Constants 



Many of the modules described below 
discriminate between algorithms or test for 
errors by comparisons involving hexadecimal 
constants; it must be realized that where 
decimal fractions are used in the descrip- 
tions the fractions are only quoted as 
convenient approximations to the hexadeci- 
mal values actually employed. 
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Terminology 

Maximum and root-mean-square values for 
the relative and (where necessary) the 
absolute errors are given for each module. 
These are defined thus: 

Let f(x) = the correct value for a func- 
tion 

g(x) = the result obtained from the 
module in question 

Then the absolute error of the result is 

ABS(f (x) - g(x)) , 
and the relative error of the result is 

ABSUf (x) - g(x))/f (x)). 

Let the number of sample results obtained 
be N; then the root-mean-square of the 
absolute error is 



SQRT(X)i(ABS(f ix\_) - g(xi))**2)/N), 

and the root-mean-square of the relative 
error is 

SQRT(]£i(ABS((f (Xj.) - g(XjL))/f (Xi))**2)/N) . 



The Library mathematical modules are 
summarized in Figures 4 and 5. 
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Figure 4. Mathematical Functions with Real 
Arguments 



Function 



j SQRT 

| EXP 

| LOG 

j SIN, COS , SINH, COSH 

| TAN, TANH 

j ATAN, ATANH 

L 



Complex Arguments 



Short Float | Long Float | 

H 



IHESQW 
IHEEXW 
IHELNW 
IHESNW 
IHETNW 
IHEATW 



IHESQZ 
IHEEXZ 
IHELNZ 
IHESNZ 
IHETNZ 
IHEATZ 



Figure 5. Mathematical Functions with Com- 
plex Arguments 



FUNCTIONS WITH REAL ARGUMENTS 

SQRT (short floating-point real) 

Module Name: IHESQS 
Entry Point: IHESQSO 
Function: 

To calculate the square root of x. 
Method: 

If x = 0, SQRT(x) = 0. Otherwise, let 

X = 16**(2*p + q)*f, 

where p is an integer, q = or 1, and 
1/16 < f < 1. Then 

SQRT(X) = 16** (p + q)*z, 
where z = SQRT(f) if q = 0, 
z = SQRT(f)/4 if q = 1. 

An initial approximation, y # is taken in 
the hyperbolic form a + b/(c + f ) with 
different sets of constants for the two 
cases: 

1. q = a = 1.80713 

b = -1.57727 
c = 0.954182 

The maximum relative error in this range 
is then less than 2** (-5. 44), with an 
exact fit at f = 1 to guard as far as 
possible against loss of the last hexa- 
decimal digit when f is nearly 1. 

2. q = 1 a = 0.428795 

b = -0.3430368 
c = 0.877552 

The maximum relative error in this range 
is less than 2** (-6) *f **(-l/8) . 
Then y x = 16** (p + q)*y . 
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Two Newton-Raphson iterations then yield: 

y 2 = <Yx + x/ yi )/2 
SQRT(x) = y 2 + (x/y 2 - y 2 )/2 

For case q = f the final relative error 
from this algorithm is less than 
2** (-24.7), and, for case q = 1, less 
than 2** (-29) . 

Effect of Argument Error: 

The relative error caused in the result 
is approximately half the relative error 
in the argument. 

Accuracy: 

r t t 

| Arguments | Relative Error | 
j j *10**6 | 

j. T + T ^ 

| Range | Distribution | RMS | Maximum | 
j. + + + .j 

| Full Range | Exponential | 0.230 j 0.924 j 

l X X X J 

Error and Exceptional Conditions: 

P : x < 
Implementation : 

• Module size: 168 bytes 

• Execution time: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

r T T T T T 1 

| Module j 30 j 40 j 50 j 65 j 75 j 

I Name I I I I I | 
|. x x x x x .| 

|IHESQS j 3140 j 793 j 227 j 68.4 | 40.7 j 

L X X X X X J 



SQRT (long floating-point real) 



An initial approximation, y , is taken by 
using (2/9 + 8/9*f) for SQRT(f). Multi- 
plication by 2** (-2) when q = 1 is accom- 
plished by using the HALVE instruction 
twice. The maximum relative error of 
this approximation is 1/9 . 

Four Newton-Raphson iterations of the 
form yn+i = (yn + x/y n )/2 are then 
applied, two in short precision and two 
in long precision,., the last being comput- 
ed as 

SQRT(x) = y 3 + (x/y 3 - v 3 )/2 

to minimize the truncation error. 

The maximum relative error in the result 
from this algorithm is 2** (-65. 7). 

Effect of an Argument Error: 

The relative error caused in the result 
is approximately half of the relative 
error in the argument. 

Accuracy: 



T 1 

Arguments | Relative Error 
I *10**15 



Range (Distribution 



RMS 



Maximum 



x > 10** (Exponential j 0.0276 j 
(-52) | | | 

X X X. 



0.124 



Error and Exceptional Conditions: 

P : x < 
Implementation : 

• Module size: 160 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 



Module Name: IHESQL 

Entry Point: IHESQL0 

Function: 

To calculate the square root of x. 

Method: 

If x = 0, SQRT(x) = 0. Otherwise, let 
x = 16**(2*p - q)*f, where p is an inte- 
ger, q = or 1, and 1/16 < f < 1. Then 

SQRT(X) = 16**p*2**(-2*q)*SQRT(f ) . 



T T T 1 

40 I 50 I 65 I 75 I 



r t -r 

j Module j 30 j 

I Name | | j j j j 

|. x 1 x x x -i 

j IHESQL j 8282 | 1733 | 376 | 97. 7 j 57.2 j 
L X X X X X J 



EXP (short floating-point real) 

Module Name: IHEEXS 
Entry Point: IHEEXS0 
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Function: To calculate e to the power x. 

Method: 

If x < -180.2, a zero result is returned 
immediately. 

Otherwise, EXP(x) is calculated as: 

2**(x*LOG2(e)) 

The calculation is performed as follows: 

x*LOG2(e) = f + N f 

where N = 4h + g, h is an integer such 
that g = 0, 1, 2 or 3, and < f < 1.0. 

Then, by subtracting 0.5, this is reduced 
to the range -0.5 < f < 0.5. Next, 2**f 
is calculated as: 

(a + b/(c + x*x)+x)/(a + b/(c + x*x)-x). EXP (long floating-point real) 



Implementation : 

• Module size: 232 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

| Module | 30 | 40 | 50 | 65 | 75 | 
I Name | I I I I I 

|. + x x + + -I 

JIHEEXS j 3847 | 1172 | 356 | 90. 0| 58.0 j 

L X X X X _X J 



This is multiplied by 2**0.5 and then 
shifted in the appropriate direction to 
give the effect of multiplication by 
2**g. Finally, the exponent of the 
result is obtained from h. 



Effect of Argument Error: 

The relative error caused in the result 
is approximately equal to the absolute 
error in the argument, i.e., to the 
argument relative error multiplied by x. 
Thus for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 



Accuracy: 



Arguments 



. T , ., 

| Relative Error 
I *io**6 



Range 



■T T T 

| Distribution | RMS | Maximum 



-1 < x < lj Uniform j 0.132 j 0.490 
x x x ;_. 

Full Range | Uniform j 1.29 | 2.61 
x . x . x . 



Error and Exceptional Conditions: 
I : OVERFLOW if x > 174.673 



Module Name: IHEEXL 

Entry Point: IHEEXL0 

Function: To calculate e to the power x. 

Method: 

If x < -180.2183, return zero result. 

Otherwise let y = x/LOG(2) 

= 4*a - b - c/16 - d 

where a, b and c are integers, < b < 3, 
< c < 15 and < d < 1/16. 

Then EXP(x) = 2**y 

= 16**a*2**(-b)*2**(-c/16) 
*2**(-d) . 

Compute 2**(-d) by using the Chebyshev 
interpolation polynomial of degree 6 over 
the range < d < 1/16, with maximum 
relative error 2** (-57). 

If c > 0, multiply 2** (-d) by 2** (-c/16). 
The constants 2** (-c/16), 1 < c < 15, are 
included in the subroutine. 

If b > 0, halve the result b times. 

Finally,, multiply by 16**a by adding a to 
the characteristic of the result. 

Effect of an Argument Error: 

The relative error caused in the result 
is approximately equal to the absolute 
error in the argument, i.e., to the 
argument relative error multiplied by x. 
Thus for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 
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Accuracy: 



j Arguments 


— T - 

1 
1 


Relative 
♦10** 


Error j 
15 | 


[._ _ T -. T _ 

| Range | Distribution! 
L _ J. _ ± 


T~ 

RMS j 

i_ 


Maximum j 


r 1 

j-1 < x < lj Uniform 

L X 


T 

1 
L_ 


0.0674| 

L_ 


0.216 | 


r T 

|Full range | Uniform 

L X . 


— T 
1 


0.867 | 


2.30 | 



Error and Exceptional Conditions: 

I : OVERFLOW if x > 174.673 
Implementation : 

• Module size: 448 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 



r t- 

j Module | 
I Name I 



30 



40 



50 



65 



h 



■T 1 

I 75 | 
I 



|IHEEXL | 12131| 2901 | 616 | 343 | 194 

L J. X J. J. X J 



LOG f LOG2, ICG10 (short f loatinq-point 
real) 



Module Name: 



Entry Points: 



IHELNS 



Mathematical 
function 

Log x to the base e 
Log x to the base 2 
Log x to the base 10 



PL/I 
name 

LOG(x) 

L0G2(x) 

LOGlO(x) 



Entry 
point 

IHELNSE 
IHELNS2 
IHELNSD 



Function: To calculate log x. 

Method: 

Let x = m*16**p where 1/16 < m < 1 and p 
is an integer. 

Two constants, a (= base point) and b (= 
-LOG2(a)), are defined as follows: 

1/16 < m < 1/8, a = 1/16, b = 4; 
1/8 < m < 1/2, a = 1/4, b = 2; 
1/2 < m < 1, a = 1, b = 0. 

Let y = (m-a)/(m+a) 



Then m = a*(l+y)/(l-y) and ABS(y) < 1/3. 

Now x = 2**(4*p-b)*(l+y)/(l-y) . 

Therefore 

LOG(x) = (4*p-b)*LOG(2) + 
LOG((l+y)/(l-y)) . 

LOG( (l+y)/(l-y) ) is computed by using the 
Chebyshev interpolation polynomial of 
degree 4 in y**2 for the range 
< y**2 < 1/9, with maximum relative 
error 2**-27.8. LOG2 (x) or LOG10(x) are 
calculated by multiplying LOG(x) by 
LOG2(e) or LOGlO(e) respectively. 



Effect of Argument Error: 

The absolute error caused 
is approximately equal 
error in the argument, 
argument is close to 1, even the round- 
off error of the argument causes a 
substantial relative error in the answer, 
since the function value there is very 
small. 



in the result 

to the relative 

Thus if the 



-I Accuracy: 









~i 


1 
1 


Arguments 


| Relative Error 
1 *10**6 


1 
1 


H- 




x T 


H 



Range 



| Distribution | RMS | Maximum 
-X x x 



IHELNSE 

r t t t 

I Excluding | | | 

| 0.5 < x (Exponential j 0.0320 j 0.577 

l< 2.0 | || 

L X X X 



IHELNS 2 



r t t T 

I Excluding | | | 

j 0.5 < x | Exponential j 0.342 | 
l< 2.0 | || 



0.754 



IHELNSD 

r t t t 

I Excluding | | | 

j 0.5 < x (Exponential | 0.170 j 0.963 

l< 2.0 | || 

L X X X 
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r t t 

| Arguments | Absolute Error | 
I | *10**6 | 
|. T + T ^ 

| Range | Distribution! RMS | Maximum | 

L X JL . X _J 

IHELNSE 

| 0.5 < x | Uniform | 0.0960 | 0.394 j 
l< 2.0 | III 

L X X X J 

IHELNS2 

|0.5 < x | Uniform | 0.177 1 0.842 | 
l< 2.0 | | | | 

L X X - X J 

IHELNSD 

| 0. 5 < x | Uniform | 0.0526| 0.164 | 
l< 2.0 | || | 

L X X X J 



Error and Exceptional Conditions: 
P : x < 

Implementation : 

• Module size: 256 bytes 

• Execution times: 



Method: 

Let x = 16**p*2**(-q)*m where p is the 
exponent, q is an integer such that 
< q < 3, and 1/2 < m < 1. 

Two constants, a (= base point) and b (= 
-LOG2(a)), are defined as follows: 

1/2 < m < 1/SQRT(2): a = 1/2, b = 1 
1/SQRT(2) < m < 1: a=l, b=0 

Let y = (m - a)/(m + a). 

Then m = a*(l + y)/(l - y) and ABS(y) 
< 0.1716. 

Now x = 2**(4*p - q - b)*(l + y)/(l - y) 

Therefore 

LOG(x) = (4*p - q - b)*LOG(2) + 
LOG((l + y)/(l - y)) . 

LOG((l + y)/(l - y)) is computed by using 
the Chebyshev interpolation polynomial of 
degree 7 in y**2 for the range 
< y**2 < 0.02944, with maximum relative 
error 2**(-59.6) . 

LOG2(x) or LOGlO(x) is calculated by 
multiplying the result by LOG2(e) or 
LOG10(e) respectively. 



Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 



30 



r T - 

| Entry | 
j Point | 

I- + 

IHELNSE | 4669 



IHELNS2I 5041 



40 



1238 



1342 



417 



|: 

h 

| IHELNSD | 5054 | 1366 | 417 | 180 | 101.3 

L X X— X X X J 



50 



385 



65 | 75 



, + .| 

173 | 95.7 | 



180 | 101.3| 



Effect of an Argument Error: 

The absolute error caused in the result 
is approximately equal to the relative 
error in the argument. Thus if the 
argument is close to 1 , even the round- 
off error of the argument causes a 
substantial relative error in the answer, 
since the function value there is very 
small. 



LOG,, LOG2, LOG10 (long floating-point real) 



Module Name; 



IHELNL 



Entry Points: 



Mathematical 
Function 



PL/I 
name 



Entry 
point 



Log x to the base e LOG(x) IHELNLE 

Log x to the base 2 LOG2(x) IHELNL2 

Log x to the base 10 LOGlO(x) IHELNLD 

Function: To calculate log x. 
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Accuracy: 



| Arguments 
L 


1 
1 


Relative Error | 
♦10**15 | 


r — 

| Range 

L _^ 


T T 

| Distribution! 

X X_ 


RMS 


i 

1 
i 


1 

Maximum | 


IHELN1E 


r ------ 

(Excluding 
| . 5 < x 
|< 2.0 

L _ 


T 
1 

(Exponential 
1 

X 


i 

1 

1 

1 
i 


T 

1 

0.0530 j 
1 

X 


0.329 | 


IHELN12 


r ^ 

| Excluding 
| 0. 5 < x 

|< 2.0 

i . 


T 
1 

(Exponential 
1 

X 


i 

1 

1 

1 
i 


0.443 


i 

1 

1 

1 
i 


2.60 | 


IHELNLD 


r~ ~ ~^ 

| Excluding 
| 0. 5 < x 
|< 2.0 

L_ 


T 

1 

(Exponential 

1 

X 


i 

1 

1 

1 
i 


0.155 


i 

1 

1 

1 
i 


0.402 | 




r~ ~ ~ *" 

| Arguirients 

L ._ 


i 

1 
1 


Absolute Error j 
♦10**15 j 


r 

| Range 
i 


T T 

| Distribution! 

X X 


RMS 


i 

1 
i 


Maximum j 


IHELN1E 


r " 

|0. 5 < x 
|< 2.0 

i _ 


T 

| Uniform 
1 

X 


i 

1 

1 
i 


0.192 


1 

1 
i 


0.507 | 


IHELN12 


r *~ 

|0.5 < x 
|< 2.0 

L 


T 

| Uniform 
1 

X 


T 

1 

1 
1 


0.245 


T 

1 

1 
1 


0.466 | 


IHELNLD 


r *~ 

|0.5 < x 
|< 2.0 


T 

| Uniform 
1 


1 

1 
1 


0.0318| 
1 


0.0625| 



Error and Exceptional Conditions: 
P : x < 



Implementation : 

• Module size: 360 bytes 

• Execution times : 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 



r — t- 

| Entry | 


T" 

30 j 


40 


"T" 


50 


-T" 


65 


- T" 
1 


75 | 


| Point| 


1 












1 




1- +- 

j IHELNLE j 


x. 

16216| 


3926 


-+- 


788 


"+- 


178 


1 


98. 5| 


1- +" 

j IHELNL2 j 


+ . 

17315| 


4196 


-+- 


834 


-+- 


190 


1 


107 | 


1- +" 

j IHELNLD j 


x. 

17284| 


4192 


-+- 


828 


-+" 


188 


-+- 
1 


105 | 


L J._ 


x_ 




1 




_1. 




_X- 





SIN, SIND, COS, COSD (short floating-point 
real) 



Module Name: IHESNS 
Entry Points: 



Mathematical 
function 


PL/I 
name 




Entry 
point 


Sin(x radians) 
Sin(x degrees) 
Cos(x radians) 
Cos(x degrees) 


SIN(x) 
SIND(x) 
COS(x) 
COSD(X) 




IHESNSS 
IHESNSZ 
IHESNSC 
IHESNSK 


Function: To calculate 


sin x or 


cos 


x. 


Method : 








Let k = pi/4 








Evaluate p = ABS(x)* 


(1/k) if 


x 


is in 



radians 
or p = ABS(x)*(l/45) if x is in 

degrees, 
using long- precision multiplication to 
safeguard accuracy. 

Separate p into integer part q and frac- 
tional part r, i.e. , p=q+r where 
£ r < 1. 

Define q± = q if SIN or SIND is required 

and x > 0; 
q ± = q + 2 if COS or COSD is 

required; 
qx = <3 + 4 if SIN or SIND is 

required and x < 0. 

Then for all values of x each case has 
been reduced to the computation of 
SIN(k*(qi+r) ) = SIN(t) say, where t £ 0. 
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Let q a 

if q 3 



If 
If 
If 
If 
If 
If 
If 



qa 

<3a 

qa 
qa 
qa 
qa 
qa 



= M0D(q 1# 8). 
= 0, SIN(t) 
= 1, 



= 2, 
= 3, 
= U f 
= 5 f 
= 6, 
= 7, 



SIN(t) 
SIN(t) 
SIN(t) 
SIN(t) 
SIN(t) 
SIN(t) 
SIN(t) 



SIN(k*r) 

COS(k*(l-r)) 

COS(k*r) 

SIN(k*(l-r)) 
-SIN(k*r) 
-COS(k*(l-r)) 
-COS(k*r) 
-SIN(k*(l-r)) 



Thus it is necessary to compute only 
SINCk+ra.) or COS(k*rj.) where x x = r or 
1 - r and < r x < 1. 

This is performed by using the Chebyshev 
interpolation polynomials of degree 3 in 
x ± **2, with maximum relative error of 
2**-28.1 in the sine polynomial and 
2**-24.6 in the cosine polynomial. 



r t 1 

| Arguments | Relative Error | 
| | *10**6 | 
y T + — T .] 

| Range (Distribution | RMS | Maximum | 
l J. x x J 

IHESNSS 

r t t t 1 

|ABS(x) < I Uniform | 0.198 j 1.40 j 

I Pi/2 I III 

L X X X J 



Error and Exceptional Conditions: 



P : IHESNSS, IHESNSC: 

ABS(x) £ 2**18*pi 

IHESNSZ, IHESNSK: 

ABS(x) > 2**18*180 



Effect of an Argument Error: 

The absolute error of the answer is 
approximately equal to the absolute error 
in the argument. Hence, the larger the 
argument, the larger its absolute error 
and the larger the absolute error of the 
result. Since the function diminishes 
periodically for both sine and cosine, no 
consistent ccntrol of the relative error 
can be maintained outside the range -pi/2 
to pi/2 radians (or -90 to +90 degrees) . 

Accuracy: 

r t 1 

| Arguments | Absolute Error | 
| | *10**6 I 
|. T 1 T -J 

| Range | Distribution! RMS | Maximum | 

L _ X X X J 

IHESNSS 

|ABS(x) < | Uniform | 0.0557| 0.126 | 

I Pi/2 | III 
|. x x x .] 

| pi/2 < | || | 

|ABS(x) j Uniform j 0.0553 j 0.148 j 
|< 10 | || | 

y x + x 4 

|10 < | || | 

|ABS(x) j Uniform j 0.0560 j 0.143 j 

| £ 100 | || | 

L X X X J 

IHESNSC 

r T T T 1 

JO < x < pi j Uniform j 0.0553 j 0.149 j 
|. x x x \ 

|-10 < x | | | | 

|< 0, | Uniform j 0.0571J 0.154 j 
iFi < x | | | | 

|< 10 | | | | 

j. x „ x x ^ 

110 < 1 || | 

|ABS(x) | Uniform | 0.0553| 0.142 | 

| < 100 | | | | 

L . X . X , X J 



Implementation : 

• Module size: 320 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 

|ABS(x) | 30 | 40 | 50 | 65 | 75 | 

L X X X X X J 

IHESNSS 



[< Pi/** I 



I* Fi/4 I 

L X. 



4091 
4386 



■T T" 

| 1120 | 



| 1190 | 
-X x. 



333 
362 



■T T 1 

j 85. 0J 50.6 



| 92.5| 53.8 

.X X 



IHESNSC 



|< pi/4 | 4078 | 1115 | 329 | 83. 6| 49.9 
|£ pi/4 | 4373 | 1184 | 357 | 91. 0| 53.1 

L X X X X X 



IHESNSZ 



|< 45 | 4026 | 1132 | 338 | 86.3 | 51.6 
|£ 45 | 4421 | 1202 | 366 | 93.7 | 54.8 

L _X X X X X 



IHESNSK 



|< 45 | 3693 | 1127 | 334 | 84. 8| 51.0 
|> 45 | 4408 | 1196 | 362 | 92. 3| 54.1 

L X X X X X 
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SIN,, SIND, COS, COSD (long floating-point 
real) 



Effect of an Argument Error: 



Module Name: 



Entry Points: 



IHESNL 



Mathematical 
function 

Sin(x radians) 
Sin(x degrees) 
Cos (x radians) 
Cos(x degrees) 



PL/I 
name 

SIN(x) 
SIND(x) 
COS(x) 
COSD(x) 



Entry 
point 

IHESNLS 
IHESNLZ 
IHESNLC 
IHESNLK 



Function: To calculate sin x or cos x. 

Method: 

Let y = ABS(x)/(pi/4) for x in radians, 
or y = AES(x)/45 for x in degrees, 
and y = q + r, q integral, < r < 1. 

Take q ± = q for SIN or SIND with positive 
or zero argument, 
q*. = q + 2 for COS or COSD, 
q x = q + t\ for SIN or SIND with 
negative argument,, 
and q 2 = M0D(qi.,8) . 

Since COS(x) = SIN(ABS(x) + pi/2) 
and SIN(-x) = SIN(ABS(x) + pi), 

it is only necessary to find 

SIN (pi/4* (q 2 + r)), for < q a < 7 . 

Therefore compute: 



SIN(pi/4*r) lf 

COS (pi/4* (1 - r)) , 

C0S(pi/4*r)„ 

SIN (pi/4* (1 - r)) 



if q 2 = or 4,, 

if q 2 = 1 or 5, 

if q a = 2 or 6, 

if qa = 3 or 7. 



SIN(pi/4*ri)/r ± , where r^ is r or 
(1 - r) , is computed by using the Cheb- 
yshev interpolation polynomial of degree 
6 in r x **2, in the range < r ± **2 £ 1, 
with maximum relative error 2** (-58). 

C0S(pi/4*ri) is computed by using the 
Chefcyshev interpolation polynomial of 
degree 7 in r 1 **2, in the range 
< r ± **2 < 1„ with maximum relative 
error 2** (-64. 3) . 



Finally, if q a £ 4 a 
given to the result. 



negative sign is 



The absolute error 
approximately equal to 
in the argument. Henc 
argument, the larger 
and the larger the abs 
result. Since the 
periodically for both 
consistent control of 
can be maintained outs 
to pi/2 radians (or -9 



of the answer is 
the absolute error 
e, the larger the 
its absolute error 
olute error of the 
function diminishes 
sine and cosine, no 
the relative error 
ide the range -pi/2 
to +90 degrees) . 



Accuracy: 



IHESNLS 



r T _ 

| Arguments | 
1 1 

L , X. 


Relative Error j 
♦10**15 | 


| Range (Distribution | 

L J. ±_ 


RMS | Maximum | 
i j 


r 1 1 

|-pi/2 < x | Uniform | 

l< Fi/2 | I 

L J. X~ 


0.0542| 0.381 | 
1 1 

_L J 


IHESNLC 


r t~ 

| Arguments | 

1 1 

L X- 


Absolute Error j 
♦10**15 j 


r t t 

j Range (Distribution | 

l I L. 


t 1 

RMS | Maximum | 
1 | 


r 1 1 

|-pi/2 < x | Uniform | 

l< Fi/2 | | 


0.0604| 0.168 | 
1 1 



Error and Exceptional Conditions: 

P : IHESNLS, IHESNLC: 

ABS(x) > 2'**50*pi 

IHESNLZ, IHESNLK: 
ABS(x) > 2**50*180 

Implementation : 

• Module size: 416 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 



Entry j 
Point | 



30 



40 



50 



IHESNLS j 13654 j 3290 | 661 
-_x 1 x 

IHESNLC| 13641| 3288 | 654 
IHESNLZ | 13689| 3302 | 665 



— +- 



IHESNLK j 13676 j 3300 j 659 
x x x 



65 | 75 



155 | 85.3 
x 

153 | 84.3 
157 | 86.3 



155 | 85.3 
x 
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TAN, TAND (short float inq-point real ) 



of pi/2, an argument error will produce a 
large relative error in the result. 



Module Name: 
Entry Points: 



IHETNS 



Mathematical 
function 

Tan(x radians) 
Tan(x degrees) 



PL/I 
name 

TAN(x) 
TAND(x) 



Entry 
point 

IHETNSR 
IHETNSD 



Function: To calculate tan x. 
Method: 



Evaluate p = 


(4/pi)*ABS(x) 
radians,, 


if 


X 


is 


or p = 


(l/45)*ABS(x) 
degrees,, 


if 


X 


is 



in 



using long-precision multiplication to 
safeguard accuracy. 

Let q and r be respectively the integral 
and fractional parts of p. 

If q is even, put s = r; 
if q is odd, put s = 1-r. 

Let q ± = M0D(q,4). Then 

If q ± = 0, TAN(ABS(x)) = TAN(pi*S/4) 

If q ± = 1, TAN(ABS(x)) = COT(pi*s/4) 

If q± = 2, TAN(ABS(x)) = -COT(pi*S/4) 

If qi = 3, TAN(ABS(x)) = -TAN (pi*S/4) 

Compute TAN(pi*s/4) and COT(pi*s/4) as 
the ratio of two polynomials: 

TAN(pi*s/4) = s*P(s** v 2)/Q(s**2) 
CCT(pi*s/4) = Q(s**2)V(s*P(s**2)) 

where P(s**2) = 212.58037 - 12.559912 
*s**2 
Q(s**2) = 270.665736 - 71.645273 
*s**2 + s**4. 

Finally, if x < 0,, put 

TAN(x) = -TAN(ABS(x)) . 

Effect of an Argument Error: 

The - absolute error of the answer is 
approximately equal to the absolute error 
of the argument multiplied by 
(1 + TAN(x)**2). Hence if x is near an 
odd multiple of pi/2, an argument error 
will produce a large absolute error in 
the answer. 

The relative error in the result is 
approximately equal to twice the absolute 
error in the argument divided by 
SIN(2*x). Hence, if x is near a multiple 



Accuracy: 



Arguments 



T 1 

| Relative Error | 
| | *10**6 | 
^ _ T x T -I 

| Range (Distribution | RMS | Maximum | 
i ,_X x x J 

IHETNSR 



|ABS(x) < 
|pi/4 

L___ 


j 


Uniform 


1 

1 

X- 


0.319 


1 


1.92 


1 

1 
1 


|pi/4 < 
|ABS(x) 
|< 1.5 

L 


t 
_X 


Uniform 


T 
1 
1 
1 

L. 


0.465 


1 
I 


1.24 


1 

1 

1 

1 
1 


r 

jpi/4 < 
|ABS(x) 
|< pi/2 

L 


t 
— X— 


Uniform 


1 
1 
1 
1 

X- 


3.14 


1 
1 


170* 


1 

1 

1 

1 
1 


|pi/2 < 
|ABS(x) 
|< 10 

1. 


t 
1 


Uniform 


— T 

1 

1 

1 
L. 


1.25 


1 
I 


70.6* 


1 

1 

1 

1 
1 


|10 < 

|ABS(x) 

|< 100 

L . 


t 


Uniform 


1 

1 

1 

1 
— 1. 


3.57 


1 
.X. 


205* 


1 

1 
1 
1 



♦These maximum errors are those encoun- 
tered in a sample of 5000 points; each 
figure depends very much on the particu- 
lar points encountered near the singular- 
ities of the function. 



Error and Exceptional Conditions: 

P : IHETNSR: ABS(x) > 2**18*pi 
IHETNSD: ABS(x) > 2**18*180 

I : IHETNSR: OVERFLOW 
IHETNSD: OVERFLOW 
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Implementation : 

• Module size: 280 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 



|ABS(x) 


1 
i 


30 


1 
i 


40 


1 
i 


50 


1 
i 


65 


1 


75 


1 
i 


IHETNSR 


|< pi/4 
|> pi/4 


I 

1 

-+- 

1 
i 


4429 
4788 


i 

1 

-+- 

1 
i 


1172 
1262 


T 
.1 

-+- 

1 
1 


336 
368 


i 

1 
-+- 

1 
i 


85. 
95. 


T 

8| 
-+- 

11 
. r_ 


51.0 
55.0 


i 
1 

H 

1 
i 


IHETNSD 


|< 45 
|> 45 


i 

1 

-+" 

1 


4464 
4823 


i 

1 

-+- 

1 
l. 


1184 
1274 


T 

1 

- + " 

1 


341 
373 


i 

1 

-+- 

1 

,„.L. 


87. 
96. 


T 

0| 

-+- 

31 


52.1 
56.0 


i 

1 

H 

1 
.j 



Compute TAN(pi*s/4) and COT(pi*s/4) as 
the ratio of two polynomials: 

TAN(pi*s/4) = s*P(s**2)/Q(s**2)) 
COT(pi*s/4) = Q(s**2)/(s*P(s**2)) 

where P(s**2) is of degree 3 and Q(s**2) 
is of degree 4 in s**2, and maximum 
relative error is 3.4*10**-19. 

Finally, if x < 0, TAN(x) = -TAN(ABS(x) ) . 



Effect of an Argument Error: 

The absolute error in the result is 
approximately equal to the absolute error 
in the argument multiplied by 
(l+TAN(x)**2) . Hence, if x is near an 
odd multiple of pi/2, an argument error 
will produce a large absolute error in 
the result. 

The relative error in the result is 
approximately equal to twice the absolute 
error in the argument divided by 
SIN(2*x). Hence, if x is near a multiple 
of pi/2, an argument error will produce a 
large relative error in the result. 



TAN, TAND (long floating-point real) 



Module Name: IHETNL 



Entry Points: 

Mathematical PL/I Entry 
function name point 

Tan(x radians) TAN(x) IHETNLR 
Tan(x degrees) TAND(x) IHETNLD 

Function: To calculate tan x. 

Method: 

Evaluate 

p = (4/pi) *ABS(x) if x is in radians 
or p = (l/45)*ABS(x) if x is in degrees. 

Let q and r be respectively the integral 
and fractional parts of p. 

If q is even, put s = r; 
If q is odd, put s = 1 - r. 

Let q ± = MOB(q,4). Then 

If q x = 0, TAN(ABSCx)) = TAN(pi*S/4) 

If q ± = 1, TAN(ABS(x)) = COT(pi*S/4) 

If q ± = 2, TAN(ABS(x)) = -COT(pi*S/4) 

If Qi = 3, TAN(ABS(x)) = -TAN(pi*s/4) 



Accuracy: 



IHETNLR 



Arguments 



T 1 

Relative Error 
*10**15 



Range j Distribution 



ABS(x) < 
pi/4 



j Uniform 



pi/4 < 
ABS(x) 
< 1.5 



Uniform 



pi/4 < 
ABS(x) 
< pi/2 



Uniform 



pi/2 < 
ABS(x) 
< 10 



j Uniform 



10 < 
ABS(x) 
£ 100 



| Uniform 

I 

-X 



RMS j Maximum 



0.091 



0.437 



7.75 



18.3 



271 



0.530 



2.31 



416* 



1140* 



13400* 



♦These maximum errors are those encoun- 
tered in a sample of 5000 points; each 
figure depends very much on the particu- 
lar points encountered near the singular- 
ities of the function. 
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Error and Exceptional Conditions: 

P : IHETNLR: ABS(x) > 2**50*pi 
IHETNLD: ABS(x) > 2**50*180 

I : IHETNLR: OVERFLOW 
IHETNLD: OVERFLOW 

Implementation : 

• Module size: 352 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/36 models given 
below are obtained from the appropriate 
entry point in the table: 

r T T T T T n 

|ABS(x) | 30 | 40 | 50 | 65 | 75 | 

L J. X J. J. X J 

IHETNLR 

|< pi/4 | 15440| 3622 | 687 | 154 | 87.2 | 
j. x x x x x .| 

|> pi/4 | 16130 J 3817 j 747 j 169 j 93.9 j 

L . X X X X X J 

IHETNLD 

|< 45 | 15475| 3634 | 691 | 155 | 88.2 | 
|. 1 x x 1 x -| 

j> 45 j 16165J 3829 j 751 j 170 j 94.9 j 

L X X X X X J 



ATAN(X),, ATAND(X), ATAN (Y,X), ATANP (Y,X) 
(short floating-point real) 



Module Name: 



IHEATS 



Method: 

1. ATAN(y,x) 

If x = or ABS(y/x) > 2**24, the answer 
SIGN(y)*pi/2 is returned except for the 
error case x = y = 0. Otherwise 

ATAN(y f x) = ATAN(y/x) if x > 
or ATAN(y f x) = ATAN(y/x) + SIGN(y)*pi 
if x < 0. 

Hence the computation is now reduced to 
the single argument case. 

2. ATAN(x) 

The general case may be reduced to the 
range < x < 1 since 

ATAN(-x) = -ATAN(x), and 
ATAN(1/ABS(X)) = pi/2 - ATAN(ABS (x) ) . 

A further reduction to the range ABS(x) < 
TAN (pi/12) is made by using 

ATAN(X) = pi/6 + ATAN((SQRT(3)*X - 1)/ 
(x + SQRT(3))) . 

Care is taken to avoid the loss of 
significant digits in computing 

SQRT(3)*x - 1. 

For the basic range ABS(x) < TAN(pi/12), 
use an approximation formula of the form 

ATAN(x)/x = a + b*x**2 + c/(d + x**2) 

with relative error less than 2**-27.1. 

3. ATAND(x) and ATAND(y,x) 



Entry Points: 



Mathematical 
function 



PL/I 
name 



Entry 
point 



The treatment is as above with the addi- 
tion of a final conversion of the result 
to degrees. 



Arctan x (radians) ATAN(x) IHEATS1 

Arctan(y/x) (radians) ATAN(y,x) IHEATS2 

Arctan x (degrees) ATAND(x) IHEATS3 

Arctan (y/x) (degrees) ATAND(y,x) IHEATS4 

Function: 

To calculate arctan x or arctan (y/x). 
The result range is: 

Arctan x (radians) ± pi/2 

Arctan (y/x) (radians) ± pi 

Arctan x (degrees) + 90° 

Arctan(y/x) (degrees)! 180° 



Effect of an Argument Error: 

Let t = x or y/x; then the absolute error 
of the answer approximates to the abso- 
lute error in t divided by (1 + t**2) . 
Hence, for small values of t f the two 
errors are approximately the same; howev- 
er,, as t becomes larger the effect of the 
argument error on the answer error dimin- 
ishes. 
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Accuracy: 

T T 1 

| Arguments | Relative Error | 
| | *10**6 | 

j. T + T ^ 

| Range | Distribution | RMS | Maximum | 
x x x x J 

IHEATS1 

r t t t — 1 

I Full range j Tangents of | 0.443 | 0.958 | 
| | numbers | j 

j j uniformly j j j 

| j distributed | | j 

j j between j | j 

| |-pi/2 and | | | 

I lFi/2 | | | 

L X X X J 

IHEATS2 

r t t t 1 

JFull range | y = sin and j 0.449 j 1.42 | 

j | x = cos of | | | 

| j numbers j j j 

| j uniformly j | j 

I j distributed | | j 

j | between j | j 

| | -pi/2 and | | | 

I lFi/2 | | | 

l X X X J 



Error and Exceptional Conditions: 
P : IHEATS2, IHEATS4: x = y = 

Implementation : 

• Module size: 408 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the tables: 



IHEATS2 



l< k 
Ik < 


1 

-+- 

1 


4193 


1 

-+- 

1 


1138 


1 

■+- 

1 


363 


1 VO 1 
1 O 1 
1 H 1 

1 1 

1 1 


1 

-+- 

1 


67.3 | 










|ABS(x) 


1 


5376 


1 


1416 


1 


431 


|126 


1 


77.0 | 


l< 1 


1 




1 




1 




1 


1 




1- 

|1 * 


-+- 
1 




-+- 
1 




-+- 
1 




_X 

1 


-+- 
1 




|ABS(x) 


1 


6082 


1 


1581 


1 


465 


|136 


1 


83.1 | 


|< 1/k 


1 




1 




1 




1 


1 




F 

|£ 1/k 


-+- 
1 


4899 


-+- 
1 


1303 


-X- 

1 


496 


_X 

|117 


-+- 
1 


73.4 | 


L___ 


i 




i 




i 




± 


i 




IHEATS3 


|£ k 

L ___. 


i 

1 
i 


3521 


i 

1 

I 


948 


i 

1 

1 


305 


| 83. 
_x 


T 

9| 
-X- 


51.3 | 


|k < 


1 
1 




1 
1 




1 
1 




T 

1 


T 
1 




|ABS(x) 


1 


4704 


1 


1226 


1 


374 


|103 


1 


61.0 | 


l< 1 

L 


1 

i 




1 
1 




1 
1 




1 
_x 


1 
I 




r 

|1 £ 


l 
1 




1 
1 




1 
1 




T 

1 


— T 
1 




|ABS(X) 


1 


5410 


1 


1391 


1 


408 


|114 


1 


67.1 | 


|< 1/k 
i 


1 

i 




1 
i 




1 
I 




1 
_i 


1 
| 




|S 1/k 


1 
1 


4227 


1 
1 


1114 


1 

1 


339 


| 94. 


— T 

2| 


57.4 | 


L_ — __ __ _ 


i 




i 




i 




j. 


i 




IHEATS4 


1* * 

L_ 


1 
1 


4552 


_ T . 

1 
1 


1228 


i 

1 
I 


389 


|112 
_j. 


1 
1 


70.2 | 


|k < 


T 
I 




1 
1 




1 
1 




T 

1 


T 
1 




|ABS(X) 


1 


5735 


1 


1506 


1 


458 


|131 


1 


69.9 | 


l< 1 

F 

|1 <, 


1 
1 




1 

-+- 

1 




1 

-+- 

1 




1 

-+- 

1 


1 

-+■ 

1 












|ABS(x) 


1 


6441 


1 


1671 


1 


492 


1 142 


1 


86.0 | 


|< 1/k 
i 


1 
1 




1 
i 




1 
1 




1 
_j. 


1 
I 




|£ 1/k 


T 
I 


5358 


1 
1 


1393 


1 
1 


323 


|122 


— T 
1 


76.3 | 























k = TAND(15) 



|ABS(X) 


1 


30 


1 
i 


40 


1 
i _ 


50 


| 65 


1 
. .1 


75 | 
j 


IHEATS1 


1* * 

i ____. 


1 


3162 


1 

I 


858 


1 

1 


279 


| 78. 
_j._ 


j._ 


48.4 | 
j 


r 

|k < 
|ABS(x) 


1 
1 


4345 


1 
1 


1136 


1 
1 


347 


1 

| 97. 


1 

8| 


58.0 | 


l< 1 


1 




1 




1 




1 


1 




J. 

|1 < 


-+- 
1 




-+- 
1 




-+- 
1 




-+- 

1 


1 


H 


|ABS(X) 


1 


5051 


1 


1301 


1 


381 


|108 


1 


64.1 | 


|< 1/k 

l ___. 


1 
i 




1 
i 




1 
1 




1 


1 
X- 


j 


|s 1/k 


T 
1 


3868 


1 
1 


1023 


1 


313 


| 88. 


7| 


54.5 | 























ATAN(X). ATAND(X), ATAN (Y,X), ATAND (Y,X) 
(long floating-point real) 



Module Name: IHEATL 




Entry Points: 




Mathematical 


PL/I 


function 


name 



Entry 
point 



Arctan x (radians) ATAN(x) IHEATL1 

Arctan(y/x) (radians) ATAN(y f x) IHEATL2 

Arctan x (degrees) ATAND (x) IHEATL3 

Arctan (y/x) (degrees) ATAND (Y f X) IHEATL4 
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Function: 



Effect of an Argument Error: 



To calculate arctan x or arctan (y/x) 
The result range is: 

Arctan x (radians) ± pi/2 
Arctan (y/x) (radians)! pi 
Arctan x (degrees) ± 90° 
Arctan (y/x) (degrees) ±180° 



Method : 

1. ATAN(y,x) 

If x = or AES(y/x) £ 2**56, the answer 
SIGN(y)*pi/2 is returned except for the 
error case x = y - 0. Otherwise 

ATAN(y,x) = ATAN(y/x) if x > 
or ATAN(y,x) = ATAN(y/x) + SIGN(y)*pi 
if x < 0. 

Hence the computation is now reduced to 
the single argument case. 

2. ATAN(x) 

The general case may be reduced to the 
range £ x £ 1 since 

ATAN(-x) = - ATAN(x), and 
ATAN(l/ABS(x)) = pi/ 2 - ATAN(ABS (x) ) . 

A further reduction to the range ABS(x) 
< TAN (pi/12) is made by using 

ATAN(x) = pi/6 + ATAN((SQRT(3)*x - 1)/ 
(x + SORT (3))) 

Care is taken to avoid the loss of 
significant digits in computing 

SQRT(3)*x - 1 

For the basic range ABS(x) < TAN dpi/12), 
use a continued fraction of the form 

ATAN(x)/x = 1 + a a .*x*x/(b i ♦ x*x + a 2 / 
(b 2 + x*x + a 3 /(b 3 + x*x + 
a^/tb* + x*x)))) 

with relative error less than 2** (-57. 9). 

3. ATAND(x) and ATAND(y„x) 

The treatment is as above with the addi- 
tion of a final conversion of the result 
to degrees. 



Let t = x or y/x; then the absolute error 
of the answer approximates to the abso- 
lute error in t divided by (1 + t**2) . 
Hence, for small values of t, the two 
errors are approximately the same; howev- 
er„ as t becomes larger the effect of the 
argument error on the answer error dimin- 
ishes. 



Accuracy: 

IHEATL1 

r t 1 

| Arguments | Relative Error | 
| | *10**15 | 
|. T + T .] 

| Range (Distribution | RMS | Maximum | 
j. x x x .] 

j-1 < x < 1| Uniform | 0.0438 | 0.207 | 
l x x x J 



Error and Exceptional Conditions: 
P : IHEATL2, IHEATL4: x = y = 
Implementation : 

• Module size: 544 bytes 

* Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 



| Entry | 
| Point | 

j. X. 

j IHEATL1 | 

L__ -— i_ 


30 | 

1 
X- 

20472J 
i. 


40 
4389 


-+- 

I 


50 
826 


-+- 
1 


65 
181 


1 

1 

-+- 

1 
I 


75 | 
100 | 


| IHEATL2 | 

L X. 


23523| 
i. 


5042 


1 
I 


967 


1 
i 


217 


1 

1 

1 


123 | 


| IHEATL3 | 
| IHEATL4 | 

L X- 


21575| 
24602| 


4667 
5305 


1 
-+- 


871 
1008 


1 
-+- 

.X. 


190 
226 


1 

1 

-+- 

1 


105 | 
128 | 



SINH, COSH (short floating-point real) 

Module Name: IHESHS 

Entry Points: 

Mathematical PL/I Entry 

function name point 

Hyperbolic sin x SINH(x) IHESHSS 

Hyperbolic cos x COSH(x) IHESHSC 
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Function: 

To calculate hyperbolic sin x or hyper- 
bolic cos x. 



Error and Exceptional Conditions: 

H : OVERFLOW in real EXP 
(IHEEXS) . 



routine 



Method : 

For IHESHSS, if ABS(x) < 1, use a poly- 
nomial approximation of the seventh 
degree. 

Otherwise,, 

SINH(x) = EXP(x)/2 - 0.5/EXP(x), 
COSH(x) = EXP(x)/2 + 0.5/EXP(x). 

These two versions of EXP(x)/2 ± 
0.5/EXP(x) are preferable to the equival- 
ent versions of (EXP(x) - l/EXP(x))/2 
because, in floating-point, 0.5 has three 
more significant bits than 1.0. 



Effect of Argument Error: 

The relative error caused in the result 
is approximately as follows: 

SINH: The absolute error in the argument 
divided by TANH(x), i.e., of the 
order of the absolute error in the 
argument for large x, or of the 
relative error in the argument for 
small x. 



Implementation : 

• Module size: 216 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 

r T T T T T 1 

| Entry | 30 j 40 | 50 j 65 | 75 j 
I Point | | | | | | 
h x x 1 x + ., 

j IHESHSS j j j j j j 
|ABS(x)<l|2544 | 692 | 228 | 56. 3| 33.0 | 
|ABS(x)>l|5647 | 1693 | 526 |144 | 91.8 | 
h + x 1 + + ^ 

JIHESHSC j 5500 j 1648 j 509 | 139 j 88.4 j 

L X J. X X X J 



SINH,, COSH (long floating-point real) 



COSH: The absolute error in the argument 
multiplied by TANH(x)„ i.e., of the 
order of the absolute error in the 
argument. 

Thus, for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 

Accuracy: 



IHESHSS 



Arguments 



"T 1 

| Relative Error | 
I *10**6 I 



t 1 t 1 

Range | Distribution! RMS | Maximum | 

X X X J 



o < I II I 

ABS(x) < lj Uniform | 0.200 j 0.932 | 



+ + + H 

1 < I II I 

ABS(x) < 2 j Uniform j 0.221 j 0.950 j 



IHESHSC 



T T T 1 

ABS(x) < lj Uniform j 0.367 j 0.908 | 



1 < I II I 
ABS(x) < 2| Uniform j 0.192 | 0.700 j 
x x x J 



Module Name: 
Entry Points: 



IHESHL 



Mathematical PL/I Entry 

function name point 

Hyperbolic sin x SINH(x) IHESHLS 

Hyperbolic cos x COSH(x) IHESHLC 



Function: 

To calculate hyperbolic sin x or hyper- 
bolic cos x. 



Method •- 

For IHESHSS, if ABS(x) < 0.3465736, com- 
pute SINH(x)/x using a polynomial approx- 
imation of degree 5 in x**2, with rela- 
tive error less than 2**-61.9 



Otherwise, compute s = EXP(ABS(x)); then 

COSH(x) = (s + l/s)/2 

SINH(X) = SIGN(X)*(S - l/s)/2 
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Effect of an Argument Error: 

The relative error caused in the 
is approximately as follows: 



result 



SINH: The absolute error in the argument 
divided by TANH(x), i.e., of the 
order of the absolute error in the 
argument for large x, or of the 
relative error in the argument for 
small x. 

COSH: The absolute error in the argument 

argument multiplied by TANH(x), 

i.e., of the order of the absolute 
error in the argument. 

Thus, for large values of x, even the 
round-off error of the argument causes a 
substantial relative error in the answer. 



Accuracy: 

r ■ t 

| Arguments | Relative Error 

j | *10**15 

j. T „ + T ._ 

| Range | Distribution! RMS | Maximum 
l_ j. j. x 

IHESHLS 

|ABS(x) < | Uniform | 0.0530| 0.217 

|0.34657 | j | 

|. + + + 

|0.34567 < | || 

|ABS(x) < 5| Uniform j 0.0870J 0.359 

L J. J. X 

IHESHLC 

r T ■ T T — 

|ABS(x) < 5| Uniform j 0.123 j 0.429 

L X , X X , 



Error and Exceptional Conditions: 

H : OVERFLOW in real EXP routine 
(IHEEX1). 



Implementation : 

• Module size: 264 bytes 

• Execution times : 

Approximate execution times in micro- 
seconds for System/360 models given 
below are obtained from the appropriate 
entry point in the table: 

|ABS(x) | 30 | 40 | 50 | 65 | 75 | 

L X X X X X J 

IHESHLS 

|< 0.347| 9024| 2279 | 450 | 101 | 59.0 1 

|. X X X X X ^ 

|0.347 <| || || | 

|ABS(x) | 18634| 4338 | 938 | 215 | 125 | 
I* 174.6| I I I I I 

L X X X X X J 

IHESHLC 

|< 174.6| 18493| 4300 | 924 | 211 | 123 | 
L X X X X X J 



TANH (short floating-point real) 

Module Name: IHETHS 

Entry Point: IHETHS0 

Function: To calculate hyperbolic tan x. 

Method : 

1. ABS(x) < 2**-12 
Return x as result. 

2. 2**-12 < ABS(x) < 0.54931 

Use a transformed continued fraction of 
the form: 

TANH(x)/x = l-((x**2 + a)/(x**2 + b + 
c/x**2)) 

with relative error less than 2**-27. 

3. 0.54931 < x < 9.011 

Use TANH(x) = 1 - 2/(EXP(2*X) + 1). 

4. x > 9.011 
Return result 1. 

5. x < -0.54931 
TANH(x) = -TANH(-x). 
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Effect of an Argument Error: 

The relative error caused in the result 
is approximately twice the absolute error 
in the argument divided by SINH(2*x). 
Thus for small values of x it is of the 
order of the relative error in the argu- 
ment, and as x increases the effect of 
the argument error is diminished. 

Accuracy: 



Arguments 



| Relative Error 
I *10**6 



t T t 

Range J Distribution! RMS | Maximum 



-0.5 < x | Uniform | 0.174 | 0.867 
< 0.5 | | I 



|-9 < x < 9 | Uniform | 0.0720 | 
l x x x. 



0.782 



Implementation : 

• Module size: 200 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 



r -t 

1 1 30 

L J. 


■T 

| 40 
_i 


"T 

| 50 
_x 


| 65 

_x 


-T 1 

1 75 | 

_x J 


|ABS(x) < | 791 


| 263 


|102 


| 28.7 


| 21.7 | 


|2**-12 | 

L_ i 


1 
_x 


1 
_x 


1 
_x 


1 1 

_i J 


|2**-12 < |3033 


| 785 


|231 


| 64.1 


| 43.9 | 


|ABS(x)<0.5| 

L_ _x 


1 
.j. 


1 
_x 


1 
_x 


1 1 
_x J 


| 0.5 < x |5934 


|1805 


| 562 


1 152 


|117 | 


l< 9 | 

|x £ 9 J1095 


1 

-X 

| 363 


1 

-X 

|139 


1 

| 40.5 


1 1 
| 35.2 | 


L_ X 


.x . 


-X 


_x 


_x J 



TANH (long floating-point real) 

Module Name: IHETHL 

Entry Point: IHETHL0 

Function: To calculate hyperbolic tan x. 

Method: 

1. ABS(x) < 0.54931 

Compute TANH(x)/x using a rational 
approximation, with relative error less 
than 2**-64.5 



2. 0.54931 < X < 20.101 

TANH(x) = 1 - 2/(EXP(2*x) + 1). 

3. x £ 20.101 
Return result 1. 

4. x < -0.54931 
TANH(x) = -TANH(-x) 

Effect of an Argument Error: 

The relative error caused in the result 
is approximately twice the absolute error 
in the argument divided by SINH(2*x). 
Thus for small values of x it is of the 
order of the relative error in the argu- 
ment, and as x increases the effect of 
the argument error is diminished. 

Accuracy: 



j Arguments 


— t- 

1 
1 

_ X 


Relative Error j 
♦10**15 | 


r~ T~ ' — T 

| Range | Distribution | 

L X__ -L 


RMS | 

x_ 


Maximum j 


|ABS(x) < | 
|0.54931 | 

L X 


Uniform 


T 
1 
1 

X- 


0.0440| 

1 

X- 


0.211 | 


|0.54931 < | 
|ABS(x) < 5| 

L . _X— 


Uniform 


— T 
1 
1 


, — .j.- 

1 

0.0250| 


0.199 | 



Implementation : 

• Module size: 280 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/ 360 models given 
below are obtained from the table: 



|ABS(X) j 
L i_ 


30 | 
. x_ 


40 


1 
1 


50 


1 65 
_i 


1 
—X- 


75 


1 
1 


|< 0.549| 


12745| 


3030 


1 

1 


564 


|123 


T 
1 


67. 


1 
9| 


L X_ 


X- 




1 




_x 


X- 




1 


|0.549 £| 


T 

1 




1 
1 




T 

1 


t 

1 




1 
1 


|ABS(x) | 


16400| 


3918 


1 


878 


| 205 


1 


119 


1 


|< 20*1 | 


1 




1 




1 


1 




1 


h +- 

|£ 20.1 | 

L X_ 


X. 

1239| 
x_ 


372 


-+- 

1 
_x_ 


135 


_X 

| 39. 
_x 


-+- 

3| 
-X_ 


25. 


H 

5| 



ATANH (short floating-point real) 

Module Name: IHEHTS 
Entry Point: IHEHTS0 
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Function: To calculate hyperbolic arctan x. 
Method: 

1. AES(x) < 0.2 

Use a rational approximation of the form: 
ATANH (x) = x + x ** 3/ (a + b*x**2) 

2. 0.2 < ABS(X) < 1 

ATANH (x) = -SIGN (x)*0.5*LOG( (0.5 - 

ABS(x/2))/(0.5 + ABS(x/2))) 

Effect of an Argument Error: 

The absolute error caused in the result 
is approximately equal to the absolute 
error in the argument divided by 
(1 - x**2). Thus as x approaches 
+l' or -1, relative error increases rapid- 
ly. Near x = 0., the relative error in 
the result is of the order of that in the 
argument. 

Accuracy: 



Arguments 



-t 1 

| Relative Error 
I (*io**6) 



Range j Distribution! RMS | Maximum 



-0.8 < x j Uniform j 0.440 j 1.32 
S 0.8 I I I 



-0.9 < x | Uniform | 0.389 | 1.14 
< 0.9 | | | 

X . X . X 



Error and Exceptional Conditions: 
P : ABS(x) £ 1 



Implementation : 

• Module size: 192 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

r t t t "T t n 

|ABS(x) I 30 I 40 j 50 I 65 I 75 I 

^ X X X— X X ^ 

|< 0.2 | 2520 | 667 | 208 | 52. 0| 31.3 | 

^ X X X _ + X -! 

|0.2 < | | | || | 
|ABS(x) | 7091 | 1829| 606 |163 | 94.8 j 

l< 1 I II I I I 
L X X X -X „X J 



ATANH (long floating-point real) 

Module Name: IHEHTL 
Entry Point: IHEHTL0 

Function: To calculate hyperbolic arctan x. 

Method: 

1. ABS(x) < 0.25 

Use a Chebyshev polynomial of degree 8 in 
x**2 to compute ATANH (x)/x. 

2. 0.25 < ABS(x) < 1 

ATANH (x) = -SIGN (x)*0.5*LOG( (0.5 - 

ABS(x/2) )/(0.5 + ABS(x/2))) 

Effect of an Argument Error: 

The absolute error caused in the result 
is approximately equal to the absolute 
error in the argument divided by 
(1 - x**2) . Thus as x approaches +1 or 
-1,, relative error increases rapidly. 
Near x = 0, the relative error in the 
result is of the order of that in the 
argument. 

Accuracy: 



Arguments 



|. T 

| Range | Distribution 



h 

|ABS(x) <, | Uniform 

|0.25 | 



|ABS(x) < | Uniform 
|0.95 | 

L X 



Relative Error 
♦10**15 



RMS 



Maximum 



0.0650 j 0.223 
I 



0.133 | 0.397 

I 

x 



Error and Exceptional Conditions: 
P : ABS(x) £ 1 



Chapter 3. Mathematical: Functions with Real Arguments 55 



Iirpleirentation : 

• Module size: 272 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the tafcle: 

|ABS(x) | 30 | 40 | 50 | 65 | 75 | 
j. x + 1— + 1 -| 

|< 0.25 | 12252| 3037 | 562 | 121 | 68 | 
J. + + + + 1 1 

10.25 < | | | || | 
|ABS(x) j 20448 j 4900 j 1040 j 242 j 137 j 
l< 1 I I I II I 

L X X X X X J 



3. k < x < 13.306 

ERFC(x)*x*EXP(x**2) is computed by using 
a Chebyshev interpolation polynomial of 
degree 6 in x**-2, with relative error 
less than 1.2 * 2**-23. 

If x < 3.9192, ERF(x) = 1 - ERFC(x) 
If x > 3.9192, ERF(x) = 1 

4. x > 13.306 

Results 1 and are returned for ERF(x) 
and ERFC(x) respectively. 

5. x < 

ERF(x) = - ERF(-x) 
and ERFC(x) = 2 - ERFC(-x). 



ERF, ERFC (short floating-point real) 



Effect of an Argument Error: 

The absolute error caused in the result 
is approximately equal to the absolute 
error in the argument multiplied by 
EXP(-x**2) . 



Module Name: 
Entry Points: 



IHEEFS 



Mathematical 
function 

Error function (x) 
Complement of error 
function (x) 



Function: 



ERF(x) 
ERFC(x) 



Entry 
point 

IHEEFSF 
IHEEFSC 



To calculate the error function of x or 
the complement of this function. 



ERF(x): As the magnitude of the argument 
increases from 1, the effect of an argu- 
ment error diminishes rapidly. For small 
x, the relative error of the result is of 
the order of the relative error of the 
argument. 

ERFC(x): For x > 1, ERFC(x) is approxi- 
mately EXP(-x**2)/(2*x) . Thus the rela- 
tive error in the result is approximately 
equal to the relative error in the argu- 
ment multiplied by 2*x**2. For negative, 
or small positive, values of x, the 
relative error in the result is approxi- 
mately equal to the absolute error in the 
argument multiplied by EXP(-x**2). 



Method: 

1. < x < 1.317 

Compute ERF(x)/x by using a Chebyshev 
interpolation polynomial of degree 6 in 
x**2„ with relative error less than 
2**-24. 

ERFC(x) = 1 - ERF(x) 

(ERFC(x) > 1/16 in this range). 

2. 1.317 < x < k„ where k = 2.04000092 

Compute ERFC(x) by using a Chebyshev 
interpolation polynomial of degree 7 in 
(x-k) , with absolute error less than 
1.3 * 2**-30. 

ERF(x) = 1 - ERFC(x). 

(ERFC(x) > 1/256 in this range). 



Accuracy: 



Arguments 



| Relative Error 
♦10**6 



r t • 1 t 

| Range (Distribution | RMS | Maximum 
l x x x 

IHEEFSF 

r 

|ABS(x) 
|< 1.3 

^ 

1 1. 3 < 

|ABS(x) 

|2 < 
|ABS(X) 
|< 3.9 
L 



Uniform j 0.139 j 0.934 j 

I I I 
x + x -i 

I I I 

Uniform j 0.0372J 0.263 j 

I I I 
x .. x + -i 

I I I 

Uniform j 0.0347 j 0.0605 j 



x x x J 
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IHEEFSC 



ERF, ERFC (long floating-point real) 



1-3.8 < x 




Uniform 




0.297 




0.930 


l< o 

L 


—4. 




X_ 




1 




|0 < x 


T 


Uniform 


— T 


0.505 


1 


3.80 


|< 1.3 

L 


— X_ 




j.. 




1 




|1.3 < x 


T 


Uniform 


— T 


0.314 


1 


4.08 


|< 2 

L 


— 4— 




—4. 




I 




|2 < x 


T 


Uniform 


1 


0.367 


1 


1.45 


|< 3.9 

L _ 


— 4— 




x_ 




1 




| 3.9 < x 


T 


Uniform 


— -j- 


9.09 


1 


17.1 


|< 13.3 





























Implementation : 

• Module size: 376 bytes 

• Execution times: 

Approximate execution tiroes, in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 

r T T — " T 'V T 1 

|ABS(x) I 30 I 40 I 50 I 65 I 75 I 

L 4. X A_ X .4. J 

IHEEFSF 

|< 1.32 | 4354| 1191 | 392 |100 | 57.2 | 

|1.32 < | | | || | 
|ABS(x) | 4613| 1266 j 418 J110 | 62.0 j 
|£ 2.04 | | | | | | 
|. + + + 4. + ^ 

|2.04 < | | | | | 1 
|ABS(x) | 10013| 2843 | 868 | 228 |140 | 
l< 3.92 | | | | | | 
h + + x + „x 1 

|£ 3.92 | 1530| 473 | 183 | 50. 7| 32.0 | 

L X 1 4. X X J 

IHEEFSC 

r T T— T T ~T 1 

\<, 1.317J 4412| 1204 | 385 |103 | 58.8 | 

|. + + 4. _ + - + 1 

|1.32 < | | | || | 
|ABS(x) | 4582| 1262 | 405 J110 | 63.4 | 
|< 2.04 | | | | | | 
j. x + - x _ + x___ ., 

|2.04 < | | | | | | 
|ABS(x> j 9982 j 2839 j 854 j 228 |142 | 
l< 3.92 | | | | | I 

|3.92 < | | | || | 

|ABS(x) | 10168| 2897 | 879 | 236 | 147 | 

l< 13.31| | | III 
|. x x x _x____x 4_ 

|£ 13.31| 1598| 499 | 182 | 54. 5| 35.5 | 

L X X X X _X J 



Module Name: 
Entry Points: 



IHEEFL 



Mathematical 
function 

Error function (x) 
Complement of error 
function (x) 



Function: 



PL/I 
name 

ERF(x) 
ERFC(x) 



Entry 
point 

IHEEFLF 
IHEEFLC 



To calculate the error function of x or 
the complement of this function. 



Method : 

1. < x < 1 

Compute ERF(x)/x by using a Chebyshev 
interpolation polynomial of -degree 11 in 
x**2, with relative error less than 
1.07*2**-57. 

ERFC(x) = l-ERF(x) . 

(ERFC(x) > 1/16 in this range). 

2. 1 < x < 2.04000092 

Compute ERFC(x) by using a Chebyshev 
interpolation polynomial of degree 18 in 
(x - 1.999999), with absolute error less 
than 1.5*2**-61. 

ERF(x) = l-ERFC(x) . 

(ERFC(x) > 1/256 in this range). 

3. 2.04000092 < x < 13.306 

ERFC(x) is computed by using a Chebyshev 
interpolation polynomial of degree 20 in 
x**-2 for ERFC(x)*x*EXP(x**2) , with rela- 
tive error ranging from 2**- 5 3 at 
2.04000092 to 2**-51 at 13.306. 

If x < 6.092, ERF(x) = l-ERFC(x). 
If x 2: 6.092, ERF(x) = 1 

4. x k 13.306 

Results 1 and are returned for ERF(x) 
and ERFC(x) respectively. 

5. x < 

ERF(x) = -ERF(-x) 
and ERFC(x) = 2-ERFC(-x). 
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Effect of an Argument Error: 

The absolute error caused in the result 
is approximately equal to the absolute 
error in the argument multiplied by 
EXP(-x**2). 

ERF(x): As the magnitude of the argument 
increases from l, f the effect of an argu- 
ment error diminishes rapidly. For small 
x, f the relative error of the result is of 
the order of the relative error of the 
argument. 

ERFC(x): For x > 1,, ERFC(x) is approxi- 
mately EXP(-x**2)/(2*x). Thus the rela- 
tive error in the result is approximately 
equal to the relative error in the argu- 
ment multiplied by 2*x**2. For negative, 
or small positive, values of x, the 
relative error in the result is approxi- 
mately equal to the absolute error in the 
argument multiplied by EXP(-x**2). 

Accuracy: 



| Arguments 




Relative 


Error 










♦10 s 


f*15 


i 














r 

| Range 


~T~ T" 

| Distribution! 


RMS 


~T"~ 

| Maximum 


L _ 


_L 




i 




j_ 




IHEEFLF 


r 

|ABS(x) < 


1 


Uniform 


i 


0.0280 


j 


0.202 


|1.317 

L _ 


-X — 




i 




i 




|1.317 < 


1 




1 




1 




|ABS(x) 




Uniform 




0.0107 




0.0291 


|< 2.04 
L 


-X— 




I 




i 




| 2. 04 < 


1 




1 




1 




|ABS(x) 




Uniform 




0.008031 


0.0170 


|< 6.092 














L _ 


J. 




i 




i 




IHEEFLC 


r ~ 

j-6 < x 


i 


Uniform 


K 


0.0684 


i 


0.188 


l< ° 














f 

|0 £ x 


-+- 


Uniform 


-+- 


0.0762 


-+- 


0.352 


|< 1.317 

L_ 


1 




I 




l 




|1.317 < x 


1 


Uniform 


1 


0.127 


1 


0.445 


|£ 2.04 














i- 

j 2.04 < x 


-+- 


Uniform 


-+- 


1.24 


-+- 


4.02 


|< 4 

L ___ _ __ 


-X— 




i 




i 




r 

|4 < x 


1 


Uniform 


1 


1.40 


1 


5.02 


|< 13.3 





























Implementation : 

• Module size: 744 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/ 360 models given 
below are obtained from the appropriate 
entry point in the table: 

|ABS(x) | 30 | 40 | 50 | 65 | 75 | 

l j. x x x ; X J 

IHEEFLF 

|£ 1.00 | 16567| 4154| 805 |180 | 103 | 
j. x X x + x ^ 

|1.00 < | I - I I I I 
|ABS(x) | 24413| 6095| 1175 | 263 | 147 | 
I* 2.04J I I I I I 

|2.04 < | | || | | 
|ABS(x) | 45574 j 1105 j 2142 | 477 | 269 | 
l< 6.09 | | | | | | 

|£ 6.09 | 2493| 707| 222 | 58. 0| 36. 3| 

L X X X X X J 

IHEEFLC 

|< 1.00 | 16905| 4240| 837 |188 | 106 | 
j. x x x x x -| 

|1.00 < | | | || | 

|ABS(x) j 24263| 6065 j 1166 | 261 j 147 | 
|£ 2.04 || | || | 

|._- x x x x x ■! 

12.04 < | | | || | 

|ABS(x) | 45424| 10974J 2133 | 474 | 270 | 

l< 6.09 | | | | | I 

j. 1 x x x x .| 

16.09 < | | | || | 
|ABS(x) | 45624| 11040 | 2160 | 482 | 274 | 
l< 13.3 | | | | | | 
H x x x x x ^ 

|£ 13.3 | 2408| 693| 219 | 57. 4| 37. 5| 

L X X X X X J 



FUNCTIONS WITH COMPLEX ARGUMENTS 

■j SQRT (short floating-point complex) 

Module Name: IHESQW 
Entry Point: IHESQW0 



Function: 

To calculate the principal value of the 
square root of z, i.e., -pi/2 < argument 
of result < pi/2. 
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Method: 

Let z = x + yl, and 
SQRT(z) = u + vl. 

1. x = y = 
Then u = v = 0. 

2. x £ 

Then u = SQRT < (ABS (x) + AES(x + yI))/2) 
and v = y/(2*u) . 

3. x < 

Then u = y/(2*v) 
and v = S(y)*SQRT((ABS(x) + 
AES(x + yI))/2) 

where S(y) = 1 if y £ 
= -1 if y < 



Effect of an Argument Error: 

Let z = r*EXP(hI), and 
SQRT(z) = s*EXP(kI). 

Then the relative error in s is approxi- 
mately half the relative error in r, and 
the relative error in k is approximately 
equal to the relative error in h. 



Accuracy: 

r t t 

| Arguments | Relative Error | 
| | *10**6 | 
^ T + T _„ ^ 

| Range (Distribution | RMS | Maximum | 
j. + _ x x .j 

| Full range | Exponential | 0.513 | 1.51 | 
| | radially, | | | 
I | uniform j j j 
j j round origin | j j 
l x , x , x J 



Error and Exceptional Conditions: 
I : OVERFLOW 



H 



OVERFLOW 
(IHEAEW) 



in complex ABS routine 



Implementation : 

• Module size: 152 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

IModule | , 30 | 40 | 50 | 65 | 75 | 

I Name I I I I I I 
j. + x x x x ^ 

IIHESQW | 11130| 3023 | 1006 | 265 | 164 | 
L X X X X X J 



SQRT (long floating-point complex) 

Module Name: IHESQZ 
Entry Point: IHESQZ 

Function: 

To calculate the principal value of the 
square root of z, i.e. , -pi/2 < argument 
of result £ pi/2. 

Method: 

Let z = x + yl, and 
SQRT(z) = u + vl. 

1. x = y = 
Then u = v = 0. 

2. x £ 

Then u = SQRT ( (ABS (X) + ABS(x + yI))/2) 
and v = y/(2*u) . 

3. x < 

Then u = y/(2*v) 
and v = S(y)*SQRT((ABS(x) + 
ABS(x + yI))/2) 

where S(y) = 1 if y £ 
= -1 if y < 

Effect of an Argument Error: 

Let z = r*EXP(hI), and 
SQRT(z) = s*EXP(kI). 

Then the relative error in s is approxi- 
mately half the relative error in r, and 
the relative error in k is approximately 
equal to the relative error in h. 
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Accuracy: 

r t 1 

| Arguments | Relative Error | 
j | *10**15 | 
j. T x T .j 

| Range | Distribution! RMS | Maximum | 
^ + + x ^ 

j Full range j Exponential j 0.263 j 1.54 j 
| | radially, | | | 

j j uniform | j j 

| j round origin | j j 

l x j. x J 



Error and Exceptional Conditions: 

I : OVERFLOW 

H : OVERFLOW in complex ABS routine 
(IHEAEZ) 

Implementation 

• Module size: 144 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for System/360 models given 
below are obtained from the table: 

r T T T T T" 1 

| Module | 30 | 40 j 50 | 65 | 75 | 
I Name I I 1 I I I 

|. + x + + + ^ 

IIHESQZ | 26996| 5957 | 1352 | 341 | 203 | 

L X X X X X J 



Accuracy: 



Arguments 



T 1 

I Relative Error 
| *10**6 



r t T t 

| Range (Distribution | RMS | Maximum 



1- 

|ABS(x) 


__x_ 
1 


Uniform 


__x 

| 1.32 


-+- 


3.14 


|< 170 


1 




1 






|ABS(y) 


1 




1 






|< pi/2 

L 


1 

_ _j 




1 

L 


x 




r 

|ABS(x) 


— -j-- 
1 


Uniform 


| 1.31 


t — 


3.34 


|< 170 


1 




1 






|pi/2 < 
|ABS(y) 
|< 20 


1 
1 
1 




1 
1 
1 



















Error and Exceptional Conditions: 

O : ABS(y) £ 2**18*pi : error caused in 
real SIN routine (IHESNS) 

H : OVERFLOW in real EXP routine (IHEEXS) 



Implementation : 

• Module size: 136 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 



|Module | 30 | 40 | 50 | 65 | 75 | 
I Name | j I I I I 

|IHEEXW | 14240| 4058 | 1244 | 326 | 200 | 
L X X X X X J 



EXP (short floating-point complex) 



Module Name: IHEEXW 

Entry Point: IHEEXWO 

Function: To calculate e to the power z. 

Method: 

Let z = x + yl. 

Then REAL (EXP (z)) = EXP (x) *COS (y) 
and IMAG(EXP(z)) = EXP(x)*SIN(y) . 

Effect of an Argument Error: 

Let EXP(x + yl) = s*EXP(kI). 

Then k = y # and the relative error in s 
is approximately equal to the absolute 
error in x. 



EXP (long floating-point complex) 



Module Name: 
Entry Point: 



IHEEXZ 
IHEEXZ0 



Function: To calculate e to the power z. 

Method: 

Let z = x + yl. 

Then REAL(EXP(z)) = EXP(x)*COS(y) 
and IMAG(EXP(z)) = EXP(x)*SIN(y) . 
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Effect of an Argument Error: 

Let EXP(x + yl) = s*EXP(kI). 

Then k = y, and the relative error in s 
is approximately equal to the absolute 
error in x. 

Accuracy: 



Arguments 



T , n 

Relative Error 
♦10**15 



| Range 



Eistribution 



jABS(x) < 1 
|ABS(y) 
|< pi/2 



ABS(x) 
j< 20 
|ABS(y) 
|< 20 



Uniform 



Uniform 



RMJ 



Maximum 



0.136 j 0.478 



1.28 | 2.29 



Method: 

Let LOG(x + yl) = u + vl. 

Then u = LOG(ABS(x + yl)) 

= LOG(SQRT(x**2 + y**2)) 
= LOG(x**2 + y**2)/2 

and v = ATAN(y,x). 

In computing u„ the exponents of x and y 
are modified if necessary to avoid OVER- 
FLOW or UNDERFLOW, with the appropriate 
correction being applied after the logar- 
ithm has been taken. 



Effect of an Argument Error: 

Let z = r*EXP(hI) and LOG(z) = u + vl. 

Then the absolute error in u is approxi- 
mately equal to the relative error in r. 
For the absolute error in 
v(= h = ATAN(y,x)), see corresponding 
paragraph for module I HEATS. 



Error and Exceptional Conditions: 

O : ABS(y) > 2**50*pi : error caused in 
real SIN routine (IHESNL) 



Accuracy: 



Arguments 



"T T 

| Relative Error | 
*10**6 



H : OVERFLOW in real EXP routine (IHEEXL) | Range | Distribution 



RMS 



Maximum 



Implementation: 

• Module size: 136 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for System/360 models given 
below are obtained from the table: 

[Module j 30 | 40 | 50 | 65 | 75 | 
I Name | I I I I I 
j. + + + + + -I 

JIHEEXZ j 42838 j 10560 j 2174 | 505 j 287 j 

L J. .-J. X JL i. J 



|Full range | Exponential | 0.150 



| except 
| within 
|10**-6 of 
|1+0I 

L 



| radially, 

uniform | 
round origin j 



2.27 | 



Error and Exceptional Conditions : 

error in real LOG routine 



: x = y = 0, 
(IHELNS) 



LOG (short floating-point complex) 

Module Name: IHELNW 

Entry Point: IHELNW0 

Function: 

To calculate the principal value of the 
natural log of z,, i.e.,, -pi < imaginary 
part of result < pi. 
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Implementation : 

• Module size: 272 bytes 

» Execution tines: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the table: 

r T T T T T 1 

I I 30 | 40 | 50 | 65 | 75 | 
j. + + + + x ^ 

| (i) | 11511| 3414 j 1078 | 308 | 183 j 
j. + + + + x H 

| (ii) | 11688| 3489 | 1104 | 318 | 190 | 
|. + + + + + ^ 

| (iii) | 11800 1 3520 | 1117 J 321 | 193 | 

L _X X X X X J 

(i) ABS(x) and ABS(y) < SQRT(8)*16**31 
and either ABS(x) or ABS(y) Z 16**-30 

(ii) Either ABS(x) or ABS(y) 
2 SQRT(8)*16**31 

(iii) ABS(x) and ABS(y) < 16**-30 



LOG (long floating-point complex) 

Module Name: IHELNZ 

Entry Point: IHELNZ0 

Function: 

To calculate the principal value of natu- 
ral log of z, i.e., -pi < imaginary part 
of result 3 pi. 

Method: 

Let LOG(x + yl) = u + vl. 

Then u = LCG(ABS(x + yl)) 

= LOG(SQRT(x**2 + y**2)) 
= LOG(x**2 + y**2)/2 

and v = ATAN(y # x). 

In computing u, the exponents of x and y 
are modified if necessary to avoid OVER- 
FLOW or UNDERFLOW, with the appropriate 
correction being applied after the logar- 
ithm has been taken. 

Effect of an Argument Error: 

Let z = r*EXP(hI) and LOG(z) = u + vl. 

Then the absolute error in u is approxi- 
mately equal to the relative error in r. 
For the absolute error in 
v(= h = ATAN(y,x)) see the corresponding 
paragraph for module IHEATL. 



Accuracy: 



Arguments 



I" 



Range 



T n 

| Relative Error | 
j *10**15 | 
. T + T .| 

| Distribution | R.M.S. | Maximum | 



|. x + x 

j Full range | Exponential j 0.05581 1.46 
j except j radially, j 
| within | uniform j 
|10**-6 of j round origin j 
|1+0I | | 

L X X 



Error and Exceptional Conditions: 

: x = y = 0, error caused in log rou- 
tine (IHELNL) 



Implementation : 

• Module size: 288 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for System/360 models given 
below are obtained from the table: 

r t t t t t 1 

| | 30 | 40 | 50 | 65 | 75 | 
j. x x 1 x x -i 

| (i) | 44101| 10166| 2086 j 480 | 274 | 

| (ii) | 44398| 10316| 2155 | 507 | 290 | 
j. x x x x x 1 

| (iii) | 44438J 10325J 2156 | 507 j 290 | 

L X X X X X J 

(i) ABS(x) and ABS(y) < SQRT(8)*16**31 
and either ABS(x) or ABS(y) 2: 16**-26 

(ii) Either ABS(x) or ABS(y) 
£ SQRT(8)*16**31 

(iii) ABS(x) and ABS(y) < 16**-26 



SIN, SINH, COS, COSH (short floating-point 
complex) 



Module Name: IHESNW 
Entry Points: 



Mathematical 
function 

Sin z 

Hyperbolic sin z 
Cos z 
Hyperbolic cos z 



PL/I 
name 

SIN(z) 
SINH(z) 
COS(z) 
COSH(z) 



Entry 
point 

IHESNWS 
IHESNWZ 
IHESNWC 
IHESNWK 
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Function: 

To calculate sin z or hyperbolic sin z, 
or cos z or hyperbolic cos z. 



Method: 

Let z = x + yi. 



Then REAL(SIN(z)> 
and IMAG(SIN(z)) 

REAL(COS(z)) 
and IMAG(COS(z)) 



= SIN(x)*COSH(y) 
= COS(x)*SINH(y); 

= COS(x)*COSHCy) 
= -SIN(x)*SINH(y); 



REAL (SINH (z)) = COS(y)*SINHCx) 
and IMAG(SINH(z)) = SIN(y)*COSH(x) ; 

REAL(COSH(z)) = COS(y)*COSH(x) 
and IMAG(COSH(z)) = SIN(y)*SINH(x) . 

To avoid making calls to evaluate SINH 
and COSH separately, and thus frequently 
having to evaluate EXP twice for the same 
argument, SINH(u) is computed as follows: 

1. u > 0.3465736 

SINH(u) = (EXP(u) - l/EXP(u))/2. 

2. £ u £ 0.3465736 

SINH(u)/u is approximated by a polynomial 
of the form a + a ± *u**2 + "a a *u**4 (which 
has a relative error of less than 
2**-26.4). 

3. u £ 

SINH(u) = -SINH(-u). Then 

COSH(u) = SINH(ABS(u)) + 1/EXPCABS (u) ) . 



Effect of an Argument Error: 

Combine the effects on SIN, COS, SINH and 
COSH according to the method of evalua- 
tion described in the above paragraph. 



Accuracy: 

r t t 

| Arguments | Relative Error | 
| | *10**6 | 
j. T + _— T + 

| Range | Distribution! RMS | Maximum | 

l x x __x J 

IHESNWS 

r T T T T 

|ABS(x) j Uniform j 0.721 j 2.07 j 

|£ 10, | || I 

|ABS(y) <S 1| || | 

L X X X J 

IHESNWZ 

|ABS(x) | Uniform | 0.561 | 1.86 | 

I* 10, | | | | 

|ABS(y) £ 1| || | 

L X X X J 

IHESNWC 

r t t ■ — T 1 

|ABS(x) I Uniform | 0.546 | 2.00 j 

I* 10, | || | 

|ABS<y> | III 

I* 20 | || | 

L X * X X J 

IHESNWK 

|ABS(x) | Uniform | 0.558 | 2.35 | 

I* 10, | III 

|ABS<y> | III 

1^ 20 | || | 

l X X X J 



Error and Exceptional Conditions: 

O : IHESNWS, IHESNWC: 

ABS(x) £ 2**18*pi : error caused in 
real SIN routine (IHESNS) 

IHESNWZ, IHESNWK: 

ABS(y) £ 2**18*pi : error caused in 
real SIN routine (IHESNS) 

H : OVERFLOW in real EXP routine (IHEEXS) 
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Implementation : 

• Module size: 320 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the tables: 

SIN, COS : AES(y) > 0.3465736 
SINH,COSH: AES(x> > 0.3465736 



r — t t 

| Entry j 30 j 40 
| Point | j 

|IHESNWS| 15826| 1508 

1 | L « 


1 

1 

-+- 

1 

I 


50 
648 


1 65 
1 

| 363 
_x_ 


1 

1 

■+- 

1 
l 


75 | 
223 | 


r 1 t 

JIHESNWCJ 15898 | 1518 
f. + + 

JIHESNWZJ 15930 j 1520 
|. + X 

JIHESNWKJ 15900 | 1519 

I L X _ — 


r 

1 
-+- 

1 

-+- 

1 
i 


653 
653 
655 


| 366 
-X 

| 366 
-X 

| 367 
j. 


1 
1 

1 

-+- 

1 
i 


225 | 

226 | 

227 | 


SIN,, COS : ABS(y) < 0.3465736 
SINH, COSH: ABS (x) < o. 3465736 


r t t ~~ 

JIHESNWSJ 16896 j 1585 
1 X I 


.1 

1 
I 


674 


| 381 
i 


i 

1 
I 


232 | 


r t 1 

JIHESNWCJ 16968| 1595 

L i X 


1 

1 

I 


679 


| 384 
_x 


1 

1 
1 


234 | 


|IHESNWZ| 17000| 1596 
|. X X 

JIHESNWKJ 16970 j 1595 


1 

1 

-+- 

1 


679 
681 


| 384 

-X 

| 384 


1 

1 

-+- 

1 


234 | 

235 | 



Method: 

Let z = x + yl. 

Then REAL(SIN(z)> = SIN (x) *COSH(y) 
and IMAG(SIN(z)) = COS (x) *SINH(y> ; 

REAL (COS (z)) = COS(x)*COSH(y) 
and IMAG(COS(z)) = -SIN(x) *SINH(y) ; 

REAL(SINHCz)) = COS (y) *SINH(x) 
and IMAG(SINH(z)> = SIN(y) *COSH(x) ; 

REAL (COSH (z)) = COS(y)*COSH(x) 
and IMAG(COSH(z)) = SIN(y) *SINH(x) . 

To avoid making calls to evaluate SINH 
and COSH separately, and thus frequently 
having to evaluate EXP twice for the same 
argument, SINH(u) is computed as follows: 

1. u > 0.3465736 

SINH(u) = (EXP(u) - l/EXP(u))/2 

2. < u < 0.3465736 

SINH(u)/u is approximated by a polynomial 
of the fifth degree in u**2 which has a 
relative error of less than 2**-61.8 

3. u < 

SINH(u) = -SINH(-u). Then 

COSH(u) = SINH (ABS (u) ) + 1/EXP (ABS (u) ) . 



Effect of an Argument Error: 

Combine the effects on SIN, COS, SINH and 
COSH according to the method of evalua- 
tion described in the above paragraph. 



SIN, SINH, COS it COSH (long floating-point 
complex) 



Module Name: IHESNZ 

Entry Points: 

Mathematical 
function 

Sin z 

Hyperbolic sin z 
Cos z 
Hyperbolic cos z 



PL/I 
name 

SIN(z) 
SINH(Z) 
COS(z) 
COSH(Z) 



Entry 
point 

IHESNZS 
IHESNZZ 
IHESNZC 
IHESNZK 



Function: 

To calculate sin z or hyperbolic sin z, 
or cos z or hyperbolic cos z. 
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Accuracy: 



Arguments 



Relative Error 
♦10**15 



Range 



| Distribution | RMS | Maximum | 



IHESNZS 




._X_ 




—X- 










|ABS(x) 

1* 10,, 

|ABS(y) 
i __ 


< 


— T~ 

1 
1 

11 

± 


Uniform 


— T" 

1 

1 

1 
_ j. 


2, 


11 


"T — 

1 
1 
1 


82.4 | 


IHESNZZ 


|ABS(x) 

1* io, 

|ABS(y) 
t _ 


< 


T 

1 
1 

1| 

JL 


Uniform 


T 
1 
1 

1 

J. 


0. 


180 


T 
1 
1 

1 

_L 


2.31 \ 


IHESNZC 


|ABS(X) 

|s io, 

|ABS(y) 

L 


< 


1 

1 

1| 

x_ 


Uniform 


T 

1 

1 

1 
X. 





389 


T 

1 
1 
1 

_x 


6.33 | 


IHESNZK 



















|ABS(x) 
|£ 10„ 
|ABS(y) 
|< 20 

L 



| Uniform | 1.11 | 19.4 | 

I III 

I III 

I III 

X X X , J 



Error and Exceptional Conditions: 

O : IHESNZS, IHESNZC: 

ABS(x) £ 2**50*pi: error caused in 
real SIN routine (IHESNL) 

IHESNZZ,, IHESNZK: 

ABS(y) £ 2**50*pi: error caused in 
real SIN routine (IHESNL) 

H : OVERFLOW in real EXP routine (IHEEXL) 



Implementation : 

• Module size: 368 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/ 360 models given 
below are obtained from the appropriate 
entry point in the tables: 

SIN, COS : ABS(y) > 0.3465736 
SINK, COSH: ABS(x) > 0.3465736 



40 



■T 

| 50 



65 



■T 

I 75 



r t T" 

j Entry j 30 j 
| Point] | 
F + x + x x 

j IHESNZS | 46584 j 11267| 2363 | 552 | 313 



| IHESNZC | 46656| 11294 | 2373 | 555 | 315 



| IHESNZZ | 46726| 11317 | 2378 | 557 | 316 



|IHESNZK| 46640| 11304| 2374 | 556 | 316 

L X X X X X J 



SIN, COS : ABS(Y) < 0.3465736 
SINH,COSH: ABS(X) <, 0.3465736 

r T - -- T T T T 1 

j IHESNZS | 54173| 13141 j 2656 j 612 | 345 | 
h X X X X X ^ 

|IHESNZC| 54245| 13168J 2666 | 615 | 347 | 



H 



I IHESNZZ I 54325| 13191 | 2671 | 617 | 347 | 

|IHESNZK| 54247| 13177 | 2667 | 61S | 348 | 
L X X X X ; X J 



TAN,, TANH (short f loatinq-point complex) 



Module Name: IHETNW 






Entry Points: 








Mathematical 




PL/I 


Entry 


function 




name 


point 


Tan z 




TAN(z) 


IHETNWN 


Hyperbolic tan 


z 


TANH(Z) 


IHETNWH 



Function : 

To calculate tan z or hyperbolic tan z. 
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Method: 

Let z = x + yl. 

Then REAL < TAN (z)) = 

TAN(x)*(l - TANH(y)**2)/ 
(1 + (TAN(x)*TANH(y))**2),, 

IMAG(TAN(z)) = 

TANH(y)*(l + TAN(x)**2)/ 
(1 + (TAN(x)*TANH(y))**2) . 

TANH(z) = - (TAN(zD)I. 



Effect of an Argument Error: 

The absolute error caused in the result 
is approximately equal to the absolute 
error in the argument divided by 
AES(C0S(z)**2) for IHETNWN, or divided by 
ABS(C0SH(z)**2) for IHETNWH. The rela- 
tive error caused in the result is 
approximately twice the absolute error in 
the argument divided by AES(SIN(2*z) ) for 
IHETNWN, or divided by ABS(SINH(2*z) ) for 
IHETNWH. 



Accuracy: 



Arguments 



Range j Distribution! 



-t n 

| Relative Error | 
j *io**6 | 



RMS 



Maximum 



L . 

IHETNWN 




. X- 




._x. 




_l_. 




r _ _ _ . 

|ABS(x) 

|ABS(y) 

1 


< 
< 


— T~ 

11 
9| 


Uniform 


— T" 
1 
1 


0.430 


- T — 

1 
1 


1 . 67 | 


IHETNWH 


|ABS(x) 
|AES(y) 
L 


< 
< 


T 

91 
1| 

-_X_ 


Uniform 


T 

1 

1 
X. 


0.430 


T 

1 

1 
-1 


1.45 | 



Error and Exceptional Conditions: 

I : OVERFLOW 

: ABS(u) > 2**18*pi, where 
u = x for IHETNWN, 
u = y for IHETNWH. 



H 



OVERFLOW in real TAN routine (IHETNS) 



Implementation : 

• Module size: 184 bytes 

• Execution times : 

Approximate execution times in micro- 
seconds for the System/ 360 models given 
below are obtained from the formulas: 

IHETNWN: a + time for execution of 
IHETHS with argument y 

IHETNWH: b + time for execution of 
IHETHS with argument x 

r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
|. X + + + + .| 

| a j 9094 | 2310 | 696 | 186 | 111 | 
| b | 9197 | 2454 | 716 | 191 | 115 | 

L X X X X X J 



TAN, TANH (long f loatinq-point complex) 

Module Name: IHETNZ 
Entry Points: 



Mathematical 
function 

Tan z 
Hyperbolic tan z 



PL/I 
name 

TAN(z) 
TANH(z) 



Entry 
point 

IHETNZN 
IHETNZH 



Function: 

To calculate tan z or hyperbolic tan z. 

Method: 

Let z = x + yl. 

Then REAL ( TAN (z)) = 

TAN(x)*(l - TANH(y)**2)/ 
(1 + (TAN(x)*TANH(y))**2) , 

IMAG(TANCz)) = 

TANH(y)*(l + TAN(x)**2)/ 
(1 + (TAN(x)*TANH(y))**2). 

TANH(z) = - (TAN(zD)I. 
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Accuracy: 



ATAN, ATANH (short f loatinq-point complex) 



IHETNZN 



Arguments 



Relative Error 
*10**15 



t • T t — • 1 

Range | Distribution | RMS | Maximum | 
J. . x x , J 



ABS(x) < 1| Uniform | 0.139 | 
ABS(y) < 9 | J | 
x x . x. 



IHETNZH 



- T T * T 

ABS(x) < 9 | Uniform j 0.137 j 
ABS(y) < 1| | | 
x x . x. 



1.11 | 



0.980 | 



Module Name: IHEATW 

Entry Points: 

Mathematical 
function 

Arctan z 
Hyperbolic arctan z 

Function: 



To calculate arctan z or hyperbolic arc- 
tan z. 



PL/I 

name 

ATAN(z) 
ATANH (z) 



Entry 
point 

IHEATWN 
IHEATWH 



Error and EXCEPTIONAL Conditions: 

I : OVERFLOW 

O : ABS(u) > 2**50*pi,, where 
u = x for IHETNZN, 
u = y for IHETNZH. 

H : OVERFLOW in real TAN routine (IHETNL) 

Implementation : 

• Module size: 184 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table: 

r T ~T T •' T. 'T 1 

|ABS(y) | 30 j 40 | 50 j 65 | 75 j 
l x x x x ,x J 

IHETNZN 

r T T — ■ "T T T 1 

|< 0.549| 40843J 9625 j 1844 j 411 j 233 j 

j. x x x x _ + H 

|0.549 <| | | | | | 
|ABS(y) | 44498J10513 J10365 j 493 j 284 j 

l< 20.1 | | | | | | 
|. + + x x x 1 

|> 20.1 | 29337| 6967 | 1415 | 327 | 190 | 

L X_ X . X <_X ,X J 

|ABS(x) I 30 I 40 I 50 I 65 I 75 I 

L X X X X .X J 

IHETNZH 

|< 0.549| 41122| 9709 | 1871 | 419 | 236 | 
|. x__- x_„ x +— ~- + ^ 

|0.549 S| | | | | | 
|ABS(x) | 44777|10597 | 2185 | 501 | 287 
l< 20.1 | | | I I I 

L i i___ X .1 _„X J 

|2: 20.1 | 29616| 7051 | 1442 | 334 | 193 | 
t x x x x— x J 



Method: 

Let z = x + yl. 

Then REAL(ATANH(z)) = 

IMAG (ATANH (z)) = 



(ATANH(2*x/ 
(l+x*x+y*y)))/2 

(ATAN(2*y, 
(l-x*x-y*y)))/2 



and ATAN(z) 



= -(ATANH(zD)I. 



Effect of an Argument Error: 

The absolute error in the result is 
approximately equal to the absolute error 
in the argument divided by (1 + z**2) in 
the case of IHEATWN, or by (1 - z**2) in 
the case of IHEATWH. Thus the effect may 
be considerable in the vicinity of 
z = ± II (IHEATWN) or ± 1 (IHEATWH). 



Accuracy: 

r ■ ■ t 1 

| Arguments | Relative Error | 
j j *10**6 j 
|. T x T ^ 

] Range [Distribution | RMS | Maximum | 

L X X X J 

IHEATWN 

r T T T 1 

II II I 

| Full range | Exponential j 0.216 j 2.88 j 
l x x x J 

IHEATWH 

r T . — T ¥ 1 

j Full range | Exponential j 0.208 | 1.18 | 

L X - X X J 



j Error and Exceptional Conditions: 

P : IHEATWN: z + ±11 
IHEATWH: z = ±1 
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Implementation : 

• Module size: 304 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table, where 

a = ABS(2*u/(l+x*x+y*y)) 
u = y for IHEATWN 
= x for IHEATWH 

r — T T T T T T 

| | 30 | 40 | 50 j 65 j 75 J 
t x x x x x J 

IHEATWN 

|a < 0.2| 12235| 3306 | 1033 | 279 | 173 | 
j. x x x x x .| 

1 0.2 < a J 16056J 4454 | 1408 j 462 j 276 j 

l< 1 I I I II I 

L X X X X X J 

IHEATWH 

fa < 0.2| 12100| 3267 | 1017 | 275 | 171 | 
j. x + x x 1 .| 

| 0.2 < a| 15921J 4415 | 1392 | 458 | 274 j 

l< 1 I I I II I 

L X X X X X J 



Effect of an Argument Error: 

The absolute error in the result is 
approximately equal to the absolute error 
in the argument divided by (1 + z**2) in 
the case of IHEATZN, or by (1 - z**2) in 
the case of IHEATZH. Thus the effect may 
be considerable in the vicinity of 
Z = ± I (IHEATZN) or ± 1 (IHEATZH). 



Accuracy: 



T 1 

Arguments | Relative Error | 
I *1Q**15 I 



r t 1 t 

I Range | Distribution | RMS | Maximum 



H 



IHEATZN 






r t 

|Full range (Exponential 
L x 


| 0.141 | 


7.93 | 


IHEATZH 


r _ _ T 

|Full range (Exponential 

L X 


| 0.0826| 
.x_ x_. 


1.20 | 



Error and Exceptional Conditions : 

P : IHEATZN: z = ±11 
IHEATZH: z = ±1 



ATAN, ATANH (long floating-point complex) 



Module Name: IHEATZ 






Entry Points: 






Mathematical 


PL/ 1 


Entry 


function 


name 


point 


Arctan z 


ATAN(z) 


IHEATZN 


Hyperbolic arctan z 


ATANH (z) 


IHEATZH 



Function: 

To calculate arctan z or hyperbolic arc- 
tan z. 

Method : 

Let z = x + yl. 

Then REAL (ATANH (z)) = (ATANH(2*x/ 

(l+x*x+y*y)))/2 

IMAG (ATANH (z) ) = (ATAN(2*y, 

(l-x*x-y*y)))/2 



Implementation : 

• Module size: 296 bytes 

• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
entry point in the table, where 

a = ABS(2*u/(l+x*x+y*y)) 
u = y for IHEATZN 
= x for IHEATZH 

r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 

L X X . X X X J 

IHEATZN 

|a£0.25 | 43477| 9977 | 2006 | 455 | 260 | 

|0.25<a | 51673|11840 | 2406 | 576 | 329 | 
l< 1 I I I I I I 

L X X J X X J 

IHEATZH 

|a£0.25 | 43293| 9923 | 1987 | 450 | 258 | 
|. x x x x x .| 

|0.25<a | 51489|11786 | 2466 | 571 | 327 | 

l< 1 I I I II I 

L X X X X X J 



and ATANfz) 



= -(ATANH(zD)I. 
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CHAPTER Hi ARRAY FUNCTIONS 



The Library supports the array built-in 
functions SUM, PROD, f POLY, ALL and ANY, and 
also provides indexing routines for han- 
dling simple (i.e., consecutively stored) 
and interleaved arrays. 



Input Data 



The array function modules are distingu- 
ished from the other Library modules in 
that they all accept array arguments and 
perform their own indexing, whereas the 
other modules require that indexing should 
be handled by compiled code. Calls to 
conversion routines are included in the 
SUN,, PROD and POLY modules with fixed-point 
arguments, so that these arguments are 
converted to floating-point as they are 
accessed (it should be noted that it is a 
requirement of the language that the 
results from these modules be in 
floating-point) . On the other hand, the 
conversions necessary for the ALL and ANY 
modules (the arguments must be converted to 
bit string arrays) are not part of the 
modules and must be carried out before the 
modules are invoked. 

Any restrictions on the admissibility of 
arguments are noted under the headings 
•Range' and 'Error and Exceptional 
Conditions* . 

Range ; This states any ranges of arguments 
which a module assumes to have been exclud- 
ed prior to its being called. 

Error and Exceptional Conditions ;, These 
cover conditions which may result from the 
use of a routine; they are listed in four 
categories: 

P — Programmed conditions in the module 
concerned. Programmed tests are 
made where this is not too costly 
and, if an invalid argument is 
found, a branch is taken to the 
entry point IHEERRC of the execution 



error package (EXEP) . This results 
in the printing of an appropriate 
message and in the ERROR condition 
being raised. 

I — Interrupt conditions in the module 
concerned. For those routines where 
SIZE and FIXEDOVERFLOW are detected 
by programmed tests or where hard- 
ware interruptions may occur, the 
OVERFLOW, UNDERFLOW, and (when the 
conversion package is called) SIZE 
conditions pass to the ON handler 
(IHEERR) and are treated in the 
normal way. The machine is assumed 
to be enabled for all interruptions 
except significance, which is masked 
off. 

O — Programmed conditions in modules 
called by the module concerned. 
These occur when invalid arguments 
are detected in the module called. 

H -- As I, but the interrupt conditions 
occur in the modules called by the 
module concerned. 



Effect of Hexadecimal Truncation 



See the corresponding section in the 
introduction to Chapter 3 for guidance to 
the accuracy of SUM, PROD, and POLY. If 
fixed-point arguments are passed to these 
functions, further errors may be introduced 
by conversions. 



Speed 



The average execution times given are 
based on IBM System/360 Instruction Timing 
Information , Form A22-6825. 

A summary of the Library array modules 
is given in Figures 6 and 7. 
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"T T" 

| Simple arrays, and | 
| interleaved arrays of | 
j variable length strings | 



h 

| Indexers | 
| ALL,, ANY | 

L J.. 



Interleaved string j 

arrays with fixed- j 

length elements j 



IHEJXS 
IHENL1 



IHEJXI 
IHENL2 



Note: IHEJXI is used for 
arithmetic arrays 



indexing through interleaved 



Figure 6. Bit String Array Functions and Array Indexers 



PL/ 1 
function 



Fixed-point 
arguments 



Floating-point arguments 

T 

Short precision | Long precision 



F t +- 

| Simple | Interleaved | 



t T t- 1 

Simple | Interleaved! Simple | Interleaved! 



SUM real | IHESSF 
complex | IHESSX 

I 
PROD real j IHEPSF 
complex j IHEPSX 

j. x— 

POLY real j IHEYGF 
complex j IHEYGX 



IHESMF | 
IHESMX | 

I 

IHEPDF j 
IHEPDX j 



IHESSG 
IHESSG 

IHEPSS 
IHEPSW 



IHESMG 
IHESMG 

IHEPDS 
IHEPDW 



| IHESSH 
| IHESSH 

I 

j IHEPSL 
j IHEPSZ 



IHESMH | 
IHESMH | 

I 
IHEPDL | 
IHEPDZ | 



IHEYGS 
IHEYGW 



IHEYGL 
IHEYGZ 



Figure 7. Arithmetic Array Functions 



ARRAY INDEXERS 

Indexer for Simple Arrays 

Module Name: IHEJXS 

Entry Points: 

Element Entry 

address point 

Bit addresses IHEJXSI 
Byte addresses IHEJXSY 

Function: 

To find the first and last elements of an 
array. Their addresses are returned, in 
general registers and 1 respectively,, 
as bit addresses (IHEJXSI) or byte 
addresses (IHEJXSY). 



Method: 

The address of the virtual origin B of 
the array (i.e., the address that would 
correspond to the element A(0, ..0)) is 
obtained as a byte address for IHEJXSY, 
or a bit address for IHEJXSI, by refer- 
ring to the first word of the array dope 
vector (adv) . 



Address of first element = B + ^ M ±^i 

i=l 
n 

Address of last element = B + ^ H^Pi 

i=l 

where M-^ is the multiplier for the ith 
dimension 

L^ is the lower bound for the ith 
dimension 

Ujl is the upper bound for the ith 
dimension, and 

n is the number of dimensions. 



Range : 

< number of dimensions < 2**22 
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Implementation : 

• Module size: 104 bytes 

• Execution time: 

Let n = the number of dimensions 

Then the approximate execution tiroes in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

IHEJXSI: a + c*n 

IHEJXSY: b + c*n 

r T T 'T T 1 1 

| | 30 | 40 | 50 | 65 | 75 | 
j. + + _ + + + ., 

| a | 720 | 229 | 85.8 | 22.0 | 13.4 | 
| b | 555 | 183 | 71.3 | 18.2 | 11.4 | 
| c | 377 | 148 | 64.8 j 16.0 | 10.0 | 

L X X X X X J 



Indexer for Interleaved Arrays 

Module Name: IHEJXI 
Entry Points: 



Operation 

Initialization for bit 

addresses 
Initialization for byte 

addresses 
Elements after the first 



Entry 
point 

IHEJXII 

IHEJXIY 

IHEJXIA 



Method : 

Arrays are stored in row major order. 
Let Li be the lower bound and Ui the 
upper bound of the ith dimension, and n 
the number of dimensions. Starting with 
the element A(L if L 2f L n ) , the rou- 
tine varies the subscripts through their 

ranges to A(U ± ,U 2 , U n ) , changing the 

nth subscript most rapidly; in this way 
the elements are referenced in the order 
in which they are stored. 

The routine does not deal with actual 
subscript values but calculates the 
extent Ej_(= Ui - L^ + 1) of each dimen- 
sion and uses this as a count that varies 
from Ej_ to 1 for subscript values Li to 
Ui. A "base address' for each dimension 
is maintained and, for the ith dimension, 
is defined as the address of the element 
with ith subscript equal to its lowest 
bound Li and with all other subscripts at 
their current values. 

Thus initially the base addresses are all 
equal to the address of 
A(Li,L 2 , L n ) • Each subsequent ele- 
ment address is generated from the pre- 
vious one by adding the multiplier M n 
from the array dope vector (ADV) and 
reducing the subscript count by 1. When 
the count for the ith dimension has been 
reduced from Ei to 1 it is reset to Eij, 
M i-1 i s added to the (i-l)th dimension's 
base address and the count for this 
dimension is decreased by one. 

This new base is the starting point for 
further increments by M n . When a new 
base address is calculated, the base 
addresses for all higher dimensions 

((i+1), (i+2), n) is set equal to 

the ith base address. 



Function: 

To find the next element of an array and 
to return its bit or byte address in 
general register i. 

Entry point IHEJXII is used to initialize 
the routine for bit addresses and to 
provide the address of the first element 
in the array; IHEJXIY does the same for 
byte addresses. Entry point IHEJXIA is 
used thereafter to obtain the addresses 
of subsequent elements of the array; one 
address is returned for each entry into 
the routine. 



Range: 

< number of dimensions < 2**22 
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Implementation : 

• Module size: 328 bytes 

• Execution time: 

Let n = the number of dimensions of 
the array 

E ± = Oi-Ii+1 

n-1 

T = 7T Ej_ 

i=l 



n-x n-x 

= E 7T Ej 

i=l i=l J 



n^l n-i 
7 

y 

V x = time to get a VDA (IHESADF) 

V a = time to free a VDA (IHESAFD) 

Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 



IHEJXII 
IHEJXIY 
IHEJXIA: 
n-1: 
n > 1: 



a + V A + c*n 
b + V ± + c*n 

d + V 2 + e*E ± 

k + V 2 + h*S + 

T*(f + e*E n + g*n*(n-D) 



The total time required to index 
through the complete array is the sum 
of: 

i) the time for either IHEJXII or 
IHEJXIY, and 

ii) the time from the appropriate 
IHEJXIA formula 



1 


1 


30 


1 


40 


1 


50 


1 


65 


1 


75 


1— 


-+- 




-+- 




-+- 




-+- 




~H 





1 a 


1620 | 


1 fc 


1398 | 


1 c 


508 1 


1 d 


122 | 


1 e 


786 | 


1 f 


204 | 


1 g 


42 | 


1 h 


260 | 


1 fc 


220 | 



526 


203 


55.9 


472 


182 


50.8 


192 


75.8 


20.6 


51.5 


8.8 


2.5 


267 


107 


28.8 


66.7 


25.8 


7.7 


14.4 


4.8 


1.4 


95.3 


36.8 


11.1 


91.9 


24.0 


7.6 



34.4 

31.8 

14.3 

1.6 

19.7 

5.3 

1.1 

8.2 

5.3 



Note.: the fastest overall execution time 
will occur when the extent of the 
nth dimension is the largest of the 
subscript extents of the array 
dimensions. 



ARRAY FUNCTIONS 



AIL (X), ANY (X) 



Module Names: 



Arguments 

Simple arrays and interleaved 
arrays with variable-length 
elements 

Interleaved arrays with fixed- 
length elements 

Entry Points: 



PL/I function 



AIL(X) 1# ANY(X), byte-aligned 



ALL(X), any alignment 



ANY(X), any alignment 



Function: 



Module 
name 

IHENL1 



IHENL2 



Entry 
point 

IHENL1A 
IHENL2A 

IHENL1L 
IHENL2L 

IHENL1N 
IHELN2N 



The argument X is a bit string array (any 
necessary conversion having been per- 
formed prior to the invocation of these 
modules) . The result is a scalar bit 
string of length equal to the greatest of 
the current lengths of the elements of X. 

ALL(X) : the ith bit of the result is 1 if 
the ith bits of all the elements of X 
exist and are 1; otherwise it is 0. 

ANY(X): the ith bit of the result is 1 if 
the ith bit of any element of X exists 
and is 1; otherwise it is 0. 

Method: 

For byte-aligned string arrays, AND 
(IHEBSA) and OR (IHEBSO) are used for ALL 
and ANY respectively; for string arrays 
with any alignment BOOL (IHEBSF) is used 
with appropriate parameter bits. 

The elements of the array are passed to 
IHEBSA, IHEBSO or IHEBSF one at a time, 
and the result is developed in the target 
field. For the first call to any of 
these logical modules the first element 
of the array serves as both first and 
second source arguments. For subsequent 
calls, the result already developed in 
the target field is the first argument 
and the next element of the array is the 
second argument. 
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Range: 



SUM (X) 



Bit strings are limited to a maximum of 
32,767 bits. 



Implementation : 

• Module size: IHENLl: 280 bytes 

IHENL2: 192 bytes 

• Execution time: 

Let R = number of elements in the 
array 

T x = time required to execute 
appropriate bit string routine 
via IHEBSAO, IHEBSO0 or 
IHEBSFO 

T a = time to index via IHEJXSI 

T 3 = time to index via IHEJXSY 

1\ ■= sum of times required to index 
via IHEJXII and IHEJXIA 



Then the approximate execution times in 
microseconds for the System/360 models 
given below are obtained from the fol- 
lowing formulas: 

IHENLl 

Fixed-length elements: 

a + R*(b + T x + T 2 ) 

Varying-length elements: 

a + R*(c + T x + T 3 ) 

r T T — T T 'T 1 

I I 30 I 40 | 50 | 65 | 75 j 
j. + + + x -X .J 

| a | 1318 | 470 | 189 | 51.6 | 36.9 | 
j b | 837 | 258 | 95 | 27.3 | 18.0 j 
| c | 375 | 134 | 49.8 | 15.2 | 10.7 | 

L X X ._X X ,X J 

IHENL2 

a + T n + R*(b + T x ) 

r T T T T "T 1 

| | 30 | 40 | 50 | 65 | 75 | 

j. x _+ x + „ „x ) 

| a | 1341 | 479 | 188 | 53.5 | 36.6 | 
| b | 428 | 133 | 54.0 | 15.2 | 11.5 | 

L X X . X X- -X J 



Module Names and Entry Points: 
Simple Arrays 





Module 


Entry 


Arquments 


name 


point 


Fixed, real 


IHESSF 


IHESSF0 


Fixed, complex 


IHESSX 


IHESSX0 


Short float 






real 


IHESSG 


IHESSGR 


complex 


IHESSG 


IHESSGC 


Long float 






real 


IHESSH 


IHESSHR 


complex 


IHESSH 


IHESSHC 


Interleaved Arrays 




Module 


Entry 


Arquments 


name 


point 


Fixed, real 


IHESMF 


IHESMF0 


Fixed, complex 


IHESMX 


IHESMX0 


Short float 






real 


IHESMG 


IHESMGR 


complex 


IHESMG 


IHESMGC 


Long float 






real 


IHESMH 


IHESMHR 


complex 


IHESMH 


IHESMHC 



Function: 

To produce a scalar with a value which is 
the sum of all the elements of the array 
argument. 



Method: 

The elements of the array are added to 
the current sum in row major order. 

For fixed-point arguments each element is 
converted to floating-point by using the 
PL/I Library conversion package. 

For a complex argument, the summation of 
the real parts is performed before the 
summation of the imaginary parts is begun 
in modules IHESSG and IHESSH, while the 
two sums are developed concurrently in 
other modules. 



Error and Exceptional Conditions: 

I : OVERFLOW, UNDERFLOW 

H : IHESSF, IHESSX, IHESMF, IHESMX: 

ABS (element of the array) > 
7.2*10**75: SIZE condition caused in 
conversion package 
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Implementation : 
• Module Sizes: 



Nodule 


Bytes 


IHESSF 


168 


IHESSX 


216 


IHESSG 


104 


IHESSH 


104 


IHESMF 


136 


IHEStfX 


224 


IHESNG 


128 


IHESMH 


128 



75 



• Execution times: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
formula. 



IHESSX 

Source 

binary 

decimal 

binary 

decimal 



Target 

short a + T 2 + R*(e+2*T 3 ) 

short b + T a + R*(e+2*T 3 ) 

long c + T a + R*(f+2*T 3 ) 

long d + T 2 + R*(f+2*T 3 ) 



r— 

1 


1 


30 


1 


40 


— T _ 

1 


50 


1 


65 


1 


75 


1 


f- 


-+- 




-+- 




-+- 




-+- 




-+- 




H 



a 


1309 


492 


208 | 


54.9 


40.7 | 


b 


1395 


515 


216 | 


58.9 


41.3 | 


c 


1385 


518 


315 | 


58.0 


41.5 | 


d 


1471 


540 


224 - | 


60.0 


42.0 | 


e 


776 


259 


94.5 | 


24.0 


15.2 | 


f 


834 


267 


96.1 1 


22.9 


14.1 I 







Constants used in these formulas are: 

R = number of elements in the array 

T x = sum of times required to execute 
IHEJXI using IHEJXIY and IHEJXIA 

T 2 = time to execute IHEJXS by means of 
IHEJXSY 

T 3 = time for the appropriate conver- 
sion using IHEDMA 

The binary and decimal source data is 
always fixed-point; target data is 
short or long floating-point. 



IHESSQ 

Real a + T 2 + R*b 

Complex c + T 2 + R*d 

r T T T T T 

| | 30 | 40 I 50 | 65 | 75 



| a | 887 | 301 | 118 | 

| b | 127 | 34.9 | 12.4 | 

| c | 1033 I 347 J 134 j 

| d | 254 | 67.8 | 24.8 | 

L X X X X. 



32.0 | 21.4 | 
4.1 | 2.1 

36.7 | 24.6 
8.3 | 4.2 

x 



IHESSF 

Source 

binary 

decimal 

binary 

decimal 



Target 

short a + T a + R*(e + T 3 ) 

short b + T a + R*(e + T 3 ) 

long c + T 2 + R*(f + T 3 ) 

long d + T a + R*(f + T 3 ) 



1 


30 


1 


40 


1 


50 


1 


65 


1 


75 


1 


X- 




-+- 




-+- 




-+- 




"+- 




H 



a 


1281 


446 | 177 | 


46.6 


31.6 | 


b 


1367 


469 | 186 | 


48.7 


32.2 | 


c 


1308 


454 | 179 | 


46.6 


31.4 | 


d 


1394 


477 | 188 | 


48.6 


32.0 | 


e 


119 


35.5 | 13.4 | 


12.4 


8.3 | 


f 


137 


35.2 | 12.2 | 


11.2 


7.2 | 




L J 


L X x_ 







IHESSH 

Real a + T 2 + R*b 

Complex c + T a + R*d 

r t t t t t 1 

| | 30 | 40 | 50 | 65 | 75 | 
j. + x x _ x x 1 

| a | 935 | 314 | 121 | 32.0 | 21.4 | 
| b | 167 | 43.9 | 15.2 | 4.2 | 2.1 | 
| c | 1129 | 372 | 142 | 36.8 | 24.6 | 
| d | 334 | 87.8 | 30.4 | 8.3 \ 4.2 | 
l x x x x x J 
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IHESMF 

Target 

short a + T ± + R*(e + T 3 ) 

long c + T ± + R*(f + T 3 ) 

r T T -T T— T 1 

I I 30 I 40 | 50 | 65 I 75 | 
J. x— + _ + + (. 1 

| a | 1074 | 363 | 141 | 37.3 | 24.6 | 
| c | 1123 | 381 | 147 | 38.5 | 25.6 | 
| e | 400 | 139 | 54.4 | 16.2 | 10.9 | 
| f | 418 | 139 | 53.2 | 15.1 | 9.8 | 
l X X X X X J 



IHESMH 

Real a + T ± + R*b 

Complex c + T ± + R*d 

r t t t t t T 

| | 30 | 40 | 50 | 65 | 75 | 



I a | 
I b I 
I c | 

I a | 

L X. 



887 
366 j 
993 | 
514 | 
x. 



298 
116 
337 
157 



| 111 | 
| 42.4 | 
I 125 | 
| 55.9 | 
-X x. 



30.9 | 20.4 | 

12.3 | 7.9 | 

34.4 j 23.5 | 
16.0 | 9.8 | 

x J 



IHESMX 




Source 


Tarqet 


binary 


short 


decimal 


short 


binary 


long 


decimal 


long 



a + T ± + R*(e+2*T 3 ) 
b + T x + R*(e+2*T 3 ) 
c + Ti + R*(f+2*T 3 ) 
d + T x + R*(f+2*T 3 ) 



1 


1 


30 


1 


40 


1 


50 


1 


65 


1 


75 


1 


K- 


--+- 




-+- 




-+- 




-+- 




--+- 




H 



1 a 


1327 | 


1 b 


1413 | 


1 c 


1319 | 


1 cl 


1405 | 


1 e 


712 | 


| f 


770 | 



511 


229 


63.8 


533 


238 


65.9 


515 


227 


63.7 


537 


235 


65.8 


232 


87 


25.1 


240 


89 


24.0 



47.8 
48.4 
47.4 
48.0 
16.4 
15.2 



IHESMG 

Real a + T ± + R*b 

Complex c + Tj. + R*d 

r T T T T T 1 

I I 30 j 40 | 50 | 65 | 75 | 
j. X X -X X + ^ 

| a | 839 | 286 | 107 | 30.3 | 20.4 | 
| b | 294 | 97.4 | 35.6 | 11.9 | 7.9 | 
I c | 929 j 320 | 119 | 33.9 | 23.5 | 
| d j 402 | 129 | 46.3 | 15.6 | 9.8 | 

L X X .-X X X J 



PROD (X) 

Module Names and Entry Points: 
Simple Arrays 



Arguments 



Module 
name 



Entry 
point 



Fixed, real 


IHEPSF 


IHEPSF0 


Fixed complex 


IHEPSX 


IHEPSX0 


Short float 






real 


IHEPSS 


IHEPSS0 


complex 


IHEPSW 


IHEPSW0 


Long float 






real 


IHEPSL 


IHEPSL0 


complex 


IHEPSZ 


IHEPSZ0 


Interleaved Arrays 




Module 


Entry 


Arguments 


name 


point 



Fixed, real IHEPDF IHEPDF0 

Fixed, complex IHEPDX IHEPDX0 

Short float 

real IHEPDS IHEPDS0 

complex IHEPDW IHEPDW0 

Long float 

real IHEPDL IHEPDL0 

complex IHEPDZ IHEPDZ0 

Function: 

To produce a scalar with a value which is 
the product of all the elements in the 
array argument. 

Method : 

The elements of the array are used in row 
major order to multiply the current pro- 
duct. 

For fixed-point arguments, each element 
is converted to floating-point by using 
the PL/I Library conversion package. 



Chapter 4. Arrays: Array Functions 75 



Error and Exceptional Conditions: 

I : OVERFLOW, UNDERFLOW 

H : IHEPSF, IHEPSX, IHEPDF, IHEPDX: 

ABS(eleirent of the array) > 
7.2*10**75: SIZE condition caused in 
conversion package 



Implementation : 




• Module sizes: 




Module 


Bytes 


IHEPSF 


160 


IHEPSS 


72 


IHEPSL 


72 


IHEPSX 


256 


IHEPSW 


96 


IHEPSZ 


96 


IHEPDF 


144 


IHEPDS 


88 


IHEPDL 


88 


IHEPDX 


272 


IHEPEW 


120 


IHEPDZ 


120 



• Execution tiroes: 

Approximate execution times in micro- 
seconds for the System/360 models given 
below are obtained from the appropriate 
formula. 



IHEPSF 




Source 


Tarqet 


binary 


short 


decimal 


short 


binary 


long 


decimal 


long. 



a + T 2 + R*(e+T 3 ) 
b + T 2 + R*(e+T 3 ) 
c + T 2 + R*(f+T 3 ) 
d + T 2 + R*(f+T 3 ) 



«■-■ 

1 


— T 

1 


30 


T - 
1 


40 


1 


50 


_ T 
1 


65 


- T 

1 


75 


F- 


"+- 




"+- 




--+- 




-+- 




"+- 





a 


1208 


416 


163 | 


43.8 


29.3 | 


b 


1294 


439 


172 | 


45.9 


29.9 | 


c 


1257 


434 


169 | 


45.0 


30.2 | 


d 


1343 


456 


178 | 


47.1 


30.8 | 


e 


613 


18 8 


62.3 | 


16.1 


10.3 | 


f 


1331 


357 


74.8 I 


18.1 


11.2 1 


X X X x_ 


X J 



IHEPSS 

a + T a + R*b 

r T T T T T 1 

| J 30 | 40 | 50 | 65 | 75 | 

|. x x x x + H 

| a | 456 | 191 | 82.3 | 24.6 | 17.7 | 
| b | 372 | 96.9 | 27.0 | 5.9 | 3.3 | 
l. X X X X X J 



Constants used in the formulas are: 

R = number of elements in the array 

Ta. = sum of times required to execute 
IHEJXI using IHEJXIY and IHEJXIA 

T 2 = time to execute IHEJXS via 
IHEJXSY 

T 3 = time for the appropriate conver- 
sion using IHEDMA 

The binary and decimal source data is 
always fixed-point; target data is 
short or long floating-point. 



IHEPSL 

a + T 2 " + R*b 

r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
|. x x x x x -J 

| a | -252 j 22.5 j 41.5 j 16.0 j 12.9 j 
| b | 1112 | 276 | 43.5 | 9.1 | 5.3 | 

L X X X X X J 
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IHEPSX 

Source 

binary 

decimal 

binary 

decimal 



Target 

short a + T 2 + R*(e+2*T 3 ) 

short b + T 2 + R* (e+2*T 3 ) 

long c + T 2 + R*(f+2*T 3 ) 

long d + T 2 + R*(f+2*T 3 ) 



30 



40 



50 



65 



I 75 | 



Y — + +- 




-+- 




-+- 


X ., 


| a | 1405 | 


494 




192 




52.1 | 35.5 | 


| b | 1491 | 


517 




200 




54.1 | 36.0 | 


| c | 1481 | 


520 




199 




53.2 | 36.2 | 


| d | 1567 | 


542 




208 




55.2 | 36.8 | 


| e | 1993 | 


553 




173 




41.6 | 24.1 | 


| f | 5043 | 


1285 




241 




54.1 | 31.0 | 






_.x_ 




...x_ 





IHEPDS 

a + T x + R*b 

r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
(. x x x x x ^ 

| a | 418 | 178 | 80.0 | 23.7 | 17.5 | 
| b | 492 | 143 | 43.3 | 11.4 | 7.2 | 

L X X — . X X X J 

IHEPDL 

a + Ti + R*b 
r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
,. x x x x x .| 

| a J -290 | 8.8 j 67 . 5 j 20.7 | 15.5 | 
I b j 1264 | 331 | 63.8 | 14.8 | 9.2 | 
l x x x x x J 



IHEPSW 



a + T 2 + R*b 



r — t- 

1 1 

L X- 


T~ 

30 | 
_x_ 


40 | 

x_ 


50 | 

X— 


T ~ T 

65 || 75 j 

JL_ _ J 


r — t 
1 a | 


-614 | 


-69.8 | 


T 

7.7 | 


f- - "J 

8.9 j 10.6 j 


1 b | 


1508 | 


382 | 


108 | 


24.1 1 .12.1 | 


L X- 




. X. 




_ X -J 



IHEPSZ 



a + T 2 + R*b 



1 t 30 | 

L X -X- 


40 


1 
I 


50 | 

. L 


65 j 

X- 


75 | 
j 


r — T T 

| a |-3622 | 


-791 


1 
1 


-56.2 | 


f- 

-4.3 j 


2.6 | 


| b | 4580 | 


1124 


1 


180 | 


37.7 | 


20.1 | 


L X _X. 






. X 


— X- 





IHEPDX 

Source 

binary 

decimal 

binary 

decimal 



Target 

short a + Ti. + R*(e+2*T 3 ) 

short b + Tj. + R*(e+2*T 3 ) 

long c + T ± + R*(f+2*T 3 ) 

long d + T ± + R*(f+2*T 3 ) 



| | 30 | 40 | 50 | 65 | 75 | 



1382 
1468 
1374 
1460 
2047 
5097 



497 
520 
501 
523 
574 
1306 



208 
216 
205 
214 
18 2 
25 



56.8 
58.9 
56.7 
58.8 
44.4 
56.9 



41.2 
41.8 
40.8 
41.4 
26.3 
33.2 



IHEPDF 

Target 

short a + Ta. + R*(e + T 3 ) 

long c + To. + R*(f + T 3 ) 

r T T T T T 1 

| | 30 | 40 | 50 | 65 | 75 | 
|. + x x x x , .| 

| a | 1075 | 365 | 141 | 37.8 | 25.0 | 

| c | 1124 | 382 | 147 | 39.0 | 25.9 | 

| e | 645 | 201 | 69.0 | 18.2 | 12.1 | 

| f | 1363 | 371 | 81.5 | 20.2 | 13.0 | 

L X X .-X X X J 



IHEPDW 



a + Ti + R*b 















1 1 


30 | 


40 


1 


50 | 


65 | 75 | 


j. X. 


+- 




-+- 


x_ 


X ^ 


1 a | 


-814 | 


-96 


1 


12.1 | 


8.9 | 12.5 | 


1 b | 


1694 | 


452 


1 


132 | 


32.0 | 17.7 | 


L _X. 
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IHEPDZ 



a + T x ■+ R*b 



r T T T T T 1 

I I 30 I 40 J 50 J 65 | 75 j 
|. + + 1 + + ^ 

| a |-3841 | -852 | -61.5 | -4.7 | 4.2 | 
| t | 4830 | 1214 | 211 | 46.1 j 25.7 j 

L JL J. J. J. J. J 



POLY (A,,X) 



Module Names and Entry Points: 





Module 


Entry 


Arquments 


name 


point 


Fixed, real 






vector X 


IHEYGF 


IHEYGFV 


scalar X 


IHEYGF 


IHEYGFS 


Fixed,, complex 






vector X 


IHEYGX 


IHEYGXV 


scalar X 


IHEYGX 


IHEYGXS 


Short float,, real 






vector X 


IHEYGS 


IHEYGSV 


scalar X 


IHEYGS 


IHEYGSS 


Short float, complex 




vector X 


IHEYGW 


IHEYGWV 


scalar X 


IHEYGW 


IHEYGWS 


Long float, real 






vector X 


IHEYGL 


IHEYGLV 


scalar X 


IHEYGL 


IHEYGLS 


Long float, complex 




vector X 


IHEYGZ 


IHEYGZV 


scalar X 


IHEYGZ 


IHEYGZS 


Function: 






Vector X: 







Let the arguments he arrays declared as 
A(m:n) and X(p:q). Then the function 
computed is: 

n-m j-1 
A(m) + J2 A(m + j) * [77" X(p + i) 

j=l i=0 

unless n = m, when result is A(m). 

If q-p<n-m-l, then, for 
p + i > q, X(p + i) = X(q) . 

Scalar X: 

This may fce interpreted as a special case 

of vector X,, that is, a vector with one 

element, X(l)„ which is equal to X. Then 
the function computed is: 

n-m 

X) A(m + j)*X**j 
j=0 



A floating-point result is obtained in 
both cases. 

Method: 

1. Vector X, (q-p^n-m-1): 

POLY (A, X) is evaluated by nested multi- 
plication and addition, i.e., 

(.. . (A(n)*X(k) + A(n-l))*X(k-l) + 

A(n-2))* ... + A(m+l))*X(p) + A(m) 

where k = p + n-m- 1. 

2. Vector X, (q - p < n - m - 1) : 

In the expression above, the terms in X 
with subscript ranging from k down to q 
are all made equal to X(q) . The evalua- 
tion is treated as for scalar X until 
sufficient terms in X have been made 
equal to X(q), when the computation con- 
tinues as in (1. ) . 

3. Scalar X: 

Terms in X with subscript ranging from k 
to p are equal to X. 

For fixed-point arguments each element is 
converted to floating-point, by using the 
PL/I Library conversion package. 

Error and Exceptional Conditions: 

I : OVERFLOW, UNDERFLOW 

H : IHEYGF, IHEYGX: 

ABS( element of the array) > 
7.2*10**75: SIZE condition caused in 
conversion package 



Implementation : 




• Module sizes: 




Module 


Bytes 


IHEYGF 


432 


IHEYGS 


240 


IHEYGL 


240 


IHEYGX 


688 


IHEYGW 


280 


IHEYGZ 


280 



• Execution times: 

Let the arguments be declared as A(m:n) 
and X(p:q), or X, and T be the time for 
one conversion using the arithmetic 
conversion director IHEDMA. Then the 
approximate execution times in micro- 
seconds for the System/360 models given 
are obtained from the appropriate for- 
mula. 'Short' or 'long' refers to the 
floating-point result. 
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IHEYGF 

Scalar X: 

short a + 2*T + (n-ro)* (b+T) 

long c + 2*T + (n-m)* (d+T) 

Vector X f (q - p 2: n - m - 1) : 

short e + T + (n-m)* (f+2*T) 

long g + T + (n-m) * (h+2*T) 

Vector X, (q - p < n - it - 1) : 

short i + 2*T + (n-m)* (b+T) + 
(q-p+l)*(j+T) 

long k + 2*T + (n-m)*(d+T) + 
(q-p+l)*(l+T) 



1 


1 


30 | 


40 


1 


50 


1 


65 


1 


75 


1 


k~ 


--+-- 


+- 




-+- 




-+- 




-+- 




h 



a 


2297 


834 


338 


j 99.0 


67.1 | 


b 


904 


387 


100 


| 26.0 


17.1 | 


c 


2408 


870 


351 


| 102 


69.2 | 


d 


1706 


484 


124 


| 30.4 


20.2 | 


e 


2629 


910 


370 


| 97.1 


64.6 | 


f 


1480 


459 


155 


| 41.2 


26.3 | 


q 


2740 


946 


383 


| 99.4 


66.7 | 


h 


3044 


844 


200 


| 50.0 


32.6 | 


i 


3370 


1197 


491 


j 140 


94.0 | 


i 


258 


91 


35 


| 9.5 


6.0 j 


k 


3581 


1233 


504 


| 143 


96.1 | 


1 


280 


101 


39 


| 10.7 


7.1 | 















IHEYGX 

Scalar X: 

short a + 4*T + (n-m)* (b+2*T) 

long c + 4*T + (n-m)* (d+2*T) 

Vector X„ (q - p £ n - m - 1) : 

short e + 2*1 + (n-m)* (£+4*T) 

long g + 2*T + (n-m)* (h+4*T) 

Vector X, (q-p<n-ro-l): 

short i + 4*T + (n-m)* (b+2*T) + 
(q-p+l)*( j+2*T) 

long k + 4*T + (n-m)* (d+2*T) + 
(q-p+l)*(l+2*T) 



30 



40 



50 



65 



75 



a 


3245 


1174 


478 




142 


1 95.4 | 


b 


2345 


664 


220 




58.3 


| 33.5 | 


c 


3471 


1221 


496 




145 


1 97.6 | 


d 


5519 


1433 


301 




73.2 


| 42.6 | 


e 


3187 


1114 


447 




127 


| 83.5 | 


f 


4368 


1226 


399 




103 


1 58.1 | 


q 


3533 


1161 


465 




129 


| 85.7 | 


h 


10636 


2746 


556 




133 


1 76.4 | 


i 


4087 


1459 


592 




171 


|114 | 


i 


545 


187 


72. 


8 


21.1 


| 12.7 | 


k 


4243 


1506 


610 




176 


|116 | 


1 


567 


196 


76. 


8 


22.3 


| 13.6 | 









IHEYGS, IHEYGL, IHEYGW, IHEYGZ 
Scalar X: a + (n-m) *b 
Vector X: 

(q-p2:n-m-l) : c + (n-m)*d 

(q-p<n-m-l) : e + (n-m)*b + (q-p+l)*f 
IHEYGS 



30 



40 



50 



65 



75 | 
.j 



1232 


4 30 


182 


49.5 


33.3 | 


461 


121 


37.6 


9.8 


5.3 | 


1871 


623 


259 


69.6 


45.4 | 


490 


128 


40.9 


10.5 


5.7 | 


2140 


733 


304 


83.6 


54.7 | 


29 


7.5 


3.3 


0.7 


0.4 | 











IHEYGL 










1 
X- 


30 


1 
-+- 


40 


1 
-_X_ 


50 


1 


65 


1 


75 



1320 


4 51 


189 


49.6 


33.3 | 


1241 


308 


56.9 


13.1 


7.3 | 


1959 


644 


266 


69.7 


45.4 | 


1270 


316 


60.2 


13.7 


7.7 | 


2228 


755 


311 


83.6 


54.7 | 


29 


7.5 


3.3 


0.7 


0.4 | 



l . 










IHEYGW 










1 


1 


30 


1 
-+- 


40 


1 
— +- 


50 


— T" 

1 


65 


— T 
1 


75 



a | 1396 


475 




| 198 




54.9 


36.1 | 


b | 1672 


425 




| 126 




30.5 


15.0 | 


c | 2035 


667 




| 275 




75.0 


48.1 | 


d | 1701 


432 




| 129 




31.2 


15.4 | 


e | 2304 


775 




| 320 




88.9 


57.4 | 


f | 29 


7. 


5 


1 3. 


3 1 


0.7 


0.4 | 
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IHEYGZ 



30 



40 



50 



65 



75 | 



a 


1572 


517 




| 211 




54.9 


36.1 | 


fc 


4824 


1184 




| 203 




44.2 


23.0 | 


c 


2211 


710 




| 288 




75.0 


48.1 | 


d 


4853 


1192 




| 207 




44.9 


23.4 | 


e 


2480 


821 




| 333 




89.0 


57.4 I 


f 


29 


7. 


5 


1 3 « 


3 | 


0.7 


0.4 1 
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INDEX 



ABS 

complex fixed-point 34-35 

complex floating-point 35-36 
ADD 

complex arguments 31 

real arguments 29 
ALL 72-73 
'and' operator 8-9 
ANY 72-73 
array indexers 

interleaved arrays 71-72 

simple arrays 70-71 
assignment operations 

tit string 12-13 

character string 16-17 
ATAN 

complex arguments 67-68 

real arguments 49-51 
ATAND (real arguments) 49-51. 
ATANH 

complex arguments 67-68 

real arguments 54-56 



bit string operators 
BOOL 14 



8-10 



comparison operator 
tit string 

tyte-aligned 11 
general 11-12 
character string 15-16 

concatenation operator 
bit string 10-11 
character string 14-15 

COS 

complex arguments 62-65 
real arguments 44-46 

COSD (real arguments) 44-46 

COSH 

complex arguments 62-65 
real arguments 51-53 

DIVIDE (complex fixed-point) 

division operator 

complex fixed-point 24 
complex floating-point 26 



ERF (real arguments) 56-58 
ERFC (real arguments) 56-58 
EXP 

complex arguments 60-61 
real arguments 40-42 
exponentiation operator 
complex operations 

floating-point exponents 
integer exponents 26-28 
real operations 

floating-point exponents 
integer exponents 20-22 



32-34 



28-29 



22-23 



fill operations 
bit string 12-13 
character string 16-17 



HIGH 16-17 



IHEABU 

see ; ABS (complex fixed-point) 
IHEABV 

see ; ABS (complex fixed-point) 
I HE ABM 

see ; ABS (complex floating-point) 
IHEABZ 

see ; ABS (complex floating-point) 
IHEADD 

see ; ADD (real arguments) 
IHEADV 

see ; ADD (complex arguments) 
IHEAPD 

see ; shift-and-assign, shift-and-load 
(real operations) 
IHEATL 

see ; ATAN (real arguments); ATAND (real 
arguments) 
IHEATS 

see ; ATAN (real arguments); ATAND (real 
arguments) 
IHEATW 

see ; ATAN (complex arguments) ; ATANH 
(complex arguments) 
IHEATZ 

see ; ATAN (complex arguments); ATANH 
(complex arguments) 
IHEBSA 

see ; 'and* operator 
IHEBSC 

see ; comparison operator (bit string, 
byte-aligned) 
IHEBSD 

see ; comparison operator (bit string, 
general) 
IHEBSF 

see ; BOOL 
IHEBSI 

see ; INDEX (bit string) 
IHEBSK 

see ; concatenation operator (bit string); 
REPEAT (bit string) 
IHEBSM 

see ; assignment operations (bit string) ; 
fill operations (bit string) 
IHEBSN 

see ; 'not' operator 
IHEBSO 

see ; 'or' operator 
IHEBSS 

see ; SUBSTR (bit string) 
IHECSC 

see ; comparison operator (character 
string) 



Index 
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IHECSI 

see ; INDEX (character string) 
IHECSK 

see ; concatenation operator (character 
string); REPEAT (character string) 
IHECSK 

see ; assignment operations (character 
string) ; fill operations (character 
string); HIGH; LOW 
IHECSS 

see ; SUBSTR (character string) 
IHEDVU 

see ; DIVIDE (complex fixed-point) 
IHEDVV 

see ; DIVIDE (complex fixed- point) 
IHEDZW 

see ; division operator (complex 
floating-point) 
IHEDZZ 

see ; division operator (complex 
floating-point) 
IHEEFI 

see ; ERF (real arguments) ; ERFC (real 
arguments) 
IHEEFS 

see ; ERF (real arguments); ERFC (real 
arguments) 
IHEEXL 

see ; EXP (real arguments) 
IHEEXS 

see ; EXP (real arguments) 
IHEEXW 

see ; EXP (complex arguments) 
IHEEXZ 

see ; EXP (complex arguments) 
IHEHTL 

see ; ATANH (real arguments) 
IHEHTS 

see ; ATANH (real arguments) 
IHEJXI 

see ; array indexers (interleaved arrays) 
IHEJXS 

see ; array indexers (simple arrays) 
IHELNL 

see ; LOG (real arguments); LOG2; LOG10 
IHELNS 

see ; LOG (real arguments); L0G2; LOG10 
IHELNW 

see ; LOG (complex arguments) 
IHELNZ 

see ; LOG (complex arguments) 
IHEMPU 

see ; MULTIPLY (complex fixed-point) 
IHEMFV 

see ; MULTIPLY (complex fixed-point) 
IHEMXB 

see ; MAX (real arguments); MIN (real 
arguments) 
IHEMXD 

see ; MAX (real arguments) ; MIN (real 
arguments) 
IHEMXL 

see ; MAX (real arguments) ; MIN (real 
arguments) 
IHEMXS 

see ; MAX (real arguments); MIN (real 
arguments) 



IHEMZU 

see ; multiplication operator (complex 
fixed-point) ; division operator 
(complex fixed- point) 
IHEMZV 

see ; multiplication operator (complex 
fixed-point) ; division operator 
(complex fixed-point) 
IHEMZW 

see ; multiplication operator (complex 
floating-point.) 
IHEMZZ 

see ; multiplication operator (complex 
floating-point ) 
IHENL1 

see; ALL; ANY 
IHENL2 

see ; ALL, ANY 
IHEPDF 

see ; PROD 
IHEPDL 

see ; PROD 
IHEPDS 

see ; PROD 
IHEPDW 

see : PROD 
IHEPDX 

see : PROD 
IHEPDZ 

see : PROD 
IHEPSF 

see : PROD 
IHEPSL 

see : PROD 
IHEPSS 

see : PROD 
IHEPSW 

see : PROD 
IHEPSX 

see : PROD 
IHEPSZ 

see : PROD 
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see ; SUM 
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see ; SUM 
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see ; SUM 
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IHETNI 

see : TAN (real arguments); TAND (real 
arguments) 
IHETNS 

see : TAN (real arguments) ; TAND (real 
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IHETNW 
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see : exponentiation operator (real opera- 
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